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1.0 


SUMMARY 


A  detailed  analysis  on  the  feasibility  of  efficient  CO  laser  frequency 
tripling  in  molecular  gases  has  been  carried  out.  The  results  of  the  study 
indicate  that  direct  frequency  tripling  of  multiline  CO  laser  output  leads  to 
small  conversion  efficiencies  based  on  the  nonlinear  susceptibilities  and  the 
efficiency  limiting  processes  of  the  molecules  surveyed.  However,  the  analysis 
of  the  three  resonance  enhancement  schemes  shows  that  two-photon  resonant  SFG 
and  THG  in  hydrogen  molecules  are  potentially  efficient  frequency  conversion 
processes  if  the  pump  frequencies  can  satisfy  the  resonance  condition. 

In  the  type  I  resonance  enhancement  scheme,  pump  absorption  leading  to  dynamic 
phase  mismatches  is  likely  to  be  the  prima-ry  efficiency  limiting  process. 

Since  this  scheme  requires  near-resonant  conditions  at  one-,  two-,  and  three- 
photon  resonances,  the  pump  absorption  is  always  present  and  inherently 
restricts  the  application  of  the  type  I  scheme  to  short  pump  pulses. 

The  type  II  scheme,  which  is  based  on  two-photon  resonance  enhancement,  avoids 
the  dynamic  phase  mismatch  caused  by  pump  absorption.  However,  other  competing 
nonlinear  processes  become  significant  at  the  high  pump  intensities  required 
for  efficient  conversion.  Although  two-photon  absorption,  quadratic  Kerr 
effect,  and  stimulated  Raman  scattering  ultimately  limit  the  efficiency,  near 
unity  THG  and  SFG  conversion  efficiencies  are  predicted  under  optimum  condi¬ 
tions.  The  parameters  that  govern  these  conditions  are  the  two-photon  reso¬ 
nance  frequency  detuning  and  Kerr-induced  phase  shifts. 

(3) 

The  type  III  scheme  analysis  and  the  related  x  calculations  which  were 
performed  utilizing  three-photon  resonances  in  a  number  of  infrared-active 
molecules  show  that  relatively  small  conversion  efficiencies  can  be  expected 
for  this  resonance  enhancement  scheme.  The  principal  efficiency  limitation 
in  the  type  III  scheme  is  optical  gas  breakdown  at  high  pump  intensities.  The 
high  pump  intensities  are  necessary  to  compensate  for  relatively  small  non¬ 
linear  susceptibilities. 

Efficient  multiline  frequency  up-conversion  of  CO  laser  frequencies  requires 
resonantly  enhanced  THG  and/or  SFG  processes  for  a  number  of  frequency  combi¬ 
nations.  The  nonlinear  medium  must  also  allow  simultaneous  phase-matching  for 


different  combinations  of  the  pump  frequencies.  These  requirements  suggest 
that  the  type  II  scheme  is  more  suitable  for  multiline  pumping  than  the  two 
other  schemes.  When  the  two-photon  resonance  condition  is  satisfied,  the  SFG 
nonlinear  susceptibilities  for  various  combinations  of  three  pump  frequencies 
are  very  similar.  The  type  II  candidate  molecules  such  as  H2  and  HD  have  very 
small  dispersion  at  the  pump  and  the  generated  frequencies.  This  property  is  * 

desirable  for  the  simultaneous  phase-matching  of  the  various  pump  frequency 
combinations. 

Efficient  CO  laser  output  spans  a  frequency  range  between  1920  and  2030  cm-*. 

Analysis  of  two-photon  resonances  in  H2  and  HD  indicates  that  the  sum  of  two 
pump  frequencies  lies  in  between  the  strong  Q-branch  resonances  of  H2  and  HD. 

This  leads  to  relatively  small  nonlinear  susceptibilities  for  direct  THG  and 
SFG.  Further  investigation  of  competing  effects  shows  that  stimulated  Raman 
scattering  may  be  significant  under  certain  conditions.  In  particular,  stimu¬ 
lated  rotational  Raman  scattering  can  become  an  efficient  nonlinear  process  in 
H2.  The  Raman  shifted  frequencies  may  participate  further  in  other  nonlinear  * 

interactions  such  as  induced  two-photon  absorption  and  sum  frequency  generation. 

The  Raman-shifted  frequencies  in  H2  and  the  original  pump  frequencies  closely 
match  the  Q-branch  two-photon  resonances  in  HD.  This  suggests  that  a  two-step 
frequency  conversion  involving  an  initial  Raman  shifting  in  H2  followed  by  sum 
frequency  generation  in  HD  is  a  potentially  efficient  process  for  multiline 
pumping.  The  two-step  conversion  process  leads  to  output  wavelengths  in  the 
1.77  and  1.90  region. 

These  results  indicate  that  two-photon  resonantly  enhanced  sum  frequency  gener¬ 
ation  in  HD  is  the  most  promising  approach  to  up-convert  multiline  CO  laser 
frequencies.  In  order  to  verify  this  feasibility,  an  experimental  investiga¬ 
tion  of  SFG  in  HD  is  recommended  for  future  study.  The  investigation  should 
consist  of  setting  up  an  experimental  facility  to  generate  tunable  IR  fre¬ 
quencies  in  the  5-6  p.m  region  and  performing  parametric  SFG  experiments  using 
the  generated  frequencies. 
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2.0  INTRODUCTION 

Efficient  frequency  up-conversion  of  high  power  infrared  gas  lasers  is  of 
considerable  interest  due  to  the  advent  of  very  efficient  IR  lasers.  For 
example,  the  frequency  conversion  of  CO  laser  output  to  shorter  wavelengths 
would  provide  high  power  coherent  radiation  suitable  for  atmospheric  trans¬ 
mission  with  smaller  beam  control  and  pointing  optics.  Recent  experimental 
and  theoretical  studies  indicate  that  IR  laser  frequency  tripling  in  molecular 
media  is  indeed  possible.^  High  conversion  efficiencies  at  reasonable  pump 
intensities  are  predicted  using  resonance  enhancement  of  the  third  order  non¬ 
linear  susceptibility  due  to  vibrational  states?'-6 

The  resonance  enhancement  of  frequency  up-conversion  processes  is  a  widely 
used  technique  which  is  particularly  useful  for  infrared  laser  frequency  up- 
conversion  in  molecular  media.  The  nonlinear  medium  may  consist  of  infrared- 
active  or  Raman-active  molecular  species.  In  the  IR-active  molecules,  the 
fundamental  and  overtone  transitions  provide  one-,  two-,  and  three-photon 
resonances  in  addition  to  contributions  from  the  electronic  states.  The  Raman- 
active  molecules  have  energy  level  structures  suitable  for  two-photon  resonances. 
These  resonances  may  yield  a  very  significant  enhancement  of  the  nonlinear 
susceptibility  leading  to  potentially  high  conversion  efficiencies. 

Several  recent  investigations  have  been  devoted  to  frequency  tripling  of  the 
CO2  laser  radiation.  Although  earlier  results  of  these  investigations  were 
not  promising,  some  very  recent  progress  in  the  demonstration  of  CO2  laser 

frequency  tripling  indicates  that  the  experimental  conversion  efficiencies  are 

.  2  3  4 

approaching  the  theoretically  predicted  values.  ’  ’  Analysis  of  the  particular 

resonance  enhancement  schemes  used  in  the  COj  experiments  shows  that  a  similar 

technique  is  applicable  to  CO  laser  frequency  tripling  for  generating  high 

power  coherent  radiation  at  1.7  hjh.  Since  the  CO  laser  output  consists  of 

several  frequencies,  sum  frequency  generation  may  occur  simultaneously  with 

frequency  tripling. 

Because  of  the  complexity  created  by  the  multitude  of  CO  spectral  lines,  as 
many  schemes  as  possible  should  be  considered.  For  this  reason,  three  differ¬ 
ent  resonance  enhancement  schemes  have  been  analyzed  in  terms  of  third  harmonic 


3 


and  sum  frequency  generation  in  several  different  molecular  media.  The  basic 
concept  of  two  of  these  schemes  has  been  discussed  in  the  literature.1-6 
Another  scheme  conceived  at  the  Northrop  Research  and  Technology  Center  (NRTC) 
involves  a  different  type  of  resonance  enhancement  which  may  solve  some  of  the 
problems  inherent  in  the  other  schemes. 

The  details  of  the  three  schemes  and  the  analysis  and  calculations  of  third- 
order  susceptibilities  for  each  of  the  schemes  are  described  in  Section  3. 
Section  4  discusses  the  efficiency  limiting  processes  and  phase-matching  anal¬ 
ysis,  which  are  essential  in  understanding  the  optimum  device  performance 
characteristics  and  conversion  efficiencies.  Section  5  analyzes  the  effect  of 
a  multiline  pump  source  on  the  third  harmonic  and  sum  frequency  generation 
efficiencies  and  suggests  a  potentially  efficient  two-step  frequency  conver¬ 
sion  scheme  based  on  an  intermediate  frequency  shifting  technique.  Recommen¬ 
dations  for  an  experimental  program  aimed  at  verification  of  theoretical  cal¬ 
culations  are  presented  in  Section  6. 
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THIRD  HARMONIC  AND  SUM  FREQUENCY  GENERATION:  NONLINEAR  SUSCEPTIBILITY 
ANALYSIS  AND  CALCULATION 

Optical  third  harmonic  generation  (THG)  or  frequency  tripling  is  described  by 
a  microscopic  nonlinear  susceptibility  tensor  x^  which  relates  the  cubic 
power  of  the  incident  electric  field  amplitude  to  the  generated  polarization 
in  the  medium.  More  generally,  the  quantity  x^  relates  a  product  of  three 
electric  field  amplitudes  to  the  generated  polarization  oscillating  at  the  sum 
of  the  three  field  frequencies.  This  is  referred  to  as  sum  frequency  generation 
(SFG) . 


The  theory  of  optical  frequency  mixing  has  been  well  established  and  experi¬ 
mentally  verified.  The  basic  quantity  that  enters  into  the  equations  which 
describe  the  THG  and  SFG  processes  is  the  macroscopic  nonlinear  polarization 
of  the  medium.  For  the  case  of  three  monochromatic  (plane-wave)  applied  fields, 
the  generated  polarization  at  the  sum  frequency  can  be  expressed  as:^ 

P(3)(z,t)  =  \  Ps(z)  e  1  st  +  c.  c.  (1) 

where  P«.(z)  =  |  zxj3)  (wp,  V  Ep(z)  Eq(z)  Mz) 

x  exp  [ i ( kp  +  kq  +  kr)z] 

The  factor Z  is  the  number  of  distinguishable  permutations  of  the  applied  field 

amplitudes:  Z=  6  for  three  distinct  (including  sign)  frequencies,  Z  =  3  for  a 

pair  of  equal  frequencies,  and  Z  =  1  for  w  =  w  =  u;  .  The  other  quantities 
(3}  P  P  r 

N  and  xv  ;(w  >  ju  ,  ju  )  are  the  total  number  density  and  microscopic  (molecular) 

S  P  H  ' 

third-order  susceptibility. 


The  macroscopic  nonlinear  polarization  provides  the  source  term  for  the  wave 
equation  governing  the  growth  of  the  generated  field.  In  the  slowly-varying 
amplitude  approximation,  the  wave  equation  reduces  to: 


^(z) 

dZ 


i  T|  -i  k_z 

TTT1  ps^  e 


(2) 
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v^*--  -S.f- 


where 


the  refractive  index  at  u>_ 


ns  "s 


Tt  =  foi  /e  )*,  the  impedance  of  free  space 
Substitution  of  the  expression  for  Ps(z)  from  expression  (1)  into  (2)  yields 


dE$(z)  i  T,  ujcN 


bZ 


8  n 


—  Z  xP}  En(z)  E  (z)  E  (z)  e’lAkz 


s  p'  '  q'  '  r 


(3) 


where  the  wave  vector  mismatch,  Ak,  is  defined  by 

Ak  =  ks  -  (kp  +  kq  +  kr) 


(4) 


If  the  depletion  of  the  applied  fields  is  neglected,  equation  (3)  can  be 
integrated  directly  over  the  interaction  length,  L: 


ES(L) 


i  T]  uu  NL 


Tfr1  Z  43)  EP  Eq  Er  e'“kL/2  s1"c  (“kL/2)  (5) 


The  electric  field  intensities  are  given  by 


h 


I, 


"1  lE,l2 

2 r, 

nslEsl2 


(X  =  p,q,r) 


(6a) 

(6b) 


Equations  (5)  and  (6)  are  applicable  to  third  harmonic  generation  as  well  as 
to  sum  frequency  generation.  For  example,  the  third  harmonic  field  at  ■  2a>  ^ 
is  given  by 

E,  (L)  =  ^  '"l  LN  .THG  ,  3  .-iAkL/2  sine  (AkL/2)  (7) 

Jjil  Tr^"  x  L1  e 

Accordingly,  the  field  intensity  of  the  third  harmonic  is  given  by 

1 4  O  O  O  TUP  2 


I 


THG 


=  jg  — Ou^)2  L2N2  |xTHGl  Ij3  sine2  (AkL/2)  (8) 
n_  n,J 


*s  "1 


6 


where 


IyHg  =  the  field  intensity  at 


THG 


the  field  intensity  at 

X?)  (“>1  »  "Ui  »  <".) 

3  o»  ^  i  i  i 

the  nonlinear  susceptibility  for  THG. 


All  of  the  parameters  in  Equation  (8)  are  expressed  in  MKS  units.  In  sum  fre¬ 
quency  generation  ui  may  be  either  a  sum  of  three  different  frequencies 
(ul>s  *  Wj  +  u»2  +  ‘“3)  or  a  sum  of  one  frequency  and  twice  another  frequency 
(e.g.,  o)s  =  2o» ^  +  0^).  The  generated  field  intensities  for  these  two  cases 
are  given  by 


lSFG 


4  n  n.n0n, 
s  1  2  3 


2W2(  2 
iu  N  L 


|Xg3^  (UJi»a)2’Ui3' 


,2  o 

I  Ijlglj  sine 


UkL/2)  (9) 


where  u»  =  +  w ^  2nd 


rSFG'  = 


1 

IF 


nsnl2n2 


(“s  )  N2L2  IxJ3^(«1.*1,»2)|2  I^lf 
x  sine2  (AkL/2) 


(10) 


where  =  2«^  + 

Equations  (8),  (9),  and  (10)  indicate  that  the  generated  intensities  for  THG 

(31 

and  SFG  depend  on  the  molecular  parameters  through  x  and  Ak.  Thus, 

efficient  THG  and  SFG  from  CO  laser  output  requires  molecular  species  in  which 

large  values  of  x^  can  be  obtained  at  the  CO  laser  frequencies.  Furthermore, 

the  phase  mismatch,  AkL,  must  be  minimized  in  order  to  obtain  a  value  of 
2 

sine  (AkL/2)  which  is  as  close  to  unity  as  possible.  More  detailed  analyses 
of  these  important  considerations  are  presented  in  the  following  sections. 

3. 1  Third  Order  Susceptibility  x^ 

The  derivation  of  the  nonlinear  susceptibility  tensor  is  based  on  density  matrix 

perturbation  calculations^  for  atoms  and  molecules  interacting  with  external 

(31 

electric  fields.  The  formula  for  x  is  given  bv  the  set  of  equations  in 

(31 

Table  1.  The  application  of  the  x  formula  for  molecules  involves  summation 
over  intermediate  states,  as  shown  by  the  expressions  in  Table  2.  In  principle. 
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all  of  the  rotational,  vibrational,  and  electronic  states  of  the  molecule  are 
required  for  the  calculation.  However,  various  approximations  can  be  used  to 
simplify  the  summations  under  certain  conditions.  For  example,  if  the  applied 
field  frequencies  or  their  combinations  are  nearly  resonant  with  the  molecular 
transition  frequencies  of  the  nonlinear  medium,  a  large  enhancement  of  the  fre¬ 
quency  factor  results  from  the  intermediate  states  associated  with  the  resonances 
This  resonance  enhancement  not  only  gives  large  susceptibilities  but  also  reduces 
the  number  of  intermediate  states  which  are  necessary  for  accurate  computation. 


TABLE  1.  NONLINEAR  SUSCEPTIBILITY  x 


(3) 


Process : 


,(3) 


10  +  U)  +  <0  -•  uo 

p  q  r  s 


X'  '  (to  ,tu  ,to  )  * 
s  '  p  q  r; 


s&)£ 


paa  ^ab  ^bc  ^cd  ^da  Fabcd  ^l,u,2,a)3^ 


^  =  sum  over  all  permutations  of  for 

pff)  s  population  distribution  of  molecules 

aa 

3  electric  dipole  matrix  element 
Fabcd  ta'i*u,2*u,3^  s  frequency  factor 


(u,ab+  1  rab  +  "l  +  ^2  +  *3)(uJac  +  1  Fac+  *2  +  ^Kd  +  1  rad  +  a,3) 


(iuab  •  1  Fab-  “V^ac  +  1  rac  +  *2  +  *3)(*ad  +  1  rad  +  *3] 


(uJab  ’  1  rab  •  "l^ac  -  1  rac  ‘  *1  ‘  *2)(uJad  +  1  rad  +  ^ 


(*ab  -  1  rab"  "PKc  -  1  rac  *  *1  *  a,2)(a'ad  *  1  rad  "  '"l  ‘  “2  "  ^ 
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rt.tM’jMt*-  ■  Hw'-VW 


TABLE  2.  x(3)  FORMULA 


x(3>  . 

S(^)£  £»v£fj£»r  w  %bcd 

LY]  [v]  J  [K] 

*v 

hab  14 be  14  cd  **da  with  va  =  0 

II 

o 

H- 

Q'1  ( 2J  +  1)  e_E(J)/'cT  _  partit^on  faction) 

VR 

1  R  R  R  R 

2  j+1  ^  ^ab  ^bc  ^cd  ^da 

L  Kj 

LK]  - 

set  of  (Ja,  Jb,  Jc,  J^)  allowed  by  the  selection  rules 

with  J,  =  0 

a 

^abcd 

frequency  factor  determined  by  (u^,  oj2,  uu3)  and 
molecular  transition  frequencies 

2  Resonance 

Enhancement  Schemes 

The  resonance  enhancement  schemes  considered  for  CO  laser  frequency  tripling 
are  based  on  one,  two,  and/or  three  photon  resonances  in  infrared-acti ve  and 
Raman-active  diatomic  molecules.  Figure  1  illustrates  the  three  different 
schemes.  The  difference  between  the  transition  resonance  frequency  and  the 
pump  frequency  or  the  combination  of  the  pump  frequencies  is  the  frequency 
detuning  which  determines  the  resonance  enhancement  of  the  frequency  factor. 

The  type  I  scheme  utilizes  near-resonant  frequency  detunings  at  the  fundamental 
and  the  first  two  overtone  transitions  corresponding  to  the  vibrational  levels 
v  =  1,  2,  and  3.  The  candidate  molecules  are  IR  active  species  with  fundamental 
frequencies  near  2000  cm'1.  Table  3  lists  the  molecules  and  the  optimum  THG 
pump  frequency  ranges  for  two  and  three  photon  resonance  enhancement. 
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TABLE  3,  TYPE  I  SCHEME  OPTIMUM  PUMP  FREQUENCIES 


Molecule 

THg  Optimum  bump  Frequency  Ranqe  cnrll  f 

2  Photon  Resonant 

3-Photon  Resonant 

CO18 

2050  -  2065 

2055  -  2060 

Dcx 35, 37 

1990  -  2025 

2005  -  2020 

QBr79,8I 

1860  -  1890 

1805  -  1810 

NO 

1840  -  1890 

1830  -  1860 

The  type  II  scheme  is  based  on  two-photon  resonance  in  Raman-active  molecules. 
The  candidate  molecules  and  the  optimum  two-photon  resonant  pump  frequencies 
are  listed  in  Table  4. 


TABLE  4.  TYPE  II  SCHEME  MOLECULES 
Two-Photon  Resonance  Frequency 


Molecule 

Two-Photon 
Branch  (J) 

Transition 
u>  (cm"*) 

Optimum  Pump 
Frequency  (cm-1) 

H2 

Q(0) 

4161.1 

2080.55 

Q(l) 

4155.2 

2077.60 

C\J 

O' 

4143.3 

2071.65 

Q(3) 

4125.3 

2062.65 

O 

4101.2 

2050.60 

0(2) 

3806.8 

1903.40 

0(3) 

3568.4 

1784.20 

HD 

Q(0) 

3632.1 

1816.05 

Q(l) 

3628.2 

1814.10 

Q(2) 

3620.5 

1810.25 

Q(3) 

3608.9 

1804.45 

Q(4) 

3593.3 

1796.65 

S(0) 

3887.5 

1943.75 

S(l) 

4051.8 

2025.90 

HF 

Q(0) 

3961.54 

1980.77 

Q(D 

3960.00 

1980.00 

Q(2) 

3956.90 

1978.45 

Q(3) 

3952.26 

1976.13 

Q(4) 

3946.06 

1973.03 

Q(5) 

3938.30 

1969.15 

The  type  III  scheme  derives  three-photon  resonance  enhancement  from  small  fre¬ 
quency  detunings  at  the  v  =  2  excited  vibrational  state  of  the  molecule.  The 
candidate  molecules  are  OF  and  HG£,  both  of  which  have  first  overtone  transi¬ 
tions  in  the  5500  to  5800  cm’*  range. 

1 3\ 

3.3  Nonlinear  Susceptibility  xv  Code 

A  general  third  order  nonlinear  susceptibility  code  for  infrared-active  mole¬ 
cules  has  been  developed  using  vibration-rotation  transitions  in  the  ground 

(31 

electronic  state  of  the  molecule.  The  code  is  based  on  the  x  formula  given 
in  Table  2,  and  is  capable  of  calculating  nonlinear  susceptibilities  for  THG, 
SFG,  quadratic  Kerr  effect,  two-photon  absorption,  and  stimulated  Raman  effect 
for  the  type  I  and  type  III  scheme  molecules.  The  necessary  input  parameters 
are  the  equilibrium  molecular  energy  level  constants,  gas  temperature,  the 
number  of  vibrational  and  rotational  intermediate  states,  the  dipole  matrix 
elements,  and  the  three  pump  frequency  values.  For  the  type  III  calculation, 
the  code  computes  only  the  near-resonant  contributions  from  the  R-branch  over¬ 
tone  transitions.'  TheX^  code  consists  of  the  pump  and  the  generated  fre¬ 
quencies  as  well  as  the  real  part,  the  complex  part,  and  the  modulus  of  the 
nonlinear  susceptibility  tensor. 


(31 

The  x  calculation  for  the  type  II  resonance  enhancement  scheme  involves  vir¬ 
tual  vibronic  transitions  in  Raman-active  molecules.  According  to  the  literature, 

the  perturbation  calculation  requires  extensive  computation  involving  vioronic 

(31 

matrix  elements  which  are  not  generally  available.  However,  x  values  may  be 
obtained  from  differential  Raman  scattering  cross  section  data  and  appropriate 

c 

frequency  factor  scaling.  This  procedure  is  believed  to  be  a  reliable  method. 

(31 

The  details  of  the  x  program  structure  and  the  procedures  used  in  the  com¬ 
putation  are  described  in  the  following  sections. 

3.3.1  Type  I  Scheme 

The  calculation  of  nonlinear  susceptibilities  for  molecular  media  in  principle 
requires  a  knowledge  of  all  of  the  energy  levels  and  of  all  of  the  transition 
matrix  elements  for  the  molecules  of  the  nonlinear  medium.  When  resonance 
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enhancement  occurs,  however,  a  knowledge  of  only  a  small  group  of  intermediate 
states  may  be  sufficient  to  determine  good  approximate  values  of  The  accu¬ 
racy  of  this  approximation  depends  on  the  frequency  detunings.  Thus  a  reliable 
resonance  enhancement  calculation  requires  a  precise  knowledge  of  the  molecular 
energy  level  scheme.  The  molecular  constants  taken  from  spectroscopic  data  are 
generally  adequate  for  computing  frequency  detunings  to  within  approximately 
1  cm”  . 

The  equations  for  computing  the  molecular  energy  levels  are  summarized  below. 
Each  level  is  characterized  by  vibrational  and  rotational  quantum  numbers.  The 


total  energy,  T(v,J) 
the  rotational  term. 

is  expressed  as  a  sum  of  the  vibrational  term,  G(v), 

FV(J). 

and 

T(v,  J) 

= 

G(v)  +  Fy(J) 

(ID 

G(v) 

= 

d)e(v+j)  -  -ugxe(v+s)2  +  weye(v+s>)3  +  WgZe(v4t)4  +  ... 

(12) 

FV(J) 

= 

BVJ(J+1)  -  DvLJ(J+l)j2  +  HV[J(J+1)]3  +  ...  (A-0) 

(13) 

FV(J) 

s 

-  V(-l)1  V2+Bv.(J+t)2-  Dvi(J^)4+Hvi(J+i)6  +  •• 

(A?  0) 

(14) 

where 

Bv 

= 

Be  “  ae(v4t)  +  Vg(v-tf)2  +  ... 

(15) 

°v 

= 

Dg-  9e(v*i)  +  ...  ' 

(16) 

Hv 

H  +  ... 
e 

(17) 

Av 

* 

%  Bv 

(.18) 

J 

v 

3 

Ay  -  2By  -  2Bv2/Av  +  ... 

(19) 

Bvi 

S 

Bv  ♦  0V  +  (-1)1  Bv:%+  2Xy2  +  ...] 

(20) 

°vi 

* 

Dv  -  (-1)1  Bv[\y3  ♦  6Xy4  +  ...] 

(21) 

Hvi 

3 

Hv  +  (-1)1  BvL2av5  +  20\y6  +  ...] 

(22) 

The  transition  frequency  between  any  two  levels  can  be  calculated  using 
Equations  (15)  through  (22).  The  calculated  frequencies  can  then  be  used 

to 

compute  the  frequency  denominator  term  in  Table  1. 
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Direct  calculation  of  the  nonlinear  susceptibility  requires  a  variety  of  sum¬ 
mations  according  to  Table  2.  Since  the  molecular  states  are  labeled  by 
vibrational  (v)  and  rotational  (J)  quantum  numbers,  an  explicit  form  of  the 
expression  in  Table  2  requires  the  calculation  of 

#R  V  OTT  EEEE  [<Jama|  coselJ^HJ^IcoseU^) 

ma  mc  md  (23) 

< Jcmcl  cose  1  Jdmd>  <Jdmdl  cose | Jama> ] 

In  Equation  (23)  the  indices  b,  c,  and  d  denote  the  vibrational  or  vibronic 
quantum  numbers  of  the  intermediate  states.  The  variables  (J^,  Jc>  Jd)  are  the 
rotational  quantum  numbers  of  the  intermediate  states.  The  frequency  denominator 
term  Fa5C(j(Ja»  Jb»  Jc»  Jd),  indicates  the  dependence  of  on  vibrational  and 
rotational  energy  levels.  Equation  (23)  gives  the  product  of  the  rotational  matrix 
elements.  The  selection  rules  for  linearly  polarized  fields  are  AJ  =  ±1  and 
Am»  0  for  a  molecular  state  with  A  =  0,  and  AJ  =  0,  ±1  and  Am  =  0  for  a  molecular 
state  with  A  =  1.  Using  these  selection  rules,  ^lR(Ja»Jb»Jc»JCj)  can  be  calculated 
for  allowed  values  of  Jb>  Jc» The  results  for  the  A  =  0  case  are  summarized 
in  Table  5. 


TABLE  5.  PRODUCT  OF  ROTATIONAL  MATRIX  ELEMENTS 
FOR  A  =  0  MOLECULES 


Path 

Jb 

Jc 

Jd 

^a,Jb,Jc,JcP  Ja=J 

1 

J-l 

J-2 

J-l 

2J  (J-l) 

15  (2j-l)  (2J+1) 

2 

J-l 

J 

J-l 

J  (4J2  +  l)  . 

15  (23-1)  (2J+1)2 

3 

J+l 

J 

J-l 

2J( J+l) 

TTWiy2 

4 

J-l 

J 

J+l 

jj^j+ij. 

5 

J+l 

J 

J+l 

(J+l) (4J2+8J+5) _ 

15  (2J+3)(2J+1)2 

6 

J+l 

J+2 

J+l 

2( J+l) ( J+2) 
13"(2J+3)(2J+1) 

V 
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The  experimental  data  on  vibrational  matrix  elements  are  available  for  a  number 
of  infrared-active  diatomic  molecules.  In  most  cases  Av  *  ±1  transitions  are 
well  approximated  by  the  harmonic  oscillator  matrix  elements  which  are  given  by 

<v  lp.[  v-1)  *  yfiT  (1  {m*  1  0)  (24) 

The  matrix  elements  for  overtone  transitions  (Av  =  2  and  3)  may  be  expressed  as 

^v,0  ~  ^  0  ^*1  0 

where  Sy  g  is  the  strength  relative  to  the  fundamental  (0  -  1)  transition. 
Typically,  S2  Q  «  10"1  and  S3  Q  »  5  x  10“3. 

The  population  distribution  in  the  rotational  states  is  represented  by  a  factor 
fj.  For  a  gaseous  medium  in  thermal  equilibrium  the  distribution  is  considered 
Boltzmann.  The  quantity  QRot  is  the  rotational  partition  function  which  is 
defined  as 


QRot  =  2  (2J+1)  exp  (_hc  F0(J)/kT) 


For  he  B0«  kT,  QR  t  may  be  approximated  by 


The  summations  involving  vibrational  quantum  numbers  and  frequency  permutations 
which  are  needed  for  the  A  =  1  case  are  exactly  the  same  as  the  summations  that 
are  required  for  the  A  =  0  case.  However,  since  the  rotational  states  have  two 
sets  of  levels  (Q=§,  f ) ,  two  distinct  summations  over  J  must  be  performed. 


.(3)“EE/fe  <--*i (3-ii 

v  bed  Lj=i  j=a 


In  each  summation  over  0,  there  are  nineteen  different  paths  to  be  summed  as  a 
result  of  allowed  values  of  (J^,  Jc>  J^). 

2  (**  *  ^  2  ?(Ja»Jb,Jc,Jd^  Fabcd 

j=s  VcJd 
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In  each  path  of  rotational  quantum  numbers,  the  Q  values  of  the  intermediate 
states  can  be  either  £  or  j,  making  eight  different  contributions  to  each  set 
of 


)  1  0 ( J. » ^b »*^c ’^d^ 


V'W'M'  'abed  "  % 

abed  path  (J) 


,S  "<VbJcJd> 


(-1>  I"°  Z  tCa.b.c.d)  Fabcd 

4  d,  +  1  trs  \ 


a  1  path  (n) 

where  s  “  Ja  +  db  +  ^c  +  Jd‘  w^Ja,Jb,Jc,Jd^  is  defined  by; 


(30) 


«w«"V  -ZZZZ/**  1  Jb\/Jb  i  i 

ma  mc  ny  \-ma  0  0  mc/\rmc  0  m^A-m.4  0  m 


Using  the  properties  of  the  summations  involving  products  of  3j  symbols' 
expression  (31)  reduces  to 
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W(Ja.Jb’Jc+)d) 


=  L11 


(-Jb+2JC+Jd) 


(-1) 


Jb+Jd 


1  Ja  Jb 


1  °c  Jd 


(31) 


2 

k=0 


J  1  Jh 
c  b 


1  J  k 

a 


Jc  1  Jd 


1  J  k 

a 


h(k) 


(32) 


wtere  h(k) .  [3k2  (k+1)2  -  21  k  (k+1)  +  20] 


(33) 


By  specifying  a  particular  set  of  rotational  quantum  numbers,  expression  (32) 
can  be  evaluated  using  6  j  symbol  formulas,^  and  the  results  are  summarized 
in  Table  6. 


The  frequency  factor  is  computed  in  the  same  way  as  was  done  for  the  A  =  0 
case.  The  only  difference  is  the  formulae  used  for  calculating  molecular  levels 
with  quantum  numbers  {v,  J  ,  n}. 

The  summation  over  intermediate  a  states  can  be  written 

JLd  +(a»b,c,d)  FabC<j  ZZZ  S(a,b)  §(b,c)  §(c,d)  §{d.a)  Fabcd  (34) 
path  (w)  ^b  ^c  nd 
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The  terms  of  the  form  5 (a, b)  represent  effective  reduced  matrix  elements  between 

states  a  and  b  when  both  states  are  expressed  in  the  intermediate  coupling  scheme. 

In  this  coupling  scheme,  molecular  states  are  expressed  as  linear  combinations  of 

9 

unperturbed  eigenfunctions  with  coefficients  defined  by  the  following  relations: 


and 


\\> 

2  INT 


l  2ni> 


2  INT 

where  the  mixing  coefficients  are  given  by: 


CA(v)  l2n|>  -  CB(v)  |V> 

CA(v)  l2n  >  +  CB(v)  (2tti> 

a  2 


and 


CA(v) 

CB(  v) 
X 


Xv  -_2  *  \ 

2~r 


xv  +  2  •  Xv 


2Xv 

L^v^v-4^  + 


VBv 


(35) 

(36) 


(37) 


The  rotational  reduced  matrix  elements  are  also  necessary  to  compute  §.  For 


Hr  0  cases,  these  matrix  elements  are  given  by 
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Ra  ) 


K^r 

for  J'  =  J+l 

for  J'  =  J 

r,2  ~2~\% 

N  -  a 

for  J'  =  J-l 

L  J  J 

(38) 


The  three  different  forms  of  I  in  terms  of  the  expressions  given  in  Equations 
(35)  through  (38)  are  given  in  Table  7. 
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TABLE  7.  Sa(a,b)  IN  TERMS  OF  MIXING  COEFFICIENTS 
AND  REDUCED  ROTATIONAL  MATRIX  ELEMENTS 


CA(a)  CA(b)  Ri  (J  ,J.)  +  CB(a)  CB(b)  Rjl  (Ja,JK) 

2  a  D  z  fl  O 

CA(a)  CA(b)  R|.  (Ja>0b)+  CB(a)  CB(b)  R^ 

CA(a)  CB(b)  Ri  (J  ,JJ-  CA(b)  CB(a)  R^  (Ja,JK) 

2  a  D  z  «  D 


In  Table  8,  f(a,b,c,d)  is  defined  in  terms  of  appropriate  products  of  the  effective 
reduced  matrix  elements,  determined  by  the  possible  combinations  (paths)  of  the 
quantum  numbers,  a  ,  Q.  ,  a  ,  and  a  .. 

a  D  C  G 

TABLE  8.  TERMS  IN  THE  .1-PATH  SUMMATION 


1 

2 

2 

2 

2 

l 

1 

1 

1 

2 

2 

I 

2 

■» 

1 
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The  expressions  presented  in  this  section  constitute  a  general  method  for  calcu 

lating  third-order  susceptibilities  in  IR-active  molecules.  Figure  2  shows  a 

(31 

flow  chart  of  the  x  code  based  on  these  expressions.  A  listing  of  the  com¬ 
puter  program  is  given  in  Appendix  A. 
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3.3.2  Type  II  Scheme 

(3} 

In  the  type  II  scheme,  analysis  of  the  x  formula  shows  that  THG  and  SFG, 

stimulated  Raman  scattering  (SRS),  two-photon  absorption  (TPA) ,  and  quadratic 

Kerr  effect  share  the  same  two-photon  resonance  and  vibronic  intermediate  states. 

Consequently,  the  susceptibilities  are  related  and  may  be  evaluated  from  meas- 

SRS 

ured  data  on  one  of  the  processes.  A  direct  relation  between  x  and  differ¬ 
ential  Raman  scattering  cross  section,  da  /da  ,  enables  the  calculation  of  the 

various  susceptibilities.  The  relation  has  been  derived  using  the  definitions 
'  12  13 

of  the  scattering  cross  sections.  ’  A  derivation  consistent  with  the  definition 
f  31 

of  x  of  this  report  is  outlined  below. 


The  starting  point  is  Placzek’s  definition  of  the  "effective  cross  section"  for 

scattered  light  due  to  a  transition  k  -  n  in  both  scattering  polarization 
12 

directions: 


Qkn  =  Z! 


3 


pa ;knl 


(39) 


where 


and 


[Placzek's  Equation  (5.5a)] 


<U„'H  £ 


(Mo  )kr  (lip) 


rn 


(Mp) 


kr  (Mo  >rn 


UJ  ,  - 

rk 


U)  +11) 

rn 


[Placzek's  Equation  (5.2)] 


(40) 


a  =  1,  2,  or  3  corresponding  to  the  x,  y,  or  z  directions  of 

polarization  for  the  incident  radiation,  respectively 

P  =  1,  2,  or  3  corresponding  to  the  x,  y,  or  z  directions  of 

polarization  for  the  scattered  radiation,  respectively 

=  the  wavelength  of  the  scattered  radiation 

<s>kr-  /  ’*  K  h  * 

A 

M  =  the  electric  dipole  operator  of  polarization  specified  by  y. 

*  i 

| 

i 
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In  the  equation  for  xv  »  the  states  are  labeled  as  (a,b,c,d),  so  that  k=|a,Ja,m> 
and  n=|c,Jcmc>.  Consider  the  case  of  parallel  polarization  (p  =  a) 


27  n5 

^a,Ja,ma  -c,Jc,mc)||  *  ,  x  4  I  (Cacj)a,J!i,nia;cfJ,.,mJ 


'll  3  A 


a’  a’  c’  c 


Summing  over  the  degenerate  m  states  and  dividing  by  the  degeneracy  of  the  initial 
state  (according  to  Equation  (5.5b)  of  Placzek)  yields 

(Qa,Ja-  c,Jc),|  =  (^r)  l(Caa)a,Ja,ma;c,Jc,mc|2  (42) 


ma  mc 


where  k  =  the  wavevector  for  the  scattered  radiation  (=  2n  /\  ). 

The  differential  Raman  scattering  cross  section  (da/d&  )  at  90°  to  the  direction 
of  the  incident  radiation  and  at  90°  to  the  polarization  direction  of  the  inci¬ 
dent  radiation  has  the  maximum  value  (da/drt)  .  This  is  the  value  measured 

iTlaX 

experimentally.  The  relationship  between  die  total  cross  section,  (CT)tota-|> 
and  the  differential  cross  section,  (da /dQ  )  a  .  for  the  case  of  isotropic 

■t  o  iTlaX 

scattering  is. 

(a)total  =  H  (^)max  (43) 

where  (CT)tota^  is  the  (Qa  j  _  c  j  )  of  Equation  (42).  Hence,  the  polarized 
(da/d.'i)  is  given  by:  a  c  II 


da  \  s 


(to) 


ks  9a  ^Caa^a,J  ,m  ;c,J  ,m 

d  m  m  a  d  C  C 

ma  mc 


where  the  subscript  s  on  (da/dn  )s  indicates  the  cross  section  for  the  scat¬ 
tered  (Stokes)  radiation. 

(3) 

Analysis  of  the  form  of  xv  which  is  appropriate  for  Raman  scattering  shows 
that  for  a,J.  -*  c,J„  resonance: 

d  C 


*(3)<V~p-V  ■  ITT  E  E 


a,c  ma,mc 


2  ""ah  ^bc  Aab 
b 
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where 


ab  ui  .  + 

ab  p 


“ab  '  “s 


(46) 


and  F  is  the  half-width  at  half-maximum  corresponding  to  the  a,J.  -»  c,J 

a  c 

resonance.  Note  that 

ab  ba 


and  =  hi.  +  uu  ,  or  uu.  =  uu,  -  uun 

be  be  ac’  be  ba  Raman 


so  that 


-  JU  .  =  vU.  +  dtf_ 

ab  be  Raman 


and 


UU  .  -  U)  =  -  UU.  -  (JU  . 

ab  s  be  p 


Rewriting  Equation  (46)  using  the  last  relation 
A 


ab 


i  +  1 


“ba  *  *p  "be  +  “p 


Therefore, 

V  "P 


x(3)(“p,^p,“s)  «  zsr 

(0) 


a,c  mfl  mc 


paa  ^  ^Ccj^a,Ja,ma;e,Jc,mc' 


The  expression  for  paa  x  '  is  independent  of  rn  and  is  given  by 

aa  a 


»„<0)  ■  q-  exp  c-  Ea/kT] 


•  t  ^ 


(47) 


(48) 


where  Q  is  the  rotational  partition  function  for  molecules  of  the  Raman  medium. 
Thus, 


(3) 


(VV^s*  =  6Hr£  fa  1(Cca} 


ma  mc 


a , J  ,m  >c , J  ,m 

a  a  C  C 


=  “1  rf  _L  /da\ 
SET  Z*  a  ^4  \  d  1 


(4S) 
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w 


For  vibrational  Raman  scattering,  the  final  state  is  essentially  empty.  Thus 
the  summations  over  v,  and  v.  reduce  to  a  particular  initial  state  (v*0,J  *J) 

a  C  a  a 


x 


(3) 


S 


) 


-i 

6nr 


(50) 


Since  xSRS  =  Im  [xv3^ (w  .-w-.^  )  ] 

P  P  s  ’ 


1 


(51) 


The  quantities  in  expression  (51)  are  defined  in  CGS  units.  By  multiplying 

9  SRS 

(51)  by  (4neo)  and  expressing  all  quantities  in  MKS  units,  x  is  obtained 

also  in  MKS  units. 


In  order  to  relate  xSRS  of  Equation  (51)  to  other  third-order  susceptibilities, 
the  x{3)  formula  given  in  Table  1  has  been  analyzed.  The  results  are  summarized 
in  Table  9.  Since  the  cross  section  is  proportional  to  the  fourth  power  of  the 
Stokes  wave  vector,  the  peak  values  of  x^  (i.e.  Auj=0)  are  nearly  independent 
of  frequency  for  radiation  in  the  IR  and  visible  regions.  Thus,  Raman  scatter¬ 
ing  data  obtained  at  a  certain  visible  wavelength  can  be  used  to  compute  the 
peak  values  of  x^  for  all  of  the  third-order  processes  associated  with  the 
type  II  scheme. 


3.3.3  Type  III  Scheme 

The  type  III  resonance  enhancement  scheme  utilizes  three-photon  resonance  via 

vibrational-rotational  and  vibronic  intermediate  states.  A  Detailed  study  of 
(3) 

the  x  formula  indicates  that  an  exact  calculation  requires  vibronic  transi¬ 
tion  moments  and  Franck-Condon  amplitudes  for  each  of  the  excited  electronic 
states  coupled  to  the  ground  electronic  state  by  electric  dipole  transitions. 

An  extensive  literature  search  provided  only  a  small  fraction  of  the  necessary 
parameters  for  HC£  and  DF.  Because  of  this  lack  of  data,  various  approximations 
were  considered  in  order  to  derive  formulas  that  can  be  evaluated  numerically 
with  the  existing  data. 
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TABLE  9.  TYPE  II  NONLINEAR  SUSCEPTIBILITY 


'  ) 
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Two-Photon  Resonance  Linewidth  (HWHM) 


Polarized  Differential  Raman  Scattering  Cross  Section 


Raman  Stokes  Radiation  Wave  Vector 


Fractional  Population  Distribution 
Raman  Frequency 
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(2) 

The  calculation  of  x  for  candidate  molecules  in  the  type  III  resonance  en¬ 
hancement  scheme  involves  three  different  kinds  of  pairs  of  the  intermediate 
states  |c>  and  (d>:  The  two  intermediate  states  may  both  be  pure  vibrational- 
rotational;  one  intermediate  state  may  be  pure  vibrational-rotational  and  the 
other  intermediate  state  may  be  vibronic,  or  both  intermediate  states  may  be 
vibronic.  Pairs  of*  intermediate  states  of  the  first  kind  have  already  been 
encountered  in  the  calculation  of  x  for  candidate  molecules  in  the  type  I 
resonance  enhancement  scheme.  Accordingly,  the  treatment  of  the  contributions 
from  pairs  of  intermediate  state  of  the  first  kind  has  been  discussed  during 
the  analysis  of  the  type  I  resonance  enhancement  scheme  and  will  nolj  be  repeated 
here.  Parts  of  the  intermediate  states  of  the  third  kind  lead  to  terms  which 
are  small  because  of  the  square  of  the  electronic  energy  occurring  in  the 
denominator  of  these  terms.  Thus,  pairs  of  intermediate  states  of  the  second 
kind  will  be  considered  in  detail  in  the  following  discussion. 


The  third-order  susceptibility  for  the  type  III  resonance  enhancement  may  be 
approximated  as 


(3)  ^  (0)5p  ^ab^bc  1 _ ^cd  ^da 

“I  h3^  “  "ab  +  “3  *  lTab  ”ac  +  “2  *  ”3  ~  “ad  *  “3 

a  be  d 


(52) 


Consider  the  special  case  of  three-photon  resonance  at  va=0,  Jg=J  and  vb=2, 
Jb=J+l  (0-2  R(J)  transition) 


v  (3)  2:Z1C  (0)  ^02  _  ^2c^cd*d0 

Hl.res  ^haa  *02  +  “ 3  +  ir02  “T  (u,0c  +tU2  +  u,3} (iu0d  +  ^ 


The  summations  over  the  Intermediate  states  |c>  and  |d>  must  now  be  examined 
in  the  case  where  one  of  the  two  intermediate  states  is  a  vibronic  state.  Let 
"e"  stand  for  a  set  of  vibronic  levels  and  "I"  stand  for  the  summations  over 
the  intermediate  states. 


I 


ry  *2e  ^ed  ^dO  ,rr  **2c  ^ce  ^eO  T  +T 

d  e(VV¥V¥  c  e  ^0c+u,2  +  ttt3Ktu0e  +  Ui3r  1  2 


(54) 
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where  |uuQe|is  the  electronic  transition  frequency,  assumed  to  be  much  greater 
than  any  of' the  input  frequencies.  Analyzing  the  two  terms  in  Equation  (54) 
separately, 


!  _  V"*l  **2®**etf  +  ^2e^el'  0 

1  ‘^|(tu0e+  iu2+u,3)(u,00'+uj3)  ^Oe+V  'B3,('"oi#+J'3J 

jjgeJWJVQ  i 

(u,0e+iu2  +  u,3)^02'  +  tu3)  ) 


(55) 


,  .y!  ^20/  ^cfe  ^eO  +  ^ 2 1' ^ l'  e  ^eO 

2  e|(tu00'+tu2+  "3)(tu0e+  “3J  (a,0l'  +  "2  +  tu3)(u)0e  +  u,3) 

+  ^  .^e^eQ  | 

("0*  +U)2+  "3)("0e  +  "3J  | 

where  the  primes  indicate  summation  over  rotational  states. 


Note  that  |w0e|  :?>  w.  This  fact  leads  to  the  approximations  ju^e+  w2+  ujj-  j>Qe 
and  u)Qe  +  iu3  a  u)0e  which,  when  applied  to  Equations  (54)-(56j  yields 
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(57) 


Next,  the  summation  over  vibronic  intermediate  states  is  carried  out  assuming 
that  the  quantity  o>ge  may  be  energy  averaged  over  the  vibrational  levels  of  the 
electronic  state  to  yield  the  quantity  r.g  which  is  independent  of  the  vibrational 
quantum  numbers,  vg. 

e  e  J  v  ( 


(58) 


Making  this  second  approximation, 
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where  Jg  and  vg  are  the  rotational  and  vibrational  quantum  numbers  of  states 
in  the  excited  electronic  level  manifold.  The  above  approximation  allows  a 
simplification  of  the  summation  over  vibrational  levels  of  the  electronic  state. 
Note  that  the  products  of  the  matrix  elements  in  Equation  (58)  have  the  general 
form  w«2e  ^ev'  or  ^2/  ^v'e  ^eO*  Thus»  separating  the  rotational  and 
vibrational  parts  of  each  matrix  element,  these  two  products  can  be  explicitly 
written  as 

^2e  y,ev/  ^Vo  *  <9.2|n.| e,ve><e,vejn|  g,v'><g,v'|ii|g,0>  R^  (59) 

^v'  %'e  ^eO  =  g.v'Xg.v'l^j e ,ve> <e ,veln.l g,°>  R2  (60) 

where  Rj  and  R2  represent  the  appropriate  product  of  rotational  integrals.  Both 
terms  are  of  the  form 

<9.v*  ln|g,v,><g,v/|n|e.v^<e.ve|jk|gfv>  R  (61) 

The  summation  over  |e,vg>  can  thus  be  evaluated  using  the  completeness  property 
of  the  vibrational  states  in  the  manifold  of  the  excited  electronic  state. 


The  product  of  the  second  and  third  factors  in  Equation  (61)  is  performed  first. 

The  electronic  part  of  the  matrix  element  must  be  integrated.  It  is  assumed 

14 

that  the  individual  matrix  elements  can  be  written  as 
<9.v'  (nje,ve>  =  <v'UgeIve> 

(62) 

<e,vg  |n|g,v>  «  <veUgelv> 

where  yge  *  u-eg  is  the  dipole  moment  operator  for  the  nuclear  coordinates.  The 

nuclear  dipole  moment  operator  can  be  expressed  as  a  Taylor  series  expansion 

1 4 

about  the  equilibrium  internuclear  separation,  re: 

dn  a2 

^ge(r)  *  ^ge<re)'  +  (r're)  +  \  ~ (r-re)2+...  (63) 

re  ar  re 

Using  the  completeness  property  yields 


E  <»'  UgJVj)  <»,>„!  *Mv'  l,-9e,e9|v) 
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Noting  that  Uge|is  real,  the  operator  in  Equation  (64)  may  be  reexpressed  as 

Ip.  j  ^  =  h-  x  +  x3  +  0^  X3  +  ...  (65' 

1Pge‘  *ge  ^eg  ge  ge  ge  ge  v 

where  x  is  the  displacement  of  the  atoms  from  their  equilibrium  positions  (i.e. 
x  =  r-rg)  and  the  0^  are  all  real  coefficients.  These  coefficients  are  re¬ 
lated  to  certain  derivatives  of  (63)  according  to: 

’  Ke’l2  <66: 

=  Mge’  "g”  *  Mge’  ^ 


M*?'  +I«y*|2  +  H*°)  »j'2j 


ge  ge 


ge  ge 


where 


The  remaining  integral  in  Equations  (64)  and  (65)  is  of  the  form 

A^}v  S  V  |  xk|  v>  ( 7C 

where  |v')and  \  v)  are  assumed  to  be  normalized  eigenfunctions  of  an  anharmonic 
oscillator.  Thus,  Equations  (59)  and  (60)  can  be  evaluated  according  to: 

2  2 

£  ^2e  ^ev'  Vo  =  2  A^?  /  £  A R,  (71 


j  r\r%  • 

ge  2,v 


gg  v  o  ,xl 


2-4  ^2v'  Ve  Vo  2-4  'gg  A2 ,v  L  Dge  Vo  R2 


Since  the  frequency  denominators  do  not  depend  on  Jg,  the  summations  over 
rotational  states  are  carried  out  first.  Summation  over  J'  must  be  performed 
explicitly  using  actual  energy  levels.  Summation  over  m  states  can  be  per¬ 
formed  by  factoring  the  rotational  part  of  the  ^g2  matrix  element  into  Rj  and 
R2.  This  procedure  yields  precisely  the  function  pR  defined  by  Equation  (23). 


3.4  Calculation  Results 

(3) 

The  x  calculations  for  the  three  resonance  enhancement  schemes  have  been  per¬ 
formed  using  the  line  selected  CO  laser  frequencies  corresponding  to  P(9),  P ( 10) , 
P(  11) ,  and  P(  12)  transitions  in  the  [6  -  5],  [5  -  4],  and  [4-3]  vibrational 
manifold.  In  addition,  survey  calculations  have  been  carried  out  for  each  molec¬ 
ular  species  in  order  to  determine  the  pump  frequencies  for  which  the  nonlinear 
susceptibility  is  optimized.  Nonlinear  susceptibility  values  have  been  computed 
in  CGS  units;  conversion  to  MKS  units  is  obtained  by  multiplying  the  CGS  value 
by  a  factor  of  1.235  x  10’25. 

The  following  sections  summarize  the  results  of  the  calculations  for  each  reso¬ 
nance  enhancement  scheme. 

3.4.1  Type  I  Scheme 

18  15  37 

The  candidate  molecules  studied  for  the  type  I  scheme  are  CO  ,  NO,  DU  *  , 

79  81 

and  DBr  '  .  The  spectroscopic  constants,  transition  dipole  moments,  and  line- 

width  data  references  for  each  molecule  are  listed  in  Table  10. 


TABLE.  10.  TYPE  I  SCHEME  MOLECULAR  PARAMETER  REFERENCES 


Molecule 

Spectroscopic 

Constants 

Transition 
Dipole  Moments 

Linewidth 

CO18 

[15] 

[16] 

[17] 

NO 

[18] 

[19] 

[20] 

oa35 

£21] 

[22] 

[23] 

[24] 

[25] 

[26] 

37  35 

The  spectroscopic  constants  for  DG&  were  calculated  from  the  DCt  parameters 

27 

using  the  isotope  effect  relations.  The  same  procedure  was  used  to  obtain  the 
81  79 

constants  for  DBr  and  DBr  .  The  transition  moments  and  linewidths  for  the 
isotopes  were  assumed  to  be  the  same. 

The  results  of  the  x^  code  calculations  using  the  line  selected  CO  laser  fre- 

THG 

quencies  are  listed  in  Tables  11  through  16.  The  tables  indicate  that  x 

and  xSFG  values  are  in  the  range  of  10~38  to  10~38  cm8/erg  for  CO*®  and  DCX38’^, 

79  81 

while  the  nonlinear  susceptibilities  of  NO  and  DBr  ’  are  on  the  order  of 
10’3®  cm6/erg  and  10~40  cm®/erg,  respectively. 
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rPE  I  xv 1  CALCULATION  RESULTS 
MOLECULE:  f  CO18  ] 


Resonance;  (2)  -  Two-Photon  Resonance;  (3)  =  Three-Photon  Resonance 


1PE  I  CALCULATION  RESULTS 
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Resonance;  (2)  =  Two-Photon  Resonance;  (3)  =  Three-Photon  Resonance 


TABLE  15.  TYPE 


One- Photon  Resonance;  ( 2 ) 


rug  -34  fi 

A  conversion  efficiency  estimate  based  on  a  x  value  of  10  cm  /erg  with 

N  *  10^  cm-3,  I  =  500  MW/cm^,  and  L  *  300  cm  shows  that  nearly  20%  efficiency 

is  feasible  if  phase-matching  can  be  maintained  over  the  interaction  length 

f  31 

during  the  entire  pulse  length.  From  this  analysis,  x  values  on  the  order 

oc  c  _  c 

of  10  cm  /erg  suggest  efficiencies  on  the  order  of  2  x  10  .  The  small 

conversion  efficiencies  are  the  direct  result  of  the  small  values  of  the  non¬ 
linear  susceptibilities. 

TUp 

The  survey  calculations  of  x  were  carried  out  to  find  the  largest  two-  and 

three-photon  resonance  enhancement  in  each  of  the  type  I  scheme  molecules. 

Using  the  optimum  pump  frequency  ranges  given  in  Table  3  as  a  guide,  several 

28 

known  CO  laser  frequencies  were  chosen  to  closely  match  two-  or  three-photon 

resonances  in  DC4,DBr,  and  NO.  The  CO  laser  frequencies  did  not  match  the 

10 

optimum  pump  frequency  ranges  of  CO  .  Table  17  lists  the  results  of  the 

TUp 

survey  calculations.  The  x  values  in  the  table  show  significant  enhance¬ 
ment  for  NO  and  DBr.  This  suggests  that  the  choice  of  CO  laser  output  fre¬ 
quencies  to  match  the  optimum  pump  frequency  range  and  the  resonances  is  an 
important  factor  in  the  overall  efficiency. 


TABLE  17.  TYPE  I  THG  SURVEY  CALCULATION 


CO  Laser  Frequency 

u)  ( cnT  1 ) 

Transition 

Molecule 

x(3) 

1838.49 

P(24)  [9-8] 

NO 

1.41  x  10‘35 

1859.80 

P ( 19 )  [9-8] 

NO 

1.28  x  10'35 

1846.87 

P ( 16 )  [10-91 

NO 

1.45  x  10" 35 

2016.871 

P(12)  [4-3] 

DW35 

1.52  x  10-37 

2003.154 

P(9)  [5-4] 

DCS35 

0.82  x  10-37 

1991.014 

P ( 12 )  [5-4] 

Da37 

0.42  x  10’37 

1808.35 

P ( 1 3)  [12-11] 

DBr79 

1.1  x  10'37 

1809.41 

P  ( 19 )  [11-10] 

DBr79 

2.1  x  10'37 

1805.29 

P  ( 20 )  [11-10] 

DBr81 

2.1  x  10'36 

1809.10 

P  ( 25 )  [10-9] 

DBr81 

4.6  x  10”37 

38 


3.4.2  Type  II  Scheme 

The  references  used  for  the  molecular  parameters  are  listed  in  Table  18. 
Differential  Raman  cross  section  data  for  HD  are  not  available;  however,  the 
similarity  of  its  electronic  structure  to  that  of  H2  implies  a  similar  value 
of  dc/dQ  for  both  HD  and  H2.  Thus,  do/d a  values  for  HD  may  be  scaled  from 
H2  data  using  the  dependence  of  Raman  cross  section. 


TABLE  18.  TYPE  II  SCHEME  MOLECULAR  PARAMETER  REFERENCES 


Molecule 

Spectroscopic 

Constants 

(ro)„an<i  k 

Linewidth 

H2 

[29] 

[30] 

[31] 

HD 

[29] 

calculated 
from  [30] 

[32] 

HF 

[33] 

-p* 

i — * 

[35] 

The  Raman  cross  sections  for  the  0  and  S  branches  of  H?  and  HD  were  calculated 
from  the  depolarization  ratio  and  the  scattering  strength  dependence  on  the 
rotational  quantum  number.  37  The  ratios  of  the  polarized  S  and  0  branch  dif¬ 
ferential  scattering  cross  sections  to  that  of  the  Q  branch  were  found  to  be 
given  by  the  following  relations: 


(da  /da  )*(J) 
(da  /dr.  )^J) 


2  (J  +  2)  (2J  -  1)  0 

J  (2J  ♦  1) 


(da  /da  )°(J) 
(da  /dQ  )j^jJ 


2  U-  1)  (2J  ±  3)  0(J) 

(J  +  1)  (2J  +1)  1 


(73) 


(74) 


where  p^(J)  is  the  Q  branch  depolarization  ratio  for  linearly  polarized  inci 
dent  light. 


Experimentally  measured  values  of  p  for  H,  give  cross  section  ratios  on  the 

-2  c  (3) 

order  of  10  for  both  S  and  0  branches.  This  implies  that  xv  values  are 

inherently  smaller  for  the  S  and  0  branch  resonances  than  those  of  the  Q  branch 

THG 

resonance  enhancement.  This  property  is  indicated  in  the  results  of  the  x 
calculations  given  in  Table  19. 
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TABLE  19.  TYPE  II  THG  SURVEY  CALCULATION 
(T  =  300  K,  P  =  1  Atm) 


CO  Laser  Free 

luency 

lx<3>l 

(cm® /erg) 

Two-Photon 
Resonance  Molecule 
Branch  (J) 

Transition 

ui  (cm  ^ ) 

P(19)  [12-11] 

1784.33 

2.0  x  10-38 

H2:  0(2) 

P(23)  [  5-4  ] 

1943.99 

1.1  x  10"37 

HD:  S(0) 

P(ll)  [  7-6  ] 

1943.52 

1.2  x  10"37 

HD:  S(0) 

P(10)  [  4-3  ]* 

2025.07 

4.0  x  10"38 

HD:  S(l) 

P(ll)  [12-11] 

1815.81 

1.7  x  10'36 

HD:  Q(0) 

P(  1 8)  [11-10] 

1813.50 

1.3  x  ID*36 

HD:  Q ( 1 ) 

P(19)  [11-10] 

1809,09 

6.6  x  10" 37 

HD:  Q(2) 

P(20)  [11-10] 

1805.29 

2.2  x  10-37 

HD:  Q(3) 

P(26)  [10-9  ] 

1804.76 

1.3  x  10~36 

HD:  Q( 3) 

P(  14)  [12-11] 

1804.28 

2.4  x  10“36 

HD:  Q( 3) 

P(  9)  [  6-5  ]* 

1977.264 

2.4  x  10“39 

HF:  Q(3) 

P(10)  [  6-5  ]* 

1973.285 

2.7  x  10'39 

HF:  Q(4) 

P(ll)  [  6-5  ]* 

1969.274 

8.8  x  10'40 

HF:  Q(5) 

★ 

Line  selected  CO  laser  frequency 


The  only  line  selected  CO  laser  frequency  which  approaches  a  two-photon  reso¬ 
nance  in  HD  is  the  P ( 10 )  [4-3]  transition.  However,  since  the  frequency  detuning 

- 1  THf 

is  approximately  1  cm  from  the  S  branch  resonance  x  is  relatively  small. 

TUP 

The  other  lines  selected  pump  frequencies  yield  smaller  x  values  in  HD  due 

to  off-resonance  conditions.  In  contrast  to  HD  and  H?,  some  of  the  pump  fre- 

^  THf 

quencies  are  close  to  the  Q-branch  resonances  in  HF.  However,  small  x  values 
are  calculated  based  on  the  large  linewidth  data.  Therefore  the  above  results 
suggest  that  Q-branch  resonances  in  HD  should  be  considered  for  optimum  x^ 
resonance  enhancement. 
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3.4.3  Type  III  Scheme 

The  molecular  parameter  references  for  H U  and  DF  are  listed  in  Table  20. 


TABLE  20.  TYPE  III  SCHEME  MOLECULAR  PARAMETER  REFERENCES 


Molecule 

Spectroscopic 

Constants 

Transition 
Dipole  Moment 

Linewidth 

HU 

OF 

[38] 

[39] 

[22] 

[40] 

[23] 

[41] 

The  results  of  the  xTH<*  and  x^  calculations  for  the  line  selected  CO  laser 

frequencies  are  summarized  in  Tables  21  and  22  for  HU  and  DF,  respectively. 

-38 

The  values  listed  in  Tables  21  and  22  are  in  the  range  of  0.8  to  1.5  x  10 
cm^/erg  for  most  of  the  pump  frequency  combinations.  Based  on  the  conversion 
efficiency  estimate  of  Section  3.4.1,  the  type  III  susceptibilities  are  pre¬ 
dicted  to  have  relatively  small  efficiencies. 


TABLE  21.  TYPE  III  x'  '  CALCULATION  RESULTS 


One-Photon  Resonance;  (2)  =  Two-Photon  Resonance;  (3)  =  Three-Photon  Resonance 


TABLE  22.  TYPE  III  x'  CALCULATION  RESULTS 


One-Photon  Resonance;  (2)  =  Two-Photon  Resonance;  (3)  =  Three-Photon  Resonance 


4.0  EFFICIENCY  LIMITATION  AND  PHASE-MATCHING  ANALYSIS 

Efficiency  limitations  in  THG  and  SFG  generally  arise  from  the  breaking  of 
phase-matching  and  the  competing  processes.  Breaking  of  phase-matching  can 
occur  through  various  causes.  These  include  population  transfer  due  to  one- 
and  two-photon  absorption  and  intensity  dependent  refractive  index  variation 
caused  by  quadratic  Kerr  effect.  These  time-dependent  processes  may  ultimately 
limit  pulse  lengths.  Examples  of  competing  processes  are  stimulated  vibrational 
Raman  scattering  of  the  pump  fields.  Another  limitation  is  imposed  by  gas  break¬ 
down  at  very  high  power  and  energy  densities.  This  section  examines  the  impor¬ 
tance  of  the  different  efficiency  limiting  processes.  The  quantitative  analysis 
yields  optimum  operating  parameters  for  efficient  frequency  up-conversion. 

The  breaking  of  phase-matching  limits  conversion  efficiency  because  of  reduced 

effective  interaction  length  for  nonlinear  coherent  processes.  For  a  process 

involving  an  interaction  length  L  and  a  wave  vector  mismatch  Ak,  the  allowed 
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phase  mismatch  is  given  by: 

AkL  £  n  (75) 

The  breaking  of  phase-matching  refers  to  a  condition  that  violates  the  above 
requirement.  Since  a  nonlinear  medium  is  usually  chosen  to  satisfy  the  phase¬ 
matching  condition  initially,  the  breaking  of  phase-matching  may  arise  from 
dynamic  wave  vector  mismatches  caused  by  refractive  index  variations  in  the 
medium  during  the  conversion  process.  In  the  general  case  of  SFG,  the  phase 
variation  is  expressed  by: 

&-[AkL]=  [>s  Sn(uus)  -  u>  5n(wp)  -  6n(wq)  -  iur  6n(u>r)J  |_  (76) 

where  6n(a».)  is  the  refractive  index  variation  at  u  ..  For  the  special  case  of 

J  J 

THG,  the  above  expression  simplifies  to: 

61  AkL]  =  ~  [  6n(3uu  )  -  6n(uu)j  L  (77) 

In  the  following  sections  various  linear  and  nonlinear  processes  which  are 
responsible  for  the  refractive  index  variation  and  the  competing  effects  are 
discussed  in  terms  of  specific  efficiency  limiting  mechanisms. 
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4. 1  Linear  Processes 

The  important  linear  processes  are  single-photon  absorption  at  the  pump  and  the 
generated  frequencies,  population  transfer  effects,  and  thermal  defocusing.  Since 
the  latter  two  effects  are  consequences  of  pump  absorption  in  the  medium,  it  is 
essential  to  choose  molecular  species  which  minimize  the  pump  absorption.  There¬ 
fore,  infrared  absorption  bands  of  the  candidate  molecules  provide  crucial 
parameters  for  the  assessment  of  linear  efficiency  limiting  processes. 

4.1.1  Absorption  Coefficient  and  Refractive  Index  at  Pump  Frequencies 

Single-photon  absorption  of  the  pump  radiation  arises  from  the  fundamental  tran¬ 
sitions  in  the  IR-active  molecules.  These  transitions  also  determine  the  refrac¬ 
tive  indices  at  frequencies  close  to  them.  Both  the  absorption  coefficient  and 
refractive  index  as  a  function  of  pump  frequency  can  be  evaluated  from  the  first- 
order  (linear)  susceptibility  per  molecule  which  is  given  by^ 


The  symbols  have  their  ususal  meanings  as  used  previously  in  tne  tnird-order 

susceptibility  formula.  The  numerical  values  of  were  obtained  from  the 

same  molecular  spectroscopic  constants  and  transition  moment  data  as  were  used 

(3) 

in  the  calculation  of  x  values. 

To  obtain  the  numerical  values  for  the  absorption  coefficient  c*(uu)  and  the 
refractive  index  n(w),  the  following  relations  were  used: 

«(*)  =  Im[Nx(1)(^)]  (79) 

n(u>)  =  [l  +  Re  [Nx(1)  (30) 

In  molecular  gases,  the  real  part  of  Nx^  is  generally  mud-  smaller  than  unity. 
Thus,  a  good  approximation  of  n(w)  is  given  by: 

n(w)  -  1  -  I  Re  [Nx(1)  («.)]  (81) 

A  computer  code  for  x^  has  been  developed,  and  a  program  listing  is  given  in 
Appendix  B. 
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The  results  of  the  <*(ui)  and  n(u>)  -  1  calculations  at  STP  conditions  are  plotted 
in  Fioures  4  through  6  for  the  candidate  molecules  DC£,  DBr,  and  NO,  respectively. 


xlO2 


oa35  n(w)-l 

FIGURE  4(a).  0C£  REFRACTIVE  INDEX  VS  FREQUENCY 


oa35  «(*>) 

FIGURE  4(b).  oa  ABSORPTION  COEFFICIENT  VS  FREQUENCY 
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4  . 


u(cm'l> 
NO  a*u) 


FIGURE  6(b).  NO  ABSORPTION  COEFFICIENT  VS  FREQUENCY 
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When  significant  absorption  is  present  at  m  and  3u>  ,  the  electric  field  equations 
governing  THG  can  be  cast  in  the  following  form: 


dE3 

-  3  E3  -  I 

IT 

1  2 

dEi 

q  r2  1 

~5z 

=  0  E1  -  2 

(82a) 

(82b) 


If  pump  depletion  and  phase  mismatches  are  considered  to  be  negligible,  the 
generated  field  amplitude,  £3(1.) ,  is  given  by 


E3(L) 


S  (Ej(0))3 


(83) 


In  the  case  of  strong  pump  attenuation  (a^L  »  1)  and  <*3  =  0,  the  above  ex¬ 
pression  states  that  the  effective  interaction  length  is  (30^  /2)"1  which  is 
much  less  than  L.  However,  if  attenuation  is  small .  (a^L  «  1) ,  an  approximate 
solution  for  £3(1)  is: 

E3(L)  -  @  (E^O))3  L  (1  -  3ot1  L/4)  (84) 


Thus,  small  attenuation  does  not  significantly  reduce  the  THG  efficiency  if 
dynamic  phase  mismatches  are  negligible. 


4.1.2  Absorption  Coefficient  and  Refractive  Index  at  the  Generated  Frequencies 

Although  the  generated  frequencies  are  far  from  the  fundamental  IR  transition 
bands,  they  nearly  coincide  with  the  overtone  transitions  in  some  of  the  candi¬ 
date  molecules.  The  calculation  of  a(uo)  and  n(uj)  -  1  is  straightforward  using 
the  x(1)  computer  code.  Figures  7  and  8  show  the  results  of  the  calculations 
for  the  first  overtone  bands  of  DF  and  HC£ ,  respectively. 


The  second  overtone  bands  of  the  type  I  scheme  molecules  are  also  close  to  the 

(3) 

generated  frequencies  and  are  important  in  the  resonance  enhancement  of  x  » 

however,  single-photon  absorption  at  the  generated  frequencies  is  exceedingly 

-5-1 

small  (typically  a  <10  cm  ).  Therefore,  a  model  discussed  in  the  previous 
section,  in  which  <*3  =  0  was  assumed,  is  a  realistic  description  of  the  type  I 
scheme. 
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5K0 


— ,  ■  i 

3cC0 


FIGURE  8(b). 
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HU  ABSORPTION  COEFFICIENT  VS  FREQUENCY 
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In  contrast  to  the  type  I  scheme,  absorption  in  the  type  III  scheme  is  signifi¬ 
cant  only  at  the  generated  frequencies.  By  setting  =  o  in  Equation  (83),  the 
generated  field  amplitude  becomes 

E3(L>  =  If  (l  -  e"*3172)  (85) 

In  the  strong  attenuation  limit  (a,L  »  1),  the  effective  interaction  length 
-1  ^ 

becomes  {a^/2)  .  On  the  other  hand,  if  a^L  «  1,  ( L)  is  proportional  to 

(E^O))3  and  L.  From  these  properties,  it  is  evident  that  a  small  attenuation 
of  the  generated  field  is  desirable  and  that  the  attenuation  only  slightly  re¬ 
duces  the  conversion  efficiency.  These  conclusions  together  with  those  of  the 
previous  section  are  valid  only  to  the  extent  in  which  dynamic  phase  mismatches 
can  be  neglected.  It  will  be  shown  in  the  next  section,  however,  that  the 
breaking  of  the  phase-matching  condition  due  to  absorption  is  much  more  sensi¬ 
tive  and  critical  than  pump  attenuation  in  the  IR-active  candidate  molecules. 

4.1.3  Population  Transfer  and  Dynamic  Phase  Mismatch 

Whenever  ground  state  molecules  absorb  pump  radiation,  a  certain  fraction  of 
the  population  is  transferred  to  the  first  excited  vibrational  level.  Since 
the  fundamental  IR  transitions  determine  the  refractive  indices  at  the  pump 
frequencies,  population  transfer  perturbs  the  refractive  indices  at  the  pump 
frequencies.  This  immediately  leads  to  phase  variation  through  expressions 
(76)  and  (77). 

Because  population  transfer  is  time-dependent,  the  resulting  phase  variation  is 
a  dynamic  process.  Although  this  dynamic  phase  mismatch  can  be  calculated  in 
principle,  such  a  calculation  would  require  an  estimate  of  population  transfer 
based  on  the  radiative  and  kinetic  properties  of  the  gas,  which  are  beyond  the 
scope  of  this  work.  Therefore,  a  quasi -steady  state  model  was  considered  for 
the  purpose  of  illustrating  the  effects  of  the  dynamic  phase  mismatch.  The 
basic  assumption  of  this  model  is  that  the  rotational  states  of  the  first  ex¬ 
cited  vibrational  state  are  in  thermal  equilibrium  at  the  gas  temperature.  The 
excited  vibrational  state  population  as  a  whole  is  considered  to  be  time-depen¬ 
dent  in  order  to  allow  for  pulse  length  effects  and  pump  intensity  fluctuations. 


V 
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The  refractive  Index  variation  calculations  using  the  code  were  accomplished 
by  specifying  a  certain  relative  population  transfer  value  (AN/N).  For  example, 
in  the  case  of  DG£  with  (aN/N)  =  10"  ^  at  STP,  different  pump  frequencies  yielded 

-4  -7  1 

wave  vector  mismatches  from  3.2  x  10  to  6.5  x  10  cm"  .  These  values  are  typi 

cal  of  the  type  I  scheme  molecules.  For  a  practical  device  having  an  interaction 

_2 

length  of  300  cm,  the  corresponding  phase  variations  are  9.6  x  10  to  19.5 
radians,  thereby  indicating  that  breaking  of  phase-matching  is  a  serious  effect 
for  some  pump  frequencies  and  significant  for  many  others.  If  (AN/N)  *  10  is 
assumed,  most  of  the  pump  frequencies  would  experience  severe  phase  mismatches. 
This  means  that  population  transfer  values  of  as  little  as  0.1%  can  lead  to  sig¬ 
nificant  efficiency  reduction  and  that  ll  transfer  levels  are  beyond  the  toler¬ 
ance  range.  Conversely,  efficient  frequency  conversion  can  be  maintained  only 
over  a  time  interval  during  which  population  transfer  can  be  held  to  less  than 
0.1%. 


An  estimate  for  the  allowable  pulse  duration  AT  can  be  obtained  from  the  follow¬ 
ing  relation: 


AN  _  al  AT 

TT~  "  UE¥ 


(86) 


where  I  is  the  pump  intensity  and  a  is  the  absorption  coefficient.  For  most  of 

-23  2 

the  type  I  scheme  molecules  a/N  is  typically  in  the  range  of  10  cm  (off- 
-19  2 

resonance)  to  10  cm  at  the  peaks  of  absorption  lines.  Using  the  lower  limit, 

2  -9 

the  maximum  value  of  AT  for  I  =  0.5  GW/cm  is  found  to  be  8  x  10  second  if 

-3 

(aN/N)  is  restricted  to  10  .  Therefore,  at  the  pump  intensities  of  interest 

for  high  conversion  efficiencies,  the  pulse  lengths  are  limited  to  the  nano¬ 
second  time  scale. 


4.1.4  Thermal  Defocus inq 

The  population  transfer  may  cause  thermal  defocusing  which  is  a  manifestation 
of  refractive  index  variation  across  the  beam  for  a  nonuniform  pump  beam  cross 
section.  In  the  case  of  a  gauss i an  beam  profile  with  an  initial  beam  radius, 
Wq,  the  beam  radius  changes  by  a  factor  of  two  when  the  refractive  index  on  the 
beam  axis  varies  by  an  amount  given  by: ^ 

5n  =  w2/2L2  (87) 
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where  L  is  the  Interaction  length.  If  wQ  is  taken  to  be  the  waist  such  that 
Wq  a  \b/2n  (b  *  confocal  parameter),  then  the  condition  (87)  becomes: 

6n  =  A.b/4n  L2  (88) 

For  the  case  of  confocal  focusing  (b  s  |_),  the  numerical  value  of  6n  is  1.3x  10”7 

for  L  =  300  era.  From  the  results  of  the  previous  section,  (aN/N)  »  10  in  DC* 

-8  -6 

yields  6n  values  in  the  range  of  2.5  x  10  to  5.2  x  10  .  This  clearly  shows 
that  thermal  defocusing  has  a  serious  effect  on  the  efficiency  due  to  reduced 
pump  intensities  in  the  nonlinear  medium.  However,  it  is  also  important  to  note 
that  for  large  beam  cross  sections  and/or  uniform  intensity  profiles,  thermal 
defocusing  effects  can  be  minimized. 


4.2  Nonlinear  Processes 

A  number  of  third-order  processes  besides  THG  and  SFG  may  be  present  at  high 
pump  intensities.  These  include  quadratic  Kerr  effect,  two-photon  absorption, 
and  stimulated  Raman  scattering.  Other  nonlinear  processes  such  as  three-photon 
absorption,  gas  breakdown,  and  pump  depletion  may  also  affect  conversion  effi¬ 
ciencies.  Quantitative  assessment  of  each  of  these  intensity-dependent  processes 
is  discussed  in  the  following  sections. 


4.2.1  Two-Photon  Absorption 

In  the  type  I  and  II  resonance  enhancement  schemes,  small  two-photon  detunings 

are  used  to  obtain  large  x^  values.  Because  of  such  near  resonances,  the 

probability  of  two-photon  absorption  (TPA)  can  be  very  significant.  As  in  the 

single-photon  absorption  case,  two-photon  absorption  can  reduce  conversion  effi 

ciency  due  to  pump  attenuation,  dynamic  phase  mismatch,  and  thermal  defocusing. 

(3) 

Thus,  it  is  clear  that  a  compromise  between  large  x  enhancement  and  two- 
photon  absorption  should  be  considered. 


(2) 

For  relatively  small  TPA,  the  attenuation  parameter,  a'  'L,  is  defined  by 


.<2>l 


1  -  Ip  U)/IP(0) 


(89) 


In  the  case  of  self-absorption,  the  attenuation  parameter  at  the  pump  frequency, 

a(2)l-  =  }»p  T,.2  [NL  Ip(0)][xJPA  («>p)]  (90a) 


(Up,  is  given  by: 
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where  xJPA(*>p)  =  Im  Xp3^  (®p.  -uy  *p)  (90b) 

If  photons  of  unequal  frequencies  cause  TPA,  the  induced  absorption  at  w  is 
characterized  by: 

a{2)L  =  3uupTl22^CNL  Ij(0)3[xJPA  (^j)3  Ola) 

J 

where  xpPA(u»j)  =  Imxp3^  (wp, -U)j »  “j)  for  j  ^  P  (91b) 

(2) 

A  numerical  estimate  for  a'  'L  at  typical  operating  conditions  can  be  obtained 

by  assuming  w  =  3.772  x  10*A  rad/sec  (i.e.,  2000  cm‘^),  N  =  102^  m-3,  L  =  1  m, 

P  ip  2 

and  Ip(0)  =  5  x  10  M/m  .  A  representative  value  for  (two-photon  resonant) 

xTPA  is  10" 34  cm®/erg  which  is  equivalent  to  1.235  x  10’^  MKS.  With  these 
(2) 

parameters,  av  'L  is  equal  to  0.498,  or  nearly  half  of  the  pump  radiation  is 

absorbed.  Therefore,  two-photon  absorption  is  a  significant  efficiency  limiting 

process  if  x™  values  are  on  the  order  of  10”3A  cm6/erg.  For  the  CO  laser 

TPA 

frequencies,  calculations  show  that  x  values  approach  this  magnitude  only  in 
the  type  II  scheme. 


4.2.2  Quadratic  Kerr  Effect 

Another  process  intimately  related  to  TPA  is  quadratic  Kerr  effect  which  causes 

intensity  dependent  refractive  index  variations  through  the  real  part  of 
(31 

x p  (o)p,  -iDp.ujp).  The  principal  effect  of  this  process  involves  phase  shifts 
along  the  direction  of  propagation,  which  can  cause  breaking  of  phase-matching. 
For  a  nonuniform  beam  intensity,  the  refractive  index  variations  also  lead  to 
beam  divergence  (Kerr  defocusing).  However,  since  beam  uniformity  can  be  con¬ 
trolled  to  a  certain  extent,  the  Kerr  phase  shifts  are  inherently  more  important. 

The  definition  of  the  Kerr  phase  shift  at  the  pump  frequency  is  given  by 

V^err(»pn  ,92.) 


where 


x 


Kerr 

P 


(uip)  *  Re 


'V  V 


(92b) 
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In  the  presence  of  other  pump  frequencies,  the  Kerr  phase  shift  at  is  de¬ 
fined  by 

(AkL)Kerr  =  7“i  £NL  ^  Cx*err  («*j )3  (93a) 

j 

where  x^6™*^)  =  Rex^(u>j,  -“>j,  u^.)  i  f  j  (93b) 

The  above  expression  also  describes  the  Kerr  phase  shift  at  the  generated  fre¬ 
quencies,  j>s,  when  uus  is  substituted  form...  In  the  case  of  SFG,  all  phase 
shifts  due  to  various  frequency  combinations  must  be  summed  using  (93a,  b). 

Since  the  quadratic  Kerr  effect  vanishes  (aside  from  small  nonresonant  contri¬ 
butions)  at  the  exact  two-photon  resonance,  the  Kerr  phase  shifts  can  be  made 
negligible,  in  principle,  by  an  appropriate  choice  of  pump  frequencies.  How¬ 
ever,  for  CO  frequency  tripling,  the  two-photon  detunings  are  generally  non¬ 
zero.  Thus,  considerable  phase  shifts  may  be  encountered.  For  example,  if  the 
operating  conditions  for  the  TPA  calculation  in  the  previous  section  are  assumed 
and  if  the  value  of  the  Kerr  nonlinearity  at  the  pump  frequency  is  assumed  to 
be  5  x  IQ'35  cm6/erg  (i.e.,  peak  of  xKerr  when  xTPA  =  10"34  cm6/erg  on  resonance), 
the  resulting  phase  shift  is  0.125  radian.  While  this  shift  is  smaller  than  tt  , 
the  predicted  phase-matched  efficiency  is  only  4.7%  (neglecting  pump  attenuation) 

TUP  [/-wm 

if  the  magnitude  of  x  were  equal  to  x  .  This  suggests  that  a  large  ratio 
of  l x  l /l x  e  1  is  desirable  in  order  to  maximize  the  THG  efficiency  without 
the  Kerr  induced  phase  mismatches. 

In  the  type  I  resonance  enhancement  scheme,  x  can  be  calculated  with  the 

x{3)  code.  For  the  line  selected  CO  laser  frequencies  |xn6rr  (u>  )|  values  are 

THGi  P  P  Kerr  i 
typically  two  orders  of  magnitude  larger  than  |x  I  .  In  contrast  |x.  («>  )| 

TUP  5  r 

values  are  an  order  of  magnitude  smaller  than  Jx  !  •  These  observations  indi¬ 
cate  that  the  Kerr  induced  phase  mismatches  at  the  pump  frequencies  are  likely 
to  limit  the  THG  efficiencies  in  the  type  I  scheme. 

The  role  of  the  Kerr  induced  phase  mismatches  in  the  type  II  scheme  is  important 
due  to  the  fact  that  all  the  third  order  processes  are  related  through  a  common 
two-photon  resonance  enhancement.  For  example,  if  the  detuning  is  large  with 
respect  to  the  linewidth,  |xKerr|  is  equal  to  |xTH<3[/3.  On  the  other  hand,  at 
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exact  two-photon  resonance  (Aw  =  0),  x  vanishes  and  lx  I  is  equal  to 
3  x  lxTPA| .  At  intermediate  detunings,  both  the  Kerr  phase  shift  and  TPA  pump 
attenuation  must  be  considered  together.  Such  a  treatment*"*  has  shown  that 

high  efficiencies  are  possible  with  phase  shifts  close  to  n/2  while  holding 

f  31 

TPA  losses  to  about  20%.  Although  the  x'  '  calculations  for  the  type  II  scheme 
indicate  phase  shifts  much  smaller  than  n/2  at  the  CO  laser  pump  frequencies 
due  to  relatively  large  detunings,  analysis  shows  that  if  a  near  two-photon 
resonance  can  be  achieved  on  a  Raman-active  Q-branch  transition,  the  desired 
phase  shifts  are  feasible.  This  point  will  be  discussed  further  in  connection 
with  multiline  SFG  in  HD. 

In  the  type  III  resonance  enhancement  scheme,  quadratic  Kerr  effect  at  the  pump 
and  the  generated  frequencies  can  be  expected  to  be  small.  This  is  one  of  the 
advantageous  features  of  the  scheme  when  it  was  considered  initially.  A  further 
analysis  has  shown  that  the  Kerr  effect  introduces  a  level  shift**  of  the  three- 
photon  resonance  according  to 

Ul2  H  Ip 

6uj  ~  -  ~  (94) 

hz  a 

where  n  is  the  vibronic  transition  frequency  and  *  is  the  associated  matrix 

between  the  vibronic  state  and  the  v  =  2  level  of  the  ground  electronic  state. 

In  principle,  all  such  transitions  must  be  included  in  calculating  the  level 

shift,  but  approximation  based  on  a  single  strong  transition  yields  an  order  of 
.  -2  -1 

magnitude  estimate  of  10  cm  at  the  typical  operating  conditions.  Since  this 
is  a  small  shift  relative  to  the  transition  linewidth,  resonance  enhancement  is 
not  affected  appreciably.  Also,  the  Kerr  induced  phase  mismatch  can  be  estimated 
using  the  following  relations: 


3  w  AnL 


(AkL)Kerr  = 

_ E _ 

n 

(95a) 

with 

An  ~ 

NIh^qI  ^  6uj 
.0 

(95b) 

h{  A  w)> 

where  represents  the  first  overtone  transition  matrix  element  and  Aw  is  the 

detuning  which  is  assumed  to  be  large  compared  to  the  linewidth.  Substituting 

representative  parameters  of  the  type  III  molecules  into  the  above  expressions 

-2  -1  -1 

yields  a  phase  shift  of  0.03  radian  if  *  10  cm  and  Aw  =  1  cm  (1.886  x 
1011  rad/sec)  are  used.  Thus,  breaking  of  the  phase-matching  condition  due  to  quad' 
ratic  Kerr  effect  is  not  believed  to  be  a  serious  limitation  in  the  type  III  scheme 
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4.2.3  Stimulated  Raman  Scattering 

Stimulated  Raman  scattering  (SRS)  in  molecular  gases  may  involve  vibrational 
and/or  rotational  Raman  transitions.  In  the  vibrational  SRS  process  IR  fre¬ 
quencies  up-converted  via  THG  or  SFG  become  the  "pump"  radiation.  In  the  rota 
tional  SRS  process,  both  the  pump  and  the  up-converted  radiation  can  act  as 
source  terms  for  amplifying  the  respective  Stokes  radiation.  In  either  SRS 
processes,  the  THG  and  SFG  efficiencies  will  be  limited  through  a  competition 
of  frequency  conversion  processes. 


A  key  parameter,  which  indicates  when  the  SRS  processes  become  important,  is 
the  stimulated  Raman  gain  coefficient  g$.  This  parameter  is  defined  by:45 


9S 


4  AN 
^R 


(96) 


where  AN 


population  density  difference 


Stokes  radiation  wavelength 

Raman  transition  linewidth  (FWHM)  in  rad/ sec 

pump  intensity 

pump  frequency  in  rad/sec 

differential  Raman  scattering  cross  section 
for  the  pump  radiation 


The  empirical  values  of  the  various  parameters  are  generally  available,  but 
the  cross  section  data  are  usually  measured  for  the  Stokes  raidation  in  the  vis 
ible  regi on ,  These  data  can  be  scaled  to  the  IR  frequencies  by  making  use 
of  the  w4  dependence  of  (do/d*,  )s  and  the  relationship 


(97) 


Substituting  the  expression  (97)  into  (96),  one  obtains 


9s 


(98) 
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where  the  quantity  inside  the  brackets  is  nearly  invariant  at  Stokes  frequencies 

in  the  visible  and  IR  regions.  Thus,  by  evaluating  the  invariant  quantity  from 

visible  Raman  scattering  data,  Raman  gain  coefficients  can  be  calculated  for  the 

IR  region.  Table  23  lists  the  results  of  the  calculations  based  on  a  AN  value 

of  1  amagat  (2.69  x  10  m  )  for  the  vibrational  SRS  gain  at  an  intensity  of 
2  13  2 

1  GW/cnr  (10  W/m4).  In  these  calculations  the  pump  frequency  is  assumed 

to  be  6000  cm"1.  Since  not  all  cross  section  data  were  available,  gs  estimates 
for  the  deuterated  molecules  (i.e.,  HD,  DCA,  DBr,  and  OF)  were  scaled  from  those 
of  the  normal  species  by  the  ratio  of  the  Stokes  frequencies. 

The  threshold  gain  for  Raman  oscillation  is  usually  considered  to  be  gsL  =  30. 

For  practical  device  lengths,  the  gs  values  in  Table  23  indicate  that  vibrational 
SRS  oscillation  can  be  regarded  as  a  small  effect  in  all  the  molecules  except  for 
H£.  The  hydrogen  gas  is  a  special  case  in  which  gs  increases  with  higher  pres¬ 
sures  (up  to  20  amagats).  For  example,  at  8  amagats,  the  linewidth  is  the  same 
31 

as  that  at  1  amagat,  thereby  increasing  the  gs  value  by  a  factor  of  8. 

Another  effect  of  vibrational  SRS  is  Raman  amplification  of  the  pump  radiation 
by  the  generated  radiation.  By  amplifying  the  pump  ''adiation  the  generated 
radiation  transfers  the  up-converted  power  back  into  the  pump  radiation.  This 
process  acts  as  a  stabilizing  mechanism  to  establish  a  steady-state  efficiency 
condition.  An  analysis  of  this  THG/SRS  interaction  shows  that  a  characteristic 
parameter  is  given  by:4"1 


When  the  value  of  *  is  close  to  or  larger  than  unity,  the  SRS  process  is  not 
important;  efficiency  will  be  limited  by  other  processes.  However,  if  k  1, 
then  the  efficiency  limit  is  equal  to  the  value  of  *.  For  the  type  I  scheme 
molecules,  '  code  calculations  yield  |x  |  values  which  are  an  order  of 
magnitude  smaller  than  |x  |.  This  indicates  that  vibrational  Raman  amplifi¬ 
cation  in  the  type  I  scheme  is  not  a  significant  limiting  process.  In  contrast, 
x  *  0.25  on  resonance  (Auu  =  0)  in  the  type  II  scheme,  which  suggests  that 
vibrational  Raman  amplification  is  an  important  factor  governing  the  maximum 
THG  efficiency. 
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Rotational  Raman  scattering  is  a  third-order  process  that  does  not  share  any  of 
the  resonance  enhancement  paths  used  in  the  THG  and  SFG  processes.  As  such, 
the  rotational  SRS  process  represents  a  potential  competing  loss  mechanism  for 
the  pump  and  the  generated  radiation.  As  in  the  vibrational  SRS  case,  a  key 
parameter  in  the  rotational  SRS  is  the  Raman  gain  coefficient.  Expression  (98) 
is  still  a  valid  formula  for  calculating  g$,  but  AN  must  be  defined  in  terms  of 
an  actual  population  difference  because  of  much  smaller  Raman  shifts  in  the 
rotational  case.  A  formula  for  AN  as  a  function  of  the  rotational  quantum 
number  of  the  initial  state  is  given  by 


where 


ANj  =  N(fj  -  fJ+2) 


(100) 


total  number  density 

fraction  of  molecules  in  the  initial  state 
fraction  of  molecules  in  the  final  state 


Under  equilibrium  conditions  at  a  temperature  T,  the  population  distribution 
can  be  expressed  by: 

fi'  -  2  +  1  -E(J')/kT  (101) 

U  -  Q  e 

where  Q  =■  kT/hc  B0  (rotational  partition  function) 

and  E(j')  -  he  BQ  J'  (J'+  1) 

It  is  evident  from  the  above  formulas  that  population  differences  become  smaller 

with  smaller  values  of  B  .  For  example,  J  =  6  -*  J  =  8  S(6)  transition  in  CO  at 

0  -3 

room  temperature  has  a  calculated  value  of  ANj_g/N  =  1.41  x  10  .  Since  the 

rotational  Raman  cross  sections  are  roughly  comparable  to  the  vibrational  Raman 

cross  sections,^  the  g$  values  for  the  rotational  SRS  are  smaller  than  the 

vibrational  Raman  gains  by  the  AN/N  ratio.  A  survey  of  the  candidate  molecules 

shows  that  only  H2  has  sufficiently  large  rotational  Raman  shifts  to  make  AN/N 

ratios  close  to  those  of  the  vibrational  Raman  shifts. 

46  47 

Recent  experimental  demonstrations  ’  of  rotational  SRS  in  H2  (SQ(0))  using 
CC>2  laser  radiation  indicate  that  rotational  Raman  shifting  can  be  a  very 
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efficient  IR  frequency  conversion  process.  Based  on  the  information  given  in 
these  experiments,  the  rotational  Raman  gain  coefficient  is  predicted  to  be  in 
the  range  of  11  to  20  m"*  at  a  CO  laser  pump  (single-line)  intensity  of  1  GW/cm^. 
The  predicted  value  of  g$  implies  efficient  rotational  Raman  shifting  for  inter¬ 
action  lengths  in  excess  of  3  m.  Furthermore,  the  Hj  ( Sq( 1 ) )  transition  may  also 
be  important  at  room  temperature  with  a  gs  value  comparable  to  that  of  the  Sq(0) 
transition. 

Although  H£  presents  a  potentially  serious  efficiency  limitation  due  to  rotational 
SRS,  other  type  II  molecules  (HD  and  HF)  appear  to  have  a  lesser  problem.  This 
is  because  of  smaller  AN  values  and  larger  Raman  linewidths  in  HD  and  HF  com¬ 
pared  to  that  of  H2,  both  of  which  reduce  the  gs  values  sufficiently  to  prevent 
Raman  oscillation. 


4.2.4  Three-Photon  Absorption 

Three-photon  absorption  (ThPA)  is  a  nonlinear  process  governed  by  a  fifth-order 
susceptibility.  Because  ThPA  is  a  higher-order  effect  its  presence  is  negli¬ 
gible  unless  exceptionally  small  frequency  detunings  or  high  pump  intensities 
are  employed.  Thus,  the  only  resonance  enhancement  scheme  amenable  to  ThPA  is 
the  type  I  scheme  with  near-resonant  detunings  at  one-,  two-,  and  three-photon 
resonances. 


In  the  near-resonant  case,  three-photon  transition  probability  per  unit  time 
can  be  written  as4® 


W(ThPA) 


*3  3 

tt  9j_  (3  (d)  Tj  I  /  ^  32  ^21  **  10  \2 
256  h6  n3  \  Au^Au^  / 


(102) 


where  gL(3w)  is  the  ThPA  line  shape  function  and  all  other  parameters  have 

their  usual  meanings.  This  expression  can  be  compared  with  two-photon  transi 
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tion  probability  per  unit  time  as  given  by: 


W(TPA) 


ngL  (2ui)  T|2  I2  /^21  m>10\  2 
32  ii4  n2  \  A®  j  / 


(103) 


If  the  two-photon  detuning  is  substantially  larger  than  the  linewidth,  g,  (2u>) 

2  L 
can  be  approximated  by  r^/hAiu^.  Similarly,  g^  (3u>)  is  approximately  equal 
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to  rynAoij  when  Aui^  »  With  these  approximations,  the  ratio  of  the  transi¬ 
tion  probabilities  may  be  expressed  as: 


(104) 


For  most  type  I  molecules,  the  r^/I^  ratio  is  close  to  unity  and  |p.,J  1S  on  the 
order  of  0.1  debye.  Using  I  =  1  GW/cm^  (10^  W/m^)  and  Aw,  =  l  cm“*  (1.886  x 
10  rad/sec),  expression  (104)  yields  a  ratio  of  0.13.  Since  three-photon 
detunings  in  the  type  I  molecules  are  generally  larger  than  1  cm-1  for  various 
combinations  of  CO  laser  frequencies,  the  above  ratio  indicates  that  ThPA  effects 
are  an  order  of  magnitude  smaller  than  those  of  TPA. 


4.2.5  Gas  Breakdown 

Laser  induced  gas  breakdown  is  a  well  known  phenomenon  with  high  peak  power 
laser  beams  and  is  a  major  consideration  in  the  analysis  of  efficiency  limita¬ 
tions.  Since  beam  propagation  is  severely  affected  when  gas  breakdown  occurs, 
the  critical  intensity  at  which  the  breakdown  is  initiated  sets  the  upper  limit 
on  the  operating  intensities  for  a  frequency  converter. 

A  literature  survey  showed  that  laser  induced  gas  breakdown  data  are  not  avail¬ 
able  for  the  case  of  CO  laser  radiation  interacting  with  the  particular  THG 
candidate  molecular  gases.  Most  of  the  experimentally  observed  breakdown 
intensities  in  the  IR  are  reported  for  the  10  p.m  radiation.  For  example,  at 

3  atmospheres  of  H,  the  reported  experimental  conditions4®  imply  a  breakdown 
c  1 

power  density  of  2.8  GW/ cm  at  10  The  corresponding  breakdown  energy  den- 

2 

sity  would  be  196  J/cm  .  If  the  quadratic  frequency  dependence  of  breakdown 

4q  2 

intensity  3  is  assumed,  the  10  p-m  radiation  data  suggest  a  value  of  11  GW/cm 

for  the  5  p.m  radiation  provided  all  other  parameters  are  unchanged.  However, 

since  semi -empirical  models4^’®0  of  gas  breakdown  predict  different  behavior  for 

short  and  long  laser  pulse  lengths,  estimates  using  scaling  arguments  based  on 

10  p.m  data  should  be  regarded  as  a  rough  indication.  Therefore,  breakdown 

Intensities  based  on  the  presently  available  information  may  be  reasonably 

2 

expected  to  be  in  the  range  of  1  to  10  GW/cm  .  These  intensity  limits  in  turn 

TUT 

imply  minimum  values  of  x  for  efficient  conversion.  At  the  lower  limit, 

TUP  fi  1 

lx  lvalues  in  excess  of  10  cm  /erg  are  necessary,  which  establishes  a 

criterion  for  selecting  particular  candidate  molecules. 
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4.2.6  Pump  Depletion 


Efficiency  limitation  analysis  so  far  has  been  presented  with  the  assumption 
that  pump  radiation  does  not  change  appreciably  in  order  to  illustrate  various 
limitations.  In  reality,  frequency  conversion  processes  do  involve  pump  de¬ 
pletion  by  definition,  but  an  analytical  treatment  becomes  complex  because  of 
nonlinear  coupled  differential  equations  governing  the  generation  and  depletion 
processes.  In  general,  these  equations  do  not  have  closed  form  solutions  and 
therefore  require  numerical  methods  for  analyzing  the  effects  of  pump  depletion. 
However,  rather  than  computing  various  numerical  solutions  encompassing  all 
aspects  of  interacting  processes,  a  simple  model  including  only  the  THG  and 
pump  depletion  processes  has  been  analyzed.  This  approach  has  an  advantage  of 
yielding  an  analytical  solution  to  aid  in  understanding  the  efficiency  when 
there  is  substantial  pump  depletion.  The  equations  of  this  model,  neglecting 
phase  mismatches,  can  be  cast  in  the  following  form: 


where  3 


i  u>,r,N  xTHG 

TT2 


-  3 


bl  j 


(105) 

(106) 


Combining  the  Equations  (105)  and  (106)  yields  the  statement  of  energy  conservation 


Ef(0)  -  Ef (Z)  =  E|(Z) 


(107) 


Substitution  of  this  relation  into  (105)  and  a  direct  integration  gives  the  con¬ 
version  efficiency  e: 


where 


.  _  E3<L>  2 

"  ^loT 

■  \ 

TTHq 

\  -  (p  e|(0)  L)2  -  ixTHG| 2  (NIL)2 


(108) 


For  small  pump  intensities,  m  <k  l,  and  the  efficiency  is  approximately  given 
by  the  value  of  Tg,  which  is  exactly  the  efficiency  one  would  obtain  by  inte¬ 
grating  (105)  assuming  no  pump  depletion.  Hence,  the  efficiency  figures 
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calculated  without  pump  depletion  are  still  useful  in  predicting  the  actual 
efficiencies  provided  that  expression  (108)  is  used.  The  relation  (108)  also 
shows  that  efficiencies  below  10%  can  be  calculated  quite  accurately  without 
including  pump  depletion  effects. 

At  very  high  pump  intensities  for  which  Tig  >  1,  THG  conversion  efficiency  is  no 
longer  increasing  quadratically  with  pump  intensity.  Thus,  pump  depletion  leads 
to  an  efficiency  saturation  behavior  as  described  by  (108).  Because  SFG  and  THG 
are  basically  related  processes,  similar  saturation  effects  can  be  expected  for 
SFG  efficiencies. 

4.3  Optimum  Scheme 

The  discussions  of  various  efficiency  limiting  processes  and  the  calculations 
of  key  parameters  in  the  previous  sections  can  be  summarized  qualitatively  as 
shown  in  Table  24.  When  a  particular  process  causes  a  substantial  decrease 
(e.g.  more  than  10%)  in  the  actual  efficiency  relative  to  the  predicted  effi¬ 
ciency  in  the  absence  of  the  process,  that  process  is  considered  to  be  signi- 
cant.  If  the  effect  leads  to  efficiency  losses  on  the  order  of  a  percent  for 
some  candidate  molecules  of  a  given  enhancement  scheme,  then  the  limiting 
process  is  considered  moderate.  Processes  which  are  predicted  to  produce  less 
than  one  percent  changes  in  the  conversion  efficiency  are  regarded  as  small, 
and  those  with  no  measurable  consequences  are  considered  negligible.  The  table 
aids  in  illustrating  the  essential  efficiency  limiting  characteristics  for  each 
of  the  three  resonance  enhancement  schemes.  Based  on  these  characteristics, 
the  following  observations  can  be  made: 

In  the  type  I  resonance  enhancement  scheme,  pump  absorption  leading  to  dynamic 
phase  mismatches  is  found  to  be  the  primary  efficiency  limiting  process.  Since 
this  scheme  requires  near-resonant  conditions  at  one-,  two-,  and  three-photon 
detunings,  the  pump  absorption  is  always  present  and  inherently  restricts  the 
application  of  the  type  I  scheme  to  short  pump  pulses. 

The  type  II  scheme  avoids  dynamic  phase  mismatch  caused  by  pump  absorption. 

Other  competing  nonlinear  processes  become  significant  at  the  high  pump  inten¬ 
sities  required  for  efficient  conversion,  however.  Although  two-photon  ab¬ 
sorption,  quadratic  Kerr  effect,  and  stimulated  Raman  scattering  ultimately 
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TABLE  24.  EFFICIENCY  LIMITATION  ANALYSIS 
SUMMARY 


NONLINEAR  EFFECTS: 

•  Two-Photon  Absorption 

•  Quadratic  Kerr  Effect 

•  Stimulated  Raman  Scattering 

(Vibrational ) 

•  Stimulated  Raman  Scattering 

(Rotational) 

•  Three-Photon  Absorption 

•  Gas  Breakdown 

•  Saturation 

LINEAR  EFFECTS: 

•  Absorption  Loss  at  Pump 

Frequency 

•  Absorption  Loss  at  Sum 

Frequency 

•  Population  Transfer- Induced 

Phase  Mismatch 

•  Thermal  Defocusing 


|  |  Negligible  / 


Small 


RESONANCE  ENHANCEMENT  SCHEME 

TYPE  I 

TYPE  II 

TYPE  III 

□ 


□  □ 


□ 


Moderate 


Significant 


v 
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limit  the  efficiency,  near  unity  THG  and  SFG  conversion  efficiencies  are  pre¬ 
dicted  under  optimum  conditions.  The  parameters  that  govern  these  conditions 
are  the  two-photon  resonance  frequency  and  the  values  of  Kerr-induced  phase 
shifts. 

(31 

The  type  III  scheme  analysis  and  the  x  calculations  show  that  relatively 
small  conversion  efficiencies  can  be  expected.  The  principal  efficiency  limita¬ 
tion  in  this  scheme  is  optical  gas  breakdown  at  high  pump  intensities.  The  high 
pump  intensities  are  necessary  to  compensate  for  the  relatively  small  nonlinear 
susceptibilities  calculated  for  candidate  molecules  of  this  scheme. 

In  view  of  the  efficiency  limiting  properties  of  the  three  enhancement  schemes, 

the  type  II  scheme  provides  a  most  promising  technique  for  up-converting  CO 

laser  frequencies.  In  particular,  H2  and  HD  have  an  added  advantage  of  extremely 

low  dispersion  in  the  infrared  region.  The  refractive  indices  for  H~  can  be 

51  c 

calculated  from  a  Sellmeir  formula,  and  the  calculated  THG  and  SFG  coherence 
length  (distance  at  which  AkL  -tt)  at  1  amagat  is  approximately  300  cm.  This 
value  is  not  very  sensitive  to  the  range  of  the  line-selected  CO  laser  frequen¬ 
cies,  which  permits  simultaneous  phase-matching  for  various  combinations  ©*  pump 
frequencies.  Since  H2  and  HD  are  essentially  identical  in  their  electronic 
properties,  dispersion  properties  can  be  expected  to  be  similar.  Thus,  both 
H?  and  HD  appear  to  be  most  suited  for  multiline  THG  or  SFG. 


5.0 


MULTILINE  THG  AND  SFG  ANALYSIS 


The  CO  laser  output  consists  of  several  frequencies,  as  elucidated  In  Section 
3.4.  Thus,  efficient  THG  and  SFG  processes  for  a  number  of  lines  are  required. 
The  phase-matching  requirement  for  various  frequency  combinations  is  an  im¬ 
portant  consideration  in  multi  frequency  THG  and  SFG.  For  example,  it  is 
desirable  to  phase-match  all  of  the  CO  laser  lines  referred  to  in  Section  3.4 
for  THG  in  a  given  gas  mixture.  The  THG  and  SFG  conversion  efficiencies  of  the 
various  frequency  combinations  may  differ  from  one  particular  set  to  another 
because  of  the  frequency  dependence  of  the  nonlinear  susceptibility.  There¬ 
fore,  a  scheme  in  which  is  near  optimum  and  relatively  insensitive  to 
changes  in  the  pump  frequencies  is  needed.  In  the  following  section,  calcula¬ 
tions  on  candidate  molecules  for  each  of  the  three  resonance  enhancement  schemes 
are  discussed  in  relation  to  the  criteria  outlined  above.  Accordingly,  an 
appropriate  scheme  for  SFG  is  suggested. 

5.1  Efficient  THG  and  SFG  Optimization 

Efficient  THG  and  SFG  of  the  line  selected  CO  laser  frequencies  requires  opti¬ 
mization  of  the  nonlinear  susceptibility  and  phase-matching  properties  and 
suppression  of  the  efficiency  limiting  processes.  The  nonlinear  susceptibilty 
calculations  and  the  efficiency  limitation  analysis  show  that  the  type  II 
resonance  enhancement  scheme  allows  efficient  multiline  frequency  up-conversion 
under  appropriate  conditions  in  Hg  or  HD.  These  molecules  have  extremely  small 
dispersion  and  absorption  in  the  infrared  and  have  relatively  high  peak  values 
of  on  resonance.  The  x^  frequency  dependence  also  allows  SFG  for  a 
number  of  different  pump  frequencies  when  one  or  more  pairs  of  frequencies  are 
two-photon  resonant. 

One  important  factor  to  consider  in  the  type  II  scheme  is  the  effect  of  com¬ 
peting  nonlinear  processes.  In  particular,  stimulated  rotational  Raman  scat¬ 
tering  can  be  a  dominant  frequency  conversion  process  under  certain  conditions. 

Several  reports  of  efficient  Raman  shifting  of  CO,  laser  radiation  in  hydrogen 

c  46  47 

indicate  the  significance  of  this  competing  effect.  ’ 

Two  important  aspects  of  stimulated  Raman  scattering  (SRS)  should  be  noted. 
First,  SRS  represents  a  loss  mechanism  for  the  pump  radiation  and  thereby  de¬ 
creases  the  THG  and  SFG  efficiency.  Second,  frequencies  shifted  by  SRS  may 
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participate  further  In  other  nonlinear  processes  Including  induced  two-photon 
absorption  and  sum  frequency  generation. 

A  survey  calculation  of  the  Raman-shifted  multiline  CO  frequencies  in  H2  and 

HO  shows  that  induced  TPA  and  SFG  may  indeed  occur.  Analysis  of  the  x^™  values 

for  these  frequencies  indicates  that  susceptibilities  in  the  neighborhood  of 
-34  6 

10  cm  /erg  are  possible  with  certain  combinations  of  the  pump  and  the  shifted 
frequencies.  Figure  9  shows  the  normalized  |x5rb|  as  a  function  of  two- 
photon  resonance  frequency  detuning  for  HD  Q(l)  transition  at  various  gas  den¬ 
sities.  Figure  10  illustrates  the  corresponding  curves  for  the  normalized 
1x1  as  a  function  of  gas  density  at  certain  detunings.  The  scale  factor 

for  the  actual  value  of  lxSF^l  is  3.14  x  10"34  cm6/erg.  Figure  11  shows  the 
K<»rr  TPA 

normalized  x  and  x  as  a  function  of  frequency  detuning,  which  are  ob¬ 
tained  from  the  real  and  imaginary  parts  according  to  the  equations  in  Table  9. 

The  actual  frequency  detunings  (Au)q)  obtainable  for  HD  resonances  from  the 
certain  combinations  of  the  pump  and  the  Raman  shifted  frequencies  in  H2 
(354.381  cm"1  shift29)  are  in  the  range  of  0.2  to  1  cm'1.  For  example,  P(ll) 
and  P(12)  of  [5-4] give  AwQ  =  0.328  cm"1  for  the  HD  0(0)  resonance, and  P(9)  and 

P(  13)  of  [4-3]  yield  Aui  *  1.012  cm"1  for  the  HD  Q(l)  resonance.  P(ll)  [6-5] 

0  -1 
and  P(12)  [4-3]also  approach  the  Q(0)  resonance  with  Au>  =  0.291  cm  .  These 
rrp  -34  6  " 

detunings  lead  to  x  values  of  10"  cm  /erg.  This  suggests  that  efficient 
SFG  may  occur  if  phase-matching  requirements  can  be  satisfied. 

5.2  Two-Step  Frequency  Conversion 

The  results  of  the  multiline  THG  and  SFG  analysis  suggest  that  an  efficient 
multiline  CO  laser  frequency  up-converter  may  be  possible  with  the  use  of  a  Raman 
shifting  as  an  intermediate  step.  The  Raman  shifting  in  H2  converts  a  part  of 
the  total  pump  radiation  into  coherent  radiation  near  6  iim.  The  frequencies  of 
the  remaining  pump  radiation  and  the  shifted  frequencies  add  to  yield 
the  two-photon  resonances  in  HD.  This  two-step  process  is  schematically 
Illustrated  in  Figure  12.  The  pump  ana  the  shifted  frequencies  are  repre¬ 
sented  by  {iii}and } ,  respectively.  The  generated  sum  frequencies  are  of  the 
form  uus  *  2w+  m  or  2u>  +  <a  and  are  denoted  by  ^  o» s } . 

1 
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FIGURE  9.  NORMALIZED  jxSFG|VS  HD  Q(l)  TWO-PHOTON  RESONANCE  FREQUENCY  DETUNING 


1 


6.0  EXPERIMENTAL  PROGRAM 

The  results  of  this  study  indicate  that  frequency  up-conversion  according  to 
the  type  II  resonance  enhancement  scheme  in  HD  may  be  potentially  efficient 
and  warrants  further  investigation.  Experimental  data  are  needed  to  verify 
theoretical  predictions  presented  in  this  report.  Of  particular  interest  is 
the  frequency  dependence  of  the  conversion  efficiency  and  of  the  phase-matching 
properties  of  HD.  This  may  be  used  to  determine  the  optimum  parameters  of  the 
nonlinear  medium  for  multi  frequency  CO  laser  frequency  up-conversion. 

6.1  Approach 

The  third-order  susceptibility  calculations  and  efficiency  limitation  analysis 
indicate  that  two-photon  resonantly  enhanced  sum  frequency  generation  in  HD  is 
the  most  promising  approach  to  up-convert  multiline  CO  laser  frequencies.  In 
order  to  verify  this  feasibility,  an  experimental  investigation  of  SFG  in  HD  is 
suggested.  The  investigation  should  consist  of  first  setting  up  an  experimental 
facility  to  generate  tunable  and  fixed  IR  frequencies  in  the  5-6  **m  region  and 
then  performing  parametric  SFG  experiments  using  the  generated  frequencies. 

These  IR  frequencies  are  relevant  to  the  two-step  frequency  conversion  process 
described  in  Section  5.2. 

6.2  Experimental  Facility 

An  IR  source  facility  based  on  a  Nd:YAG  laser  system  operating  at  ten  pulses 
per  second  would  be  a  suitable  choice  because  of  high  data  acquisition  rate. 
Figure  13  shows  a  block  diagram  of  the  experimental  setup  including  a  fixed 
and  a  tunable  IR  source.  The  fixed  IR  frequency  at  6.2  can  be  obtained  by 
vibrational  Raman  frequency  shifting  of  the  Nd:YAG  laser  output  in  H2  and  HD. 
Figure  14  shows  a  schematic  layout  of  the  fixed  IR  source.  The  tunable  IR  fre¬ 
quency  can  be  generated  by  difference  frequency  mixing5^  of  the  second  harmonic 
of  the  Nd:YAG  laser  output  and  a  tunable  dye  laser  output  near  0.596  t-m. 
Alternatively,  a  dye  output  at  0.690  ^m  may  be  Raman  shifted  three  times  in  H2 
to  yield  a  tunable  IR  output.  Figure  15  shows  the  two  possible  configurations 
of  the  tunable  IR  source.  The  tunable  IR  frequency  is  near  2015  cm"1  (4.96  rm) 
so  that  the  sum  of  the  tunable  and  the  fixed  IR  frequencies  matches  the  two- 
photon  Q-branch  resonances  in  HD.  The  operating  conditions  of  the  IR  source 
facility  are  summarized  in  Table  25. 
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FIGURE  13  SUM  FREQUENCY  GENERATION  EXPERIMENT  BLOCK  DIAGRAM 


FIGURE  14.  FIXED  IR  SOURCE  BLOCK  DIAGRAM 
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TABLE  25.  COHERENT  IR  SOURCE  FACILITY 


FIXED  IR  SOURCE 


Pump  Laser: 

Nd:YAG  1.064  p.m.  10  pps,  50  MW 
peak  power 

Frequency 

Converter: 

Vibrational  Raman  shifting  in  and 

HD  (7783.4  cm'1) 

IR  Output: 

~  5  MW  peak  power  at  6.2  pm 

TUNABLE  IR  SOURCE 

Pump  Laser: 

Nd:YAG  Second  Harmonic  0.532  p,m, 

10  pps,  20  MW  peak  power 

Frequency 

Tuning: 

Dye  laser  oscillator/amplifier  ~  7  MW 
peak  power  output  at  0.59  p.m 

Frequency 

Conversion: 

Difference  frequency  mixing  of  0.532  p.m 
and  0.596  p.m  in  LilO-j 

IR  Output: 

~0.02  MW  peak  power  at  4.96  p.m 
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An  apparatus  for  the  SF6  experiment  is  shown  in  Figure  16.  It  incorporates 
a  cell  consisting  of  a  tube  filled  with  several  atmospheres  of  HD  at  293  K, 
two  calcium  fluoride  windows (to  seal  the  ends  of  the  tube) which  transmit  the 
pump  and  the  generated  radiation,  and  a  gas  handling  system  consisting  of 
valves,  gauges,  and  a  vacuum  system  that  controls  the  cell  pressure. 


VACUUM 

SYSTEM 


FIGURE  16.  SCHEMATIC  OF  EXPERIMENTAL  APPARATUS 


6.3  Parametric  Study 

The  parametric  study  should  include  conversion  efficiency  measurements  as  a 
function  of  two-photon  resonance  frequency  detuning,  SFG  gas  pressure,  and 
pump  source  intensity.  From  these  measurements  the  SFG  nonlinear  suscepti¬ 
bilities  and  the  phase-matching  properties  of  the  nonlinear  medium  can  be 
compared  with  the  calculated  values. 


7.0  CONCLUSION 

Detailed  analysis  and  calculations  on  the  feasibility  of  efficient  CO  laser 
frequency  tripling  have  been  completed.  The  results  of  this  study  indicate 
that  direct  frequency  tripling  of  multiline  CO  laser  output  leads  to  small 
conversion  efficiencies.  These  efficiencies  are  primarily  dictated  by  the 
small  nonlinear  susceptibilities  and  the  efficiency  limiting  processes  of 
the  molecules  surveyed.  However,  the  analysis  of  the  three  resonance  en¬ 
hancement  schemes  has  yielded  a  potentially  efficient  sum  frequency  genera¬ 
tion  scheme  in  hydrogen  molecules  if  the  pump  frequencies  can  satisfy  the  two- 
photon  resonance  condition. 

Of  the  three  resonance  enhancement  schemes  investigated,  the  two-photon  reso¬ 
nance  (type  II)  scheme  has  been  found  to  meet  the  requirements  for  efficient 
multiline  frequency  up-conversion  of  the  CO  laser  frequencies.  Various  effi¬ 
ciency  limiting  processes  found  in  the  two  other  schemes  are  minimized  in  this 
scheme,  and  high  SFG  conversion  efficiencies  are  predicted  for  a  number  of 
pump  frequency  combinations. 

Unfortunately,  the  efficient  CO  laser  lines  do  not  yield  frequency  combina¬ 
tions  that  directly  satisfy  the  strong  resonance  conditions  in  H2  or  HD.  How¬ 
ever,  from  the  analysis  of  stimulated  rotational  Raman  scattering  as  a  compet¬ 
ing  efficiency  limiting  process,  the  Raman  shift  in  H2  has  been  found  to 
generate  frequencies  which  lead  to  HD  two-photon  resonances  when  combined  with 
the  original  pump  frequencies.  Because  Raman  shifting  has  been  demonstrated 
to  be  an  efficient  conversion  process,  this  resonance  coincidence  provides  a 
method  of  frequency  up-conversion  using  an  intermediate  frequency  shifting 
stage  to  optimize  the  sum  frequency  generation  efficiencies.  With  such  a  two- 
step  conversion  process,  output  wavelengths  in  the  1.77  and  1.90  p.m  regions 
can  be  obtained. 

In  order  to  verify  the  feasibility  of  efficient  sum  frequency  generation  in  HD, 
and  experimental  investigation  of  third-order  susceptibility  is  recommended 
for  future  study.  The  measurement  of  SFG  susceptibilities  and  phase-matching 
properties  of  HD  should  yield  useful  parameters  for  evaluating  this  molecular 
medium  as  a  potentially  efficient  frequency  up-converter. 
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APPENDIX  A:  x^  CODE 
Program  Structure  and  Performance 

The  program  consists  of  loops  over  the  vibrational  and  rotational  quantum 
numbers.  Inside  these  nested  loops,  subroutines  are  called  and  the  following 
steps  are  performed: 

(1)  Calculation  of  rotational  and  vibrational  terms,  and  partition 
functions 

(2)  Calculation  of  rotational  populations  and  all  transition 
frequencies  allowed  by  selection  rules 

(3)  Calculation  of  the  frequency  denominator  which  occurs  in 
the  nested  sum  over  states 

(4)  Inclusion  of  all  permutations  of  the  three  input  frequencies 

Several  subroutines  are  called  in  the  course  of  the  program.  The  structure 
and  function  of  each  of  these  is  relatively  simple,  and,  therefore,  will  only 
be  described  briefly.  Subroutine  FREQ  calculates  the  frequency  denominator 
and  sets  up  the  indices  to  be  printed  when  a  resonance  is  encountered.  Sub¬ 
routine  SUMO  provides  the  contributions  to  the  sum  over  states  of  the  total 
angular  momentum  projected  along  the  molecular  axis  (OMEGA)  when  the  molecular 
ground  state  is  spatially  degenerate  (LAMBDA  is  not  zero). 

Matrix  elements  involving  anharmonic  oscillator  eigenfunctions  which  are 

(3) 

required  for  inclusion  of  an  upper  electronic  level  in  the  x  calculation 
are  calculated  in  the  subroutine  LMENTS.  The  sum  over  state  of  the  total 
angular  momentum  projected  along  the  molecular  axis  for  a  spatially  degenerate 
upper  electronic  level  is  performed  in  the  subroutine  TUMO.  Finally,  the  fre¬ 
quency  denominators  which  involve  more  than  one  electronic  level  are  evaluated 
in  subroutine  GREQ. 

The  data  items  to  be  input,  specified  in  READ  statements  at  the  beginning  of 
the  main  program,  consist  of  alphanumeric  symbols,  indicators,  spectroscopic 
constants,  and  transition  moments.  The  spectroscopic  constants  may  be  found 
in  the  literature.  All  input  variables  should  be  expressed  in  cgs  units  unless 
otherwise  specified.  In  addition,  the  program  presets  all  unspecified  variables 
to  zero. 
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In  the  program,  the  input  variables  have  been  assigned  the  following  names: 

LABEL  is  an  alphanumeric  array  which  stores  the  name  of  the  molecule  for 
which  the  calculation  is  performed 

NV  is  the  number  of  vibrational  levels  used  in  the  calculation 

NJ  is  the  number  of  rotational  levels  used  in  the  calculation 

WE,  WEXE, 

WEYE,  WEZE  are  the  familiar  vibrational  constants: uu  ,  uu  x  ,  w  v  ,  u>  z  , 

e  e  e  er*  e  e 

respectively 

3E,  ALPHA, 

GAIWA,  OE, 

BETA  are  the  usual  rotational  constants:  B  .  a  ,  y  ,  D  ,  3  ,  respectively 

6  6  6  6  6 

AE,  DELTA  are  additional  rotational  constants  which  must  be  specified  only 
if  the  ground  electronic  state  is  spatially  degenerate 

OMEGAE  is  the  energy  of  an  upper  electronic  state;  this  is  set  to  zero 
if  it  is  desired  to  only  include  the  ground  electronic  state  in 
the  calculation 

MGE  is  a  real  array  which  stores  the  electronic  transition  moments 

U(5,5)  is  a  real  array  which  stores  the  dipole  moments  for  vibrational 

transitions  between  states  labeled  N  and  M  (in  Debye) 

FWHM  is  a  real  array  which  stores  the  linewidths  for  vibrational 

and  rotational  transitions 

TEMP  is  the  temperature  of  the  system  being  modeled 

WSIZE  is  the  size  of  the  interval  about  the  input  pump  frequencies 

inside  of  which  the  program  will  indicate  that  a  resonance  has 
been  approached. 

IWRITE  is  an  integer  flag  such  that  if  it  is  set  equal  to  one  the 

program  prints  out  the  individual  excitation  path  contributions 
or  if  it  is  set  equal  to  two  the  program  calculates  only  the 
contribution  from  the  va  =  0,  vb  =  3,  vc  =  2,  vd  *  1  excitation 
path 

IPOINT  is  an  integer  flag  such  that  if  it  is  set  equal  to  one  the  program 
calculates  results  for  frequency  tripling  for  pump  frequencies 
ranging  from  WQ  to  WP  with  a  step  size  of  WR 
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WP,  WQ,  WR  these  are  the  pump  frequencies  (in  cm"*).  The  program  terminates 
when  all  three  are  set  to  zero 

There  are  also  several  COMMON  variables  used  in  the  program.  SIZE,  AE,  DELTA, 
WE,  WEXE  are  also  input  variables  and  have  been  described.  Some  other  COMMON 
variables  are  listed  and  described  below: 

NA,  NB, 

NC,  ND  are  the  indices  for  the  appropriate  vibrational  levels 

RA,  RB 

RC,  RD  are  the  indices  for  the  appropriate  rotational  levels 

FABCD(20)  is  the  frequency  denominator 
EAB,  EAC.EAD  are  the  imaginary  half-widths 

W(2,5,50)  is  a  real  array  which  stores  the  calculated  energy  levels 
CA(5,50) 

CB(5,50)  are  the  intermediate  coupling  coefficients  (A  ±  0) 

RS(2,3,50)  is  a  real  array  which  stores  the  calculated  intermediate 
coupling  coefficients  (A  ^  0) 

X(5,2,3)  is  a  real  array  containing  the  anharmonic  oscillator  matrix 
elements  (when  an  upper  electronic  state  is  included  in  the 
calculation) 

OMEGAA.OMEGAB 

OMEGAC.OMEGAD  are  the  projections  of  the  total  angular  moments  on  the 
molecular  axis  (A  f  0) 

The  output  consists  of  various  blocks  listing  the  following: 

(1)  All  input  data 

(3)  Single  path  contributions  of  the  susceptibility  (if  IWRITE  =  1) 

(31 

(4)  The  sum  frequency,  the  real  part  of  x  ,  the  imaginary  part 
of  x(3),  and  the  modulus  of  x^ 

A  complete  calculation  for  one  molecule,  considering  up  to  4  vibrational  levels 
and  30  rotational  levels  takes  less  than  26  seconds  on  a  CDC  6600,  yielding 
results  for  10  sets  of  pump  frequencies. 
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CHITHRE 
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CHI THRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 

1?7 

pEA|)(5,10b>  AE.OELTa 

CHITHRE 
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CHITHRE 
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CHIThRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 

133 
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CHI THRE 

134 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 

150 

l«  (L  )  rc)i  -  ALPHA -»T  .GAUM  A®  T  ®»2 
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CHITHRE 
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CHITHRE 
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CHITHPE 
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CHITHRE 
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CHITHRE 
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CHITHPE 
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CHITHPE 
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WHITEI6.I22) 

CHITHRE 
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CHITHPE 
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CHITHPE 
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CHITHPE 
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CHJTMPE 

245 
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CHITHPE 

246 
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CHITHPE 
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CHITHPE 

248 

r,o  TO  ?3 

CHITHPE 

249 
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CHITHPF 
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CHITHPE 
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CHITHPE 
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CHITHPE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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IF(IPOINT.NE.I)  GO  TO  55 

CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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continue 
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CHITHRE 
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CHITHRE 
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CHITHPE 
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CHITHPE 
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CHITHRE 

293 

CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 

304 

(.u  TO  4 

CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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313 

1 

CHITHRE 

314 

*2=*:} 

CHITHRE 

335 

n3*«R 

CHITHRE 
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CHITHRE 
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CHITHRE 
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CMITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 
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CHITHRE 

375 

R=2»J-1 

CHITHRE 

376 

R="</2. 

CHITHRE 

377 

SA  =  B 

CHITHPE 

378 

pO  1 035  0=1. IN 

CHITHRE 

379 

PR  (<)=(). 

CHITHRE 

3«0 

Fl*.l)  =  <0.*0.) 

CHITHRE 

38 1 

lom 

F (K.2I a(0. >0.) 

CHITHRE 

3«2 

c 

CHITHRE 

383 

c 

SUM  0VfH  ROTATIONAL 

LEVELS  BY 

PATHS 

CHITHRE 

384 

c 

CMITHRE 

385 

lMc-1.5)  1015.1012. 

1011 

CMITHRE 

386 

101 1 

kN*1 

CHITHRE 

387 

5M*R-t. 

CHITHRE 

388 

SC=o-2. 

CHITHRF 

389 

SO=o-l . 

CHITHRE 

390 

P«(  1  )  =,’./  <  15. »  (2.»R- 

1  .  )  ) 

CHITHRE 

391 

C*LL  Sumo 

CHITHRE 

39? 

101? 

KM*® 

CHITHPE 

393 

S6*R-l. 

CHITHRE 

394 

sl=h 

CHITHRE 

305 

SU=P-1 . 

CHI THRE 

396 

PR(?)  =  (4.*m'*R.|.)/(I5.*R*(2.*P-I.)®<2.*R.1.I  ) 

CHI THRE 

397 

CaLL  S'lMO 

CHITHRE 

348 

KN=  « 

CHITHRE 

399 

5«(*o.  1  . 

CHITHPE 

400 

92 
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SC*P 

ChITmRE 

401 

so«r-i. 

chithre 

40? 

P«<3)  *?./  Ui.«  (2.«R»1  .)  1 

CMITMRE 

403 

call  suho 

CMITHRE 

404 

cmitmre 

405 

ChITmRE 

406 

tiC*V 

ChITmRE 

407 

Sl>*p.l. 

CHITMRE 

400 

P'<  t..  )  :(-•<  <31 

CHITHRE 

409 

1*Ll  b(|HO 

cmitmre 

410 

KH*  t 

cmitmre 

4|  | 

Sr'=t.' 

cmithre 

412 

S 

CHITMRE 

413 

SH=9-1 . 

CMITMRE 

414 

hX17)=(K-1.)/(1S.,H*(2.««»1.I 1 

cmitmre 

415 

C*Ll  Sl)MO 

cmitmre 

416 

kn=9 

CHITHRE 

417 

SM=R 

CMIThRE 

410 

CHITHRE 

419 

SU*u 

CHITMRE 

4?0 

PK(4)*MM(7> 

cmithre 

421 

OLl  SUMO 

cm i Tmre 

422 

1  1 

CHITHRE 

423 

S<J=»—  1  . 

cmitmre 

424 

SC=H 

CMITHRE 

425 

-,U»R 

CHITMRE 

426 

1 l ) =RR»7i 

cmithre 

427 

CaLl  SliMO 

Cmitmre 

428 

K«*13 

cmitmre 

429 

S»*«-l . 

CHITHRE 

430 

5L»R-1. 

CMITHRE 

431 

SU=«-1 . 

CMITMRE 

432 

RR  (13)  /<IS.*R»  <2.  •«-!.*» 

cmitmre 

433 

call  Sumo 

cmitmre 

4  34 

K"=15 

cmithre 

435 

SO=R 

cmithre 

436 

SC=H-1 . 

cmithre 

437 

SU*R-I. 

cmitmre 

438 

pRI 15)  sSLriT  (  lo.l .)• <R-1.)/ <R«R«  <2.»R*1  .)*  (2.»H-l  .)  >  »  /■  1 5  . 

cmithre 

439 

call  sumo 

CMITMRE 

44  0 

km=1  7 

cmitmre 

44  1 

Sn=M- 1 . 

Cmitmre 

442 

SC*R-l . 

cmitmre 

443 

bU*R 

cmitmre 

444 

pR<l  /laRHtlb) 

cmitmre 

445 

C»Ll  sumo 

Cm  I ihRE 

446 

rips 

cmitmre 

447 

b‘i=».l. 

cmithre 

448 

bC».) 

Cmitmre 

449 

SU*o. 1 . 

cmitmre 

450 

p«(S)  =  (<..*K*H.rt.«H.5.)/(lb.#tR.l.l»(2.,H»l.),(2.*R*3.l) 

CMIThRE 

45 1 

CALL  SllMO 

cmitmre 

452 

KU  =  ft 

CMITMRE 

453 

Su*W.l . 

CMITHRE 

454 

5«C  =  w.2. 

CMITMRE 

455 

so=p.i. 

CMITHRE 

456 

)"<<6)  »?./<  15.»C2.*H*3.)  ) 

CHITMRE 

457 

.TOO 


CHITMRt 
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lOtO 


CALL  SllMO 

CHI  TMPC 

450 

*N*0 

CHITHRE 

459 

SH*R 

CHITHRE 

460 

st*p 

CHITHRE 

46 1 

Sl>*P*l  . 

CHITMRC 

462 

PW»«)=(M.2.)/I15..(rt.l. )*(?.•».!.)) 

CHITMRC 

46  1 

c*ll  sumo 

CM I TMPC 

464 

ku*io 

CM I  TMPC 

465 

Stt*R 

CHITMRC 

466 

sc*p.i . 

CHITMRC 

46  7 

CM 1  TMPC 

460 

PH ( 10) =PR (8| 

CMITMPC 

469 

C“LL  SUMO 

CHITMRC 

470 

KN=i,3 

CMITMPC 

*71 

SU*R»l. 

CHITMRC 

472 

5C*o 

CHITMRC 

473 

SO*R 

CHITMRC 

*74 

pH  < 121 sPR («) 

CHITMRC 

475 

call  SijMO 

CHITHRE 

476 

kN=14 

CHITMRC 

477 

Sa3R*I . 

CHITHRE 

470 

SC*R.l. 

CHITHRE 

479 

SOaP.l . 

CHITHRE 

400 

MH(l4)a«/(l5.'MR.l.),(2.*M.3.)  ) 

CHITHRE 

40  1 

call  sumo 

CMITMPC 

40? 

KN*16 

CHITHRE 

403 

SO*R 

CHITHRE 

404 

SC*P.l. 

CHITHRE 

4H5 

so*n»i. 

CHITHRE 

406 

PHI  16 >*SOKT(R»<R. 2.1/1  (P.l.)«(p.l.)»(2.»R«l  .)*(2.*R»3.)  )  >/15. 

CHITHRE 

407 

CALL  sumo 

CHITHRE 

408 

KN*ia 

CHITHRE 

489 

sa=R*l . 

CHITHRE 

490 

SC3R  » l . 

CHITHPE 

491 

SO*R 

CHITHPE 

49? 

PR(iai=PR<i6) 

CHITHRE 

493 

call  Sumo 

CHITHRE 

494 

KU*19 

CHITHRE 

495 

5b*o 

CHITHRE 

496 

SC*o 

CHITHRE 

497 

S03R 

CHITHPE 

490 

PR  < 19) * (3.ARAR. 3,#H“1 . ) / ( l5.»R» (R*l.)*(2.»R.l.) ) 

CHITHPE 

499 

call  sumo 

CHITHRE 

500 

CHITHRE 

501 

HtOlM  SOM  UvE°  PATmS 

CHITHRE 

502 

CHITHRE 

50  3 

PR  1  12)  *R6  (01 

CHITHRE 

504 

0U  1040  6*1 .19 

CHITHPE 

505 

SU*L  ( 2 )  *  <  0  .  .  0 . 1 

CHITHRE 

506 

$UML  ( 1 )  «  ( 0  .  .  0  .  ) 

CHITHRE 

50  7 

SUML ( 1 ) *SUML ( 1  I .PR (6) *F <K. 1 ) 

CHITHRE 

500 

IF  IP. EO. 0.3)  r,o  TO  1040 

CHITMPE 

509 

SOML(2)«SUML<2).P«tAI»f (A. 2> 

CHITHRE 

510 

COM  nut 

CHITHRE 

511 

RR* (R.fl.S) 

CHITMPE 

512 

1 3«V  1 1 . 1 ) •KW  —  IV (1.1  )»RW»»2.HV( 1 .1 >»RR*» J-RV (1*11  * OV « 1*1  )-HV( 1 • 1  ) 

CHITHRE 

513 

OlS*t<R<-t/  16. 95031.-0  !• TEMP)  I/O 

CHITHRE 

514 

94 
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ctirrn*t 
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iUHj=SuHJ«^KMLU  1*1)1  i 

E*bV(i?2ll*HH-0tfM?i»t*HK»*2.MV<1.2>»HH*«J*ASC»li,<l)-HW«l.l)»0W<l*ll 

i-MV  ( 1  •  1 ) 

l)IS=f  *P<-t/(6.4S0Jt-01«T£MP>  I/O 
SUMj=SllMJ*3UML  (2)  *DlS 

105(1  C.Vt'l  !M()t 

0‘J  TO  lilt 
10M1  CONTINUE 

po  so  J*  1  »)*J 

;  H  Is  TmE  ACTUAL  ROTATIONAL  OUAMTUM  NUHBtR 


R3J-I 

RA*R 

f>U  15  6*l'« 
T>K|».)*0. 

-\1-  F  Abr.D  CK  )  -  <  0  .  •  0 . 1 


su«  OVe»  IrtE  ROTAUONAL  levels  by  paths 


jF<k-l  .) 15*12.1 1 
I  I  K*» 

PH*»-1. 

Rt=*p-2. 

1  # 

P«(l  )  =2. •H'»(R-i.)/|lS.*(2.*«-l.>*<?. *«♦»•) ) 

CALL  FhEU 
1?  K** 

P'J=P-1 . 
pL*J 
PU=W- 1 , 

p«  ( 2>  *k* ( a.*H»«2» 1 » / ( IS.* (2.*R”l • ) * (2»*W» 1 . 1 «*2) 

CALL  frfu 
K=J 

RIJAR.l . 

pc-p 

HU*P-1. 

pA  (3)  *2.»K**  <R«1.)/(1S.*<2.*R»1.>**2> 

CALL  FREU 

r  -a 

kA*P-l . 

kC*D 

RU=R.l. 

pK<kl=?.*A»(M.l.|/(15.*(2.*H*l.)**2) 

C»tL  FKEO 
IS  K*S 

PB*H»1 . 
pL  =  l> 

pAiSls(K.l.l*(A.»B»*2*tt.#H»S)/(15.*(2.*H»3.)»(2.*R«l.)**2) 

CALL  fkEu 
r  *o 

k'jAp  ♦  1  . 

rl*«.2. 

pllAU.  1  . 

=?.»(«.  I  2.)/(  1S.*(2.»H»3.1*  t2.*R*|.l  I 

C*Ll  FkEU 


CmITHRE 

SIS 

CHITHRE 

SI6 

chithre 

SI7 

ChlTHRE 

Sis 

ChIThRE 

SIR 

CmIThRE 

S?0 

CHIThRE 

S2I 

ChIThRE 

522 

CHIThRE 

S23 

ChITmrE 

524 

CHITHRE 

S2S 

CHIThRE 

S26 

CHITHRE 

527 

CHITHRE 

S2B 

CHIThRE 

529 

CHITHRE 

530 

CHITHRE 

531 

CHITHRE 

5.12 

CHITHRE 

533 

CHITHRE 

534 

CHITHRE 

535 

CHITHRE 

S36 

CHITHRE 

537 

CHIThRE 

538 

CHITHRE 

539 

CHITHRE 

540 

CHITHRE 

541 

CHITHRE 

542 

CHITHRE 

54  3 

CHITHRE 

544 

CHITHRE 

545 

CHITHRE 

546 

CHITMRE 

54  7 

CHITHRE 

548 

CHITHRE 

549 

CHI THRE 

550 

CHITHRE 

551 

CHITHRE 

552 

CHITHRE 

553 

CHITHRE 

554 

CHITHRE 

555 

CHITHRE 

556 

CHITHRE 

557 

CHITHRE 

558 

CHIThRE 

659 

CHIThRE 

560 

ChIThhE 

561 

CHITHRE 

562 

CHIThRE 

563 

CHITMRE 

564 

CHIThRE 

S65 

CHITHRE 

566 

CHIThRE 

567 

CHIThRE 

568 

CHITHRE 

569 

CHIThRE 

S70 

CHIThRE 

571 

non  n  ri  on 


CMTTHPt 
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ptGtn  the  so*  ovep  paths 


$umk  =  i  o  •  *  v .  > 

D<>  40  6  =  1*6 

4.0  SUMK3S|IM6*PR  0<(»FAHCD<K) 
fiK=p» (R»  I .  I 

t*6 ( 1 ) =  RH-0 ( 1 I •RR»»2*HE,RP,,3 
p IS =( 2.*R ♦ 1. 1 »E*P <-£/< 6.950 3E-0 l • TEMP) I/O 
SO  SOHj=S|lMJ*sUMS«0tS 

r.u  ro  1 1 1 1 

30nn  jNliF)(*l 

ll'<MC.E<*«NV‘l>  IN0E*=2 
wV*'J(NA»NU)/0«EGAE 

IjtOlN  THE  SUM  OVER  ROTATIONAL  LEVELS 
S'JMJ=(0.»V»> 

00  2050  J* 1 *NJ 

P  IS  The  ACTUAL  ROTATIONAL  QUANTUM  NUMBER 

R  =  J-i 

TA*R 

n‘J  ?03S  K  =  l,6 

rMt<l*0. 

r.<K, i)*(0. to. > 

207C  c,(K,2)  =  (0.tO.I 

SUM  OVER  MUTATIONAL  LEVELS  BT  paths 

tFIP.LT.n  GO  TO  2015 

kM*i 

Tti=P*l. 

TC*P 

TU*P-1. 

pH<ll*?./(l5.»t2.»RM.l) 

CALL  TllMO 
20 |S  kH*2 

Td*P.I . 

TC*P 
TO*R*l . 

pH(2)=<4.,H*R*8.*H*S.)/(15.,(R*l.)*(2.*RM)»(2.,R*3.>  ) 

CALL  TUMO 

rM*3 

TH*P*1 . 

tC=r*2 . 

TO*R«l. 

pK(1)3;i./U5.»(2.#P‘3.)  I 
C*LL  TllMO 

Km*4 

TH*P.l, 

tC*p 

TU*P 

pM(/,  )=-(Rt2.)/(lS.<>(Rtl.)*(2,»P«l.)) 
fA'-L  TllMO 
kM*s 
T  o  *9 ♦  1 . 

rc*«.  i . 

TO*w ♦  l  . 

pK<n)*-H/( 15.»(R‘ I .)»C2. •P*3.)| 

CALL  TllMO 


CHITHHE 

57? 

CHITHRE 

573 

CHlTHRE 

574 

CHITHRE 

575 

CHITHRE 

576 

CHITHRE 

577 

CHITHRE 

57B 

CHITHRE 

579 

CHITHRE 

580 

CHITHRE 

58) 

CHITHRE 

58? 

CHITHRE 

583 

CHITHRE 

584 

CHITHRE 

585 

CHITHRE 

586 

CHITHRE 

587 

CHITHRE 

688 

CHITHRE 

589 

CHITHRE 

590 

CHITHRE 

591 

CHITHRE 

592 

CHITHRE 

593 

CHITHRE 

59a 

CHITHRE 

595 

CHITHPE 

596 

CHITHRE 

597 

CHITHRE 

598 

CHITHRE 

500 

CHITHRE 

600 

CHITHRE 

o«l 

CHITHRF 

602 

CHITHRE 

603 

CHITHRE 

604 

CHITHRE 

60S 

CHITHRE 

606 

CHITHRE 

607 

CHITHRE 

608 

CHITHRE 

609 

CHITHRE 

610 

CHITHRE 

61  I 

CHITHRE 

612 

CHITHRE 

613 

CHITHRE 

614 

CHITHRE 

615 

CHITHRE 

616 

CHITHRE 

61  7 

CHITHRE 

618 

CHITHRE 

619 

CHITHRE 

620 

CHI  TMPC 

6?  1 

CHITHRE 

6?2 

CHITHRE 

6?3 

CHITHPE 

624 

CHITHRE 

6?5 

CHITHPE 

6?6 

CHITHRE 

627 

CHITHPE 

628 

96 
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11.36 

)<M®ft 

CHITHRE 

629 

T  ti  - 1 ♦  1  . 

CHITHRE 

630 

T  C  =  P  *  1  . 

CHITHRE 

631 

TO*l< 

CHITHRE 

632 

ftft(n)  =-S0<<  T(R®(H®2.)/((N®1.),(R®I.)#I2.®R.1.)®(2. 

•HO.)  )  )/15. 

CHITHRE 

633 

CAt-L  TUMO 

CHITHRE 

614 

C 

CHITHRE 

635 

C 

Hto jn  r«t  som  oven  maths 

CHITHRE 

636 

C 

CHITHRE 

637 

TUM(_i(O..U.) 

CHI THRE 

638 

1)0  2041)  A=1.6 

CHITHRE 

639 

2040 

TUML  =  fUML*PH (X  )  *G  <K  «  INDEX) 

CHITHRE 

640 

pR=R® (ft.  1 . ) 

CHITHRE 

641 

E  aa ( 1 ) »KR-U( 1 ) ®RH®®2.Ht®W«®,3 

CHITHRE 

642 

|>15=(2.®H»l.)  ®£XP(-£/  <o.9503E-01*TFMP>  )/Q 

CHITHRE 

04  3 

205p 

SOm  isSUMJ*  fOML“OlS 

CHITHRE 

644 

MX  sue 

CHITHRE 

645 

)1  (  INDEX. LU.2)  NX=NO 

CHITHRE 

646 

U‘J=MG£(l)®MGt(l)®X<NX.  index.  1)  *  IMGE  (  l  *  ®HGE  ( 2 )  *M6E  12)  ®MGE  ( 1 )  ) 

CHITHRE 

647 

=  ‘A <NX« INUtA,2l  .  (mo£ ( 1 ) ®MGE ( 3 ) »MGE 12) ®HGE ( 2 ) *M6E ( 3 ) ®MGE  < 1 ) ) 

CHITHRE 

648 

=»X(NX. (NOtX  « 3 ) 

CHITHRE 

649 

pv=pv*uo 

CHITHRE 

650 

ill  i 

CONTINUE 

CHITHRE 

651 

c 

f NO  THE  SUM  over  rotational  LEVELS 

CHITHRE 

652 

SOMV  (Nft.NL.NO)  =SllrtV  (NU.NC.NO)  «SUHJ®PV 

CHITHRE 

653 

qo 

COM  INI  it 

CHITHRE 

6S4 

C 

CHITHRE 

655 

c 

t Nl)  THE  Sum  OVE«  FREQUENCY  PERMUTATIONS 

CHITHRE 

656 

r 

CHITHRE 

657 

SOMV  CNM.NC.NO) =PM®<-2. 1 3E-26I »SUMV (N&tNCtNOI 

CHITHRE 

658 

|F( HR1 Tt.NE.l)  GO  TO  100 

CHITHRE 

659 

NUA=NU-1 

CHITHRE 

bftO 

NLA=NC-1 

CHITHRE 

Oft  1 

NUA=H0-1 

CHITHRE 

6ft2 

MUUSUMVSCAU5 (SUMV (Nb.NC.NO) ) 

CHITHRE 

6ft3 

WRITE (ft. H»)  NHA.NCA.NOA.SUMV (N8.NC.N0) .MOOSUMV 

CHITHRE 

6ft4 

lift 

FORMAT (22A. TS. 10X. IS.lOX. I5.8X.EIS.8.8X.E15.8.8X, 

E15.8) 

CHITHRE 

6ft5 

lOP 

CNl  3  =  ChI3«suMV(N8.I)C.NO) 

CHITHRE 

6ft6 

HO 

CONTINUE 

CHITHRE 

6ft  7 

70 

CONTINUE 

CHITHRE 

6ft8 

*0 

CONTINUE 

CHITHRE 

669 

c 

CHITHRE 

670 

c 

eng  the  sum  over  VIBRATIONAL  levels  for  A  GIVEN  PERMUTATION  OF 

THE 

CHITHRE 

671 

c 

PUMP  FftEUUtNCIES 

CHITHRE 

67  2 

c 

CHITHRE 

673 

c 

CHITHRE 

674 

c 

OuIpuT  Thl  RESULTS 

CHITHRE 

675 

c 

CHITHRE 

676 

WS=WP.ftU»WK 

CHITHRE 

677 

MOOCHIrCAUS(CHlJ) 

CMITMPE 

67ft 

WRITE (ft. 1  IS ) WS.CHl 3 'MOOCH  I 

ChitmRE 

6  7V 

1 1* 

F"RmaT(  /  i  5X »  *FOR  THE  SUM  FREflUENCT  ®«El5.8,/.5X, 

•THIRO  order 

SUS 

CHITHRE 

6ftO 

ECLPTIUlLlTr  HAS  RE  PART  *.£15. ft.®  CM6/E«0  ®,/.SX. 

•THIRD  ORDER 

SUS 

CHITHRE 

bftl 

;CtPTIBlLlTr  HAS  IM  PART  ®.ElS.A.®  CN6/EMG  ®,/.SX, 

•third  order 

SUS 

CHITHRE 

6ft2 

icEPriaiLlTY  HAS  MODULUS  ®»El5.ft.®  CM6/ER6  •) 

Crt[THH£ 

6«3 

2u<iO 

CUNT  iniiE 

CHITHRE 

6ft4 

GO  TO  25 

CHITHRE 

6ft5 

97 
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10nn 

W*lTE  I6«  U6> 

CHI  rw»E 

6«6 

116 

rUHMAT (////, 5*. *THaT*S  ALL  FOLKS* * 

CHITHWE 

sfop 

CH I THPE 

6Mfl 

f'JD 

CHITHRE 

b*9 

98 


frfo 

74//*  OPT*,?  UO  FTN  * ,b.4  33fl 

1 1/22/78 

11. 36, 

V'OPUOTINC  FREQ 

FREO 

2 

c 

FREO 

3 

c 

T«Is  SlIbNuUIt  JC  CALCULATES  Th£  FREQUENCY  DENOMINATOR  F ABCO 

FREQ 

4 

c 

FREO 

s 

INTEGER  o«CoAA.OMEgAB<UM£GAC  <Om£GAO 

FREO 

6 

CuHulEa  F  ABCO  <  20 ) <EAB<EAC<EAO<FD1  <F02<FD3<F04<  FQ5<F06<GABC0  (20 1 

FREO 

7 

f  uMmoN/CF «LU/  NA<NH<NC<NO<RA<RR<RC<RD<W1<W2<W3<K< W  <2<5<50 ) < 

FREQ 

8 

=F  AorDlpO ) <t Ad  <EAC<EAO<WSI2E<AE<OFLTA 

FREO 

9 

CUMmON/0SUMU/l,OMEGAA<OMEGAB<0m£G AC<OMEGAO<GABCO< 20) <SA<SB<SC<SO 

FREO 

10 

Cl 

OtbuG 

FREO 

11 

c* 

ARRAYS 

FREO 

12 

Cl 

CiUlOS 

FREO 

13 

IK<aE.EU.U..AMO.OElTA.EQ.O.)  GO  TO  S 

FREO 

14 

c 

FREO 

15 

c 

calculate  energy  level  differences 

FREO 

16 

/* 

FREO 

17 

ja*ga ♦  o  .5 

FREQ 

18 

JU*SH<0.5 

FREO 

19 
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SObPOUT INC  LMENTS 
COMnoH/LMuCOH / *E . Wt  XE « * ( 5  t  ? • 3 ) 
ktAt.  N(b) 

frtls  SUBROUTINE  calculates  The  expectation  value  of  x»«i  fob 

ANh/iWMONlt.  OSC  ILLATUB  E  IOEnFUNCTIOnSI  THE  RESULTS  ARE  STORED  IN 
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j=2  for  case  ? 
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jSURT (R»  (k* 1 . I  I  1/72.1 

i)L i £RMiNt  The  matrix  elements 
x  <  i .  i .  1 1 1 1 . 
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x(l,l.?)*-tPSIL»3./2./N«l>/NOl 
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Xl2,2«?>»U.»EPS|L**2.*< 17. "SORT  12. 1 ♦ 7a. *209. ‘SORT I2./3.) *14.» 
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X  <4 ,2, 3)  s-tPS IL*  ( 2. "SORT  ( J.  I  .9,»SURT  (2.  1  ^.‘SOrT  (6.  1 1 /2./N (4) /N ( 3) 
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HCL  THIMO  UHOE*  SU«re»U«lLlT’r  »C»N 


tHIS  CM.CUCAT10N  u*rs  »  1  l<l«*T  IO*»t  LEVELS  »nO  io  oof* r  I0N1L  LEVELS 


Vttt*Af lO^AL 

CONSTANTS t 

we 

.299o94oJc*04 

-EXE 

•528185006*02 

were 

.22430(1006*00 

S6ZE 

•121000006-01 

CN  RECIPROCAL 

CMI 

ROTaT  lUNAL 

CONSTANTS  I 

Bt 

.  lf)5**34l6C*a2 

ALPHA 

.307181006*00 

GAMMA 

•17724qOOE-02 

06 

•531936006-03 

S6TA 

.751000006-05 

( IN  RECIPROCAL  Cm» 

ROTATIONAL 

constants: 

1) 

A€  OttTA 

.  0. 

(IN  RECIPROCAL  CM) 

electronic 

STATE! 

EvC*‘»» 

.flOOnOOaUfOS 

NOCUl 
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NGC  (2) 
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HOC  (3! 

-.775000006-01 

I.INCKIOTH*:  ROTATTOnAL  0  TO  1  0  TO  2  a  TO  3  a  TO  4  (IN  RECIPROCAL  CHI 

.22SoonnoE*oi  ,24d«oouo£*oo  .zisooOoof.oo  ./osoooooe-oo 


tfrre r *  ture  row  this  calculation  is  .300oooooe*o3  kelv 


resonances  « i thin  a  range  of  .  ioooooooe-oo  /cm  centered  about  each  of  the  pump  frequencies  are  inoicaifo  bt  the  program 


THE  PUMA*  FREQUENCIES  ARE  . 20 2S0680C -04  .202S06doE *04  .202S0680E*04 
FOR  THE  SUM  FREQUENCY  . 40 752U40E -04 
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PFSONAHCtS  4 1 TmIn  i  u»n6£  OF  . I OOOOIOOt *00  /Cm  CEnIEREO  ABUUT  EaCh  OF  THE  PUMP  FREQUENCIES  ARE  |Nn ICATEO  BY  THE  PROGRAM 


Th£  pump  FREQUENCIES  A«t  .  ISNQlAlot'Oo  . lBGOlklot'Ok  .14n2*SOE»Ok 
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Third  OnU£«  SUSCEPTIBILITY  has  ROOULUS  .13SB85B2E-37  CH6/ERG 


RESONANCES  RIIHIN  a  range  OF  . lOOOOOOOfOO  /CM  centered  ABOUT  Each  OF  the  PUMP  FREQUENCIES  ARC  indicated  by  TmE  Program 


the  PURR  F»EOUENCIFA  ARt  *10  7 32flsi)£  *1)4  . J 9702B50E “Ok  .202506B0E“0k 
FOR  The  SUR  Frequency  . 517 1 0 J40E “Ok 

third  -.hoer  susceptibility  has  re  part  i i«oS77oe-37  ch6/e«o 
third  JROER  SUSCEPTIBILITY  has  Ir  Part  .2skJ8SlS£-k2  CM6/ER0 
ThIGO  OroeR  SUSCEPTIBILITY  has  HOU01.US  .i1So5770E-3T  CR8/ER0 


PFSQNANCES  4IIMIN  a  oanoa  OF  . I 000O000E • 00  /CR  CENIERCO  ABOUT  EACH  OF  THE  PURP  FREQUENCIES  ARE  INDICATED  By  IhE  Program 


IhE  PURR  FhtauENCIFY  ARE  . I  97 32BSOE “Ok  . 19702B5QE “Ok  . 1 999 L » l OE “Ok 
FOR  Iha  SUR  FREQUENCY  . S9kS  71  UE  •  Ok 

THIRD  OROER  SUSCEPTIBILITY  HAS  RE  PART  - . 1 24 1 7k I3E-37  CM8/ER0 
TMIRU  JMUER  SUSCEPTIBILITY  MAS  IR  PART  .2k3S23»8E-k2  CH8/ER0 
third  qruer  susceptibility  mas  "uuulus  .I28i7k9jc-jr  chs/erg 


RFSUNANCts  H I  Thin  A  hauGa  of  .  1  0  G.jOOUOf  JO  /CR  CENIEREO  ABOUT  EACH  OF  THE  PUHP  FREQUENCIES  ARE  I  NO  I CATED  BY  Tm£  progRar 


rue  pu*p  fpeoucncip <  »«f  . aozsoaaot -o*  , i»*pi*io£*o*  . i»733asoc»#» 

rc«  THE  SUH  rHEOUEwrT  .sboI.'H.oE.OL 

rHIBO  a»ut«  SUSCEPTIBILITY  hPS  «E  »*»T  -.I04B7707E-37  CM6/CM6 

TMIPO  OPUEH  SUSCEPTIBILITY  m»S  IN  »»» T  .31 7B3KMC-L2  CH»/£«0 

TH1BO  UHOEN  SUSCEPTIBILITY  h*S  "OOULUS  .104B7707E-37  CN»/£NO 


TH»r«S  *LL  FOLKS 


APPENDIX  B:  x 

CODE 

Program  Structure  and  Performance 

This  program 
complex.  Its 

is  quite  similar  to  the  code.  It  is,  however,  much  less 

essential  steps  are: 

(1) 

Calculation  of  rotational  and  vibrational  terms,  and 
partition  functions 

(2) 

Calculation  of  rotational  populations  and  all  transition 
frequencies  allowed  by  selection  rules 

(3) 

Calculation  of  the  frequency  denominator  which  occurs  in 
the  sum  over  states 

(4) 

Calculation  of  CHIONE  accounting  for  population  in  various 
vibrational  states 

(5) 

Derivation  of  the  absorption  coefficient  and  the  refractive 
index  minus  unity  from  the  value  of  CHIONE 

(6) 

Evaluation  of  the  k- vector  mismatches  if  desired 

Two  subroutines  are  called  in  the  course  erf  the  program.  Subroutine  ENERGY 
calculates  the  required  energy  levels  and  partition  function.  Subroutine  CHI 
calculates  the  real  and  imaginary  parts  of  CHI  for  a  particular  lower  state 
vibrational  level. 


The  data  is  input  according  to  READ  statements  at  the  beginning  of  the  main 
program,  in  a  format  very  similar  to  that  used  by  the  y.  ^  code.  All  input 
variables  should  be  expressed  in  cgs  units  unless  otherwise  specified.  Un¬ 
specified  variables  will  be  preset  to  zero  by  the  program. 

The  input  variables  for  this  program  are  as  follows: 

LABEL  is  an  alphanumeric  array  which  stores  the  name  of  the  molecule, 
for  which  the  calculation  is  performed 

NV  is  the  number  of  vibrational  levels  used  in  the  calculation 

NJ  is  the  number  of  rotational  levels  used  in  the  calculation 

NOVT  is  the  number  of  overtones  used  in  the  calculation 
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t 


are  the  familiar  vibrational  constants:  w  ,  u>  x  ,  »  v  ,  i#  z  , 

e  e  e  or  e  e  e 

respectively 

are  the  usual  rotational  constants:  B  ,  a  ,  v  ,  D  ,  g  , 

"  C  c  c  C 

respectively 

is  a  real  array  which  contains  the  dipole  moments  for 
vibrational  transitions  between  states  labeled  N  and  M  (in  Debye) 

is  a  real  array  which  stores  the  linewidths  for  vibrational 
and  rotational  transitions 

is  a  real  array  consisting  of  the  initial  distribution  of 
population  among  the  vibrational  levels 

is  the  temperature  of  the  system  being  modeled  (in  Kelvin) 

are  pump  frequencies;  the  k-vector  mismatches  are  calculated 
between  the  pairs  of  frequencies 

is  the  population  transfer  between  adjacent  vibrational  levels 
(used  to  compute  the  k-vector  mismatch) 

There  are  three  labeled  COMMON  blocks  in  the  program. 

VIBCOM  contains  the  vibrational  constants  WE,  WEXE,  WEYE,  WEZE 

ROTCOM  contains  the  rotational  constants  BE,  ALPHA,  GAMMA,  DE 

BETA,  TCOM  contains  the  temperature  inputted  as  TEMP 

The  output  consists  of  various  blocks  listing  the  following: 


(1) 

All  the  input  data 

(2) 

The  input  frequencies  and  the  population 
adjacent  vibrational  levels 

transferred 

between 

(3) 

For  each  input  frequency,  the  absorption 
refractive  index  minus  unity 

coefficient 

and  the 

(4) 

The  k-vector  mismatch 

A  complete  calculation  for  one  molecule,  considering  1  vibrational  level,  30 
rotational  levels,  and  1  overtone  takes  less  than  a  second  on  a  CDC  6600,  yield 
ing  results  for  42  pairs  of  input  frequencies. 


WE,  WEXE, 
WEYE,  WEZE 


BE,  ALPHA, 
GAMMA,  DE, 
BETA 


U(5,5) 


LINEW 


POPU 

TEMP 
WS,  WP 

DN 
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Program  Listing 


74/74  QPT*2  FTN  4,6*4338 

12/01/78 

OR. 32 

PROGRAM  CHIONE ( INPUT • OUTPUT • T  APES* INPUT  %  T  APE6*0UTPUT 1 

CHIONE 

2 

c 

CHI ONE 

3 

c 

THIS  PROGRAM  CALCULATES  FIRST  OROER  SUSCEPTIBILITIES  and 

RELATED 

CHIONE 

4 

c 

PROPERTIES  from  SPECTROSCOPIC  HAT A 

CHIONE 

s 

c 

CHIONE 

6 

c 

CHIONE 

7 

c 

CHIONE 

8 

c 

INPUTS 

CHIONE 

R 

c 

CHIONE 

10 

c 

CHIONF 

11 

c 

La«EL 1 1>  *1*1 .8 

( BA  1 0  I 

CHIONE 

12 

c 

CHIONE 

13 

c 

CHIONE 

1* 

c 

NV«HJ.nOVT 

(115) 

CHIONE 

15 

c 

CHIONE 

16 

c 

NV  =  Nl |MR£W  of  VIBRATIONAL  LEVELS  USEO  IN  THIS  CALCULATION 

CHIONE 

17 

c 

nJ  =  number  of  rotational  levels  used  in  this  calculation 

CHIONE 

l« 

c 

nuvt  =  number  of  overtones  useo  in  this  calculation 

CHIONE 

1R 

c 

CHIONE 

20 

c 

CHIONE 

21 

c 

w£.wexf»wete.weze 

(4E15.B) 

CHIONE 

22 

c 

CHIONE 

23 

c 

THESE  ARE  The  VIBRATIONAL  CONSTANTS  (IN  RECIPROCAL  CM) 

CHIONE 

24 

c 

CHIONE 

25 

c 

CHIONE 

26 

c 

be «A^"hA, uAMMA.OE. bet a 

(SEIS.8) 

CHIONE 

27 

c 

CHIONE 

28 

c 

THt-SE  ape  The  ROTATIONAL  CONSTANTS  (IN  RECIPROCAL  CM) 

CHIONE 

2R 

c 

CHIONF 

in 

c 

CHIONE 

31 

c 

FOR  each  VlHRATIONAL  level,  n - 

CHIONF 

32 

c 

CHIONE 

33 

c 

U<M,N>  .  Mai .NOV  T  .2 

(5E15.B) 

CHIONE 

14 

c 

CHIUNE 

15 

c 

ll(H.N)  3  UIPOLE  MOMENTS  for  VIBRATIONAL  TRANSITIONS  N  TO 

M 

CHIONE 

36 

c 

( in  oebtEi 

CHIONE 

37 

c 

CHIONE 

38 

c 

CHIONE 

1R 

c 

L  INFVI  (M»N)  •  M  =  1  .NOv  T  *  2 

(5E15.B) 

CHIONE 

40 

c 

CHIONE 

4  1 

c 

LlNEWInTHS  (IN  RECIPROCAL  CM)  FOR  ALL  OVERTONE  TRANSITIONS  NEEUED 

CHIONE 

42 

c 

CHIONE 

43 

c 

CHIONE 

44 

c 

POPij(N) 

(5EI5.B) 

CHIONE 

45 

c 

CHIONE 

46 

c 

INITIAL  DISTRIBUTION  OF  POPULATION  AMONG  The  VIBRATIONAL 

LEVELS 

CHIONE 

47 

c 

CHIONE 

48 

c 

CHIONE 

4R 

c 

TEMP 

( IE15.8) 

CHIONE 

SO 

c 

CHIONE 

51 

c 

TtMP  =  TEMPERATURE  OF  THE  MODEL  SYSTEM  (KELVIN) 

CHIONF 

52 

c 

CHIONE 

53 

c 

CHIONE 

54 

c 

wS.RP.pN 

(3FIS.B) 

CHIONE 

B5 

c 

CHIONE 

86 

c 

wS  ANO  WP  are  PUMP  FREQUENCIES!  ON  IS  THE  POPULATION  TRANSFER  HTW 

CHIONF 

57 

c 

ADJACENT  VIBRATIONAL  LFVElS.  THE  PROGRAM  ENOS  KMEN  ALL  »«E  ZERO. 

CHIONE 

58 

118 


non 


FTN  4.6.43JB 


1 2/01/78  09.32.2S 


DEBUG 

ARRAYS 

&OTOS 

REAL  LINEWlS.S) 

OlMENS[ON  U£LrAE(S<)0>  .U(5.SI  .LlNfW(S.S)  »POPU(5> 

OlMfNSlON  ChIRE (3.2>  «CMl IM(S«2>  .LABEL (M> 
CUMMON/VlBLOM/WE.WfXE.wEYE.WEZE 
CUMmON/KO  r LOM/HE . AlPhA . GAMMA *  OE . BE  T  A 
CUHMON/TCuM/TfMP 

rear  in  data 

REAO(S.lOl) (LABEL ( 1 1 « 1  =  1 .Bl 
Ini  FORMAT (8A10) 

m2  format uhi »// .sx.  •calculation  of  the  first  order  susceptibility 
=AN0  RE|_  ATtU  PROPER?  IES» , / .SX ,8Al 0»/> 

RtAp<S,10J>NV.NJ.N0VT 

103  FORMAT (3IS) 

WNITEIA. 104INV.NJ.r0VT 

1 04  FORMAT)/. SX,*THIS  CALCULATION  USED*. IS, •  VIBRATIONAL  LEVELS  ANO«. 
SIS,.  ROTATIONAL  levels  with*. IS. •  OVERTONES*) 

REAo<S,10S)we.wtxE,wEre.W£2E 
los  FORMAT (SE1S. 8) 

WR iTt ( 1 06) 

loo  F0RRAT(/.SX,»VIBRATI0NAL  CONSTANTS  I •  .6X  »  *WE*  * I2X,*WEXE*»ll X  »*WCYE* 
s ■  1 1 x.*wE2t*. 1 2x • • < IN  RECIPROCAL  CH)*> 
wRlTEtfc.lOT)  WF.WEXE.WEYE.WEZE 
10?  FURMAT(26x.seiS.a) 

REAP ( S. I  OS ) BE . ALPHA .GAMMA. 0£ .BETA 
wRlTE(ft.lUO) 

1 08  FORMAT </.5X, *roT at lONAl  CONSTANTS *•» 7X»«BE* , 1 2x . *ALPHA* . 10X . 
:»GAmMA»,HX,®0E*.12A.*BETA«.12x»»<IN  RECIPROCAL  CM)*) 

W«ITE<R.1u/)  BE  •  ALPHA,  f,  ARM  A  »0E.  BET  A 
WRITE(R»109) 

I  OR  FORMAT ( 1H1 ,//.5X,*V IBR A  T IOnAL  TRANSITION  MOMENTS  (IN  DEBYE)*./. 
=qA.»N», l l* .»N*0*. 12X.*n»1*. 1 2X , *N*2* , IZX.*N»3*» 12X.*N.4*1 
mOVT^NOVT  *£ 

DO  10  N  =  1 *NV 

REAp (5, 10S) (U(M.N) .M*l .MOVT ) 

NN=N-1 

WRITE (ft. 1 1» )NN. (U(M.N) ,H=l .MOVT) 

1)0  FORMAT (/ .SX , IS.SX «5E IS. 8 ) 

10  CONTINUE 

WRlTE(R.lll)  „  „ 

111  FORMAT (//«3x, ‘TRANSITION  LINEWIOTMS  (In  RECIPROCAL  Cm)*./. 
r9X,*N*,llX.*N»0*. 12X.*N*1». 12X,*N»2* . I 2X , *N* 3* . 1 2X . *N.4* ) 

no  20  n*i*nv 

REAOI5.10S) (LINEW(M.N) ,M»l,MOVT) 

nN*n-1 

WRITE (R. I iZ INN, (LlNtW (M.N) ,M*1 .MnVT) 

11?  FORMAT (/.5a, IS, 5x.SE IS. 81 
20  CONTINUE 

REAP (5.1  os  I  (POPU(N) ,N* 1 ,NV ) 

WKITF(R.IIJ)  .  „ 

in  FORMAT (//.Sx.*THE  INITIAL  VIBRATIONAL  LtVEL  POPULATIONS*,/. 


CHI ONE 

CHI ONE 

CHI ONE 

CHI ONE 

CHI ONE 

CHIONE 

CHI ONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIOnE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 

CHIONE 


I 


onn  non  non 


T 


*  *"-*’*•  ■  -I.  '1  "  "*-v  "V  j  ipa  li  jljHUJfll 


c 

c 

c 


CHI ONE  74/7*  OPT=2  F TN  4.6.433B  12/01/78 

OP. 32.25 

=22X.»0»*14X,»t»,l4x*,2«« 14X»*3*.14X.*4*1 

CHI ONE 

116 

WHITE  16 • 1 14)  (POPUlNl.Nsl.NVl 

CHIONE 

117 

114 

FORMAT (/. 15x.SE 15. 8) 

'HI ONE 

118 

RtAn<5. 10^>»  TEMP 

CHIONE 

1 1  9 

W«lTEl4*US)TEMP 

CHIONE 

120 

115 

FORmaT(//.5x. ‘TEMPERATURE  FOP  THIS  CALCULATION  IS  ‘.EIS.B.*  KELV*> 

CHIONE 

121 

CHIONE 

122 

NORMALIZE  the  population  distribution 

CHIONE 

123 

CHIONE 

124 

poptot.o. 

CHIONE 

125 

no  30  N*1  *NV 

CHIONE 

126 

70 

POP  TOT  atPOPTOT  .POPU  (N ) 

CHIONE 

127 

C)0  40  N*1  *HV 

CHIONE 

128  ' 

40 

POP)  1  ( N )  *P0PU  ( N )  /POP  T  0  T 

CHIONE 

129 

CHIONE 

130 

REAO  IN  WS.WP.ONI  TERMINATE  The  PROGRAM  IF  ThEy  ARE  ALL  ZERO 

CHIONE 

171 

CHIONE 

132 

4S 

REA0<5. 105lwS.WP.0N 

CHIONE 

133 

lFtwS.EO.0..Af|O.wP.EQ.0..AND.ON.EQ.0.l  GO  TO  1000 

CHIONE 

134 

WRITE  <6« 1 16)  WS.WP.ON 

CHIONE 

135 

114 

FORmaT(////.5X.»THE  PUMP  FREQUENCIES  ARE» , 2 (5X ,E 15 .8 > , / . 5X ,  ‘THE 

CHIONE 

136 

rPOPuLATlON  tRahSFEmEU  BETWEEN  ADJACENT  LEVELS  I S . »5X . E 15. B ) 

CHIONE 

137 

CHIONE 

138 

begin  the  Calculation  of  ch!  it  wS  and  wp»  Finn  acoef  and  pemion 

CHIONE 

139 

CHIONE 

140 

00  60  K  *  1 . 2 

CHIONE 

1  4  l 

w=ws 

CHIONE 

142 

lFlw.E0.21  W=WP 

CHIONE 

143 

CHlRETsO. 

CHIONE 

144 

CHI [Ml eO • 

CHIONE 

145 

00  50  IGROON=1.NV 

CHIONE 

146 

call  energy i igroun.nj.novt ,oelt ae .0) 

CHIONE 

147 

CALL  CHI < 1 GROUN. NJ. NOVT. W.O.U.LINEW.DElTAE. CHIRE (IGROUN.K) , 

CHIONE 

1*8 

:CHI IM( iGWOUN.K) ) 

CHIONE 

1*9 

CHlRETaCHlRET’POPUl IGHOUNI*CHIRE( IGHOUN.K) 

CHIONE 

ISO 

so 

CHI  [MT*CHl  IMT.POPUl  IGROUNI  »CH1  IM  (  If.ROUN.K  ) 

CHIONE 

IM 

Z=SoRT ( I ♦4.*3. 141 5927#CM I  RE  T  1 

CHIONE 

152 

ACOEF >4.* J. 141 5927» W*CHI 1MT/Z 

CHIONE 

153 

rEMi0NzZ-1. 

CHIONE 

IS* 

W«lTE  16. 1 17)W,AC0EF*REMI0N 

CHIONE 

155 

117 

format (/,5Xf»F0R  The  FREQUENCY  *.E15.8,/»5X,*THE  ABSORPTION  COEFIC 

CHIONE 

156 

sILNT  IS  ».tlS.B./.5X.»THE  REFRACTIVE  INDEX  MINUS  UNITY  IS  ••EIS.S) 

CHIONE 

157 

41* 

CONTINUE 

CHIONE 

158 

lFlNW.EO.il  GO  TO  45 

CHIONE 

159 

CHIONE 

160 

begin  the  calculation  of  the  k-vector  mismatch 

CHIONF 

161 

CHIONE 

162 

lTOp*NV”l 

CHIONE 

163 

po  70  I*1«I TOP 

CHIONE 

164 

[LOWER* l 

CHIONE 

165 

I  UPPER* I ♦ l 

CHIONE 

166 

wAVEM|*WS»ON* (CHIRE < 1UPPER. 1 1 -CHIRE 1 ILoWER. 1 > -CHI  RE < I UPPER. 2) ♦ 

CHIONE 

167 

rCHlREI IL0WLR.2) ) »2.*J. 1 415927 

CHIONE 

168 

IL* ILOwER- 1 

CHIONE 

169 

IU*  j upper-i 

CHIONE 

170 

WRITE (6» 1 10) IL. IU.WAVEMI 

CHIONE 

171 

1  1« 

FORMAT  (//.5X.  <*THE  N-VECTOR  MISMATCH  VARIATION  DUE  TO  POPULATION 

CHIONE 

172 

v 


120 


CM t ONE 


74//*»  OPT*2 


FTN  4.6.433b 


iTMAnSfE*  BETWEEN  LEVELS*.  IS.*  AND*.|S«*  IS  *.E15.B> 
70  CONTINUE 
GO  TO  45 

1000  WHlT£<6»120) 

1  pi)  FUNm*T<////,5x.*THaT*S  ALL  FOLKS*) 

STOP 

cun 


lP/01/78  09.32.25 

CHI ONE  173 
ChlONE  174 
CHI ONE  1 75 
CHI ONE  176 
CHIONE  177 
CHI ONE  178 
CHIONE  179 


nonnn 


ENERGY 


74/74  OPT  >3 


FTN  4.6.433B 


12/01/78  09.32.2S 


SO&ROU T 1 NE  ENERGY ( I GROUN ,NJ .NOV T . OEL  T  A£  *  0  > 

THIS  SUBROUTINE  CALCULATES  The  NEEDED  ENERGY  LFVELS  AND  PARTITION 

function*  Energy  level  differences  are  stored  linearly  in 
oEltae.  the  PARTITION  function  is  stored  in  0 

C*  pEbuo 

CS  ARRAYS 

C*  GOlOS 

DIMENSION  UELTAE(5(*0) 

COMMON/ V I HCOM/WE . RE  XE . RE  YE , RE  2F 

common/Ro r  com/be . al  pha . gamma . of . be  t  a 

COMMON/ TCOM/ TEMP 

ibranc.-i 

ILOreRxIGHOun-I 

IUPpeRx1GHOuN-1 

ISToRExO 

hE=2.»0E»< 12.»flE»«2-ALPHA»wE)/(3.»WE,#2) 

8  IROtxO 

X  =  I|.0MFR»0.5 

B=bE-ALPHA»x.r,AMMA»X»»2 

D=Df.BETA»X 

1 1  Ya««IRoT» ( I ROT . 1 ) -0#  < I ROT* ( I ROT ♦  1 1 ) **2.HE* ( IROT» < IROT.l ) )**3 
Z=MF»X_WEXE»X»*2.«(FYE,X»*3.MEZE#X*»4 
INOeX  =  JROT • I STORE  * l 
OELTAE ( INDEX) x-Y-Z 
|ROT  x IbOT ♦ 1 

lE( IROT.LE.NJ)  GO  TO  II 
1«OT=0 

x*iuppfR*«.5 

B=bF-ALPHA«x.GAMMA»X»»2 
d=0f«8fta*x 
18  JPOtx IrOT ♦ iRPANC 

Y=B«JRoT»( JP0T«n-0*(JH0T»( JROT‘1 > I ••2»HE* (JROT» I JROT . 1 )  )<*»3 

Zxmf*x-wexe*x**2.mevE*x**3.MEZE*X**4 

inOfx* iRO  r ♦ I STORE* l 

OELTAE ( INDEX )=OELTaE( INDEX) »Y .7 

IROt=IrOT.1 

lE <  IROT.Lt.NJ)  GO  TO  18 
IKOTxO 

24  jHOT*  tpOT ♦ IRRANC 
IN0EX*IR0T ♦ I STORE* I 
lE( JP0T.LT. 0)  OELTAE* INDEX >*0. 

|WOTx IpOT ♦ 1 

If* IR0t.LE.NJ)  go  TO  24 
I E  1 1 BPANC *EO . • 1 )  GO  TO  27 
ITEsT*IGROUn*NOVT-1 
lE( lUPPEH.LE. I  TEST )  GO  TO  28 
B*«F.-AlPHA/2.  .GAMMa/4. 

0*6.9S03E-01*TEMP/B 

RETURN 

27  I BK  ANCx 1 

lST0PExNJ*IST0RE 
GO  TO  8 

28  IHHaNCx-1 
IUPpeRxIUPPER‘1 
ISToRExNJ* ISTOPE 
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ENERGY 

3 

ENERGY 

4 

ENERGY 

s 

ENERGY 
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ENERGY 

7 

ENERGY 
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ENERGY 
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ENERGY 

10 

ENERGY 

11 

ENERGY 

12 

ENERGY 

13 

ENERGY 

14 

ENERGY 

15 

ENERGY 

16 

energy 

17 

ENERGY 

18 

ENERGY 

19 

ENERGY 

20 

ENERGY 

21 

ENERGY 

22 

ENERGY 

?3 

ENERGY 

24 

ENERGY 

25 

ENERGY 

26 

ENERGY 

27 

ENERGY 

28 

ENERGY 

29 

energy 

30 

ENERGY 

31 

ENERGY 

32 

ENERGY 

33 

ENERGY 

34 

energy 

15 

ENERGY 

36 

ENERGY 

37 

energy 

38 

ENERGY 

39 

ENERGY 

40 

ENERGY 

41 

ENERGY 

42 

ENERGY 

43 

ENERGY 

44 

ENERGY 

45 

ENERGY 

46 

ENERGY 

47 

ENERGY 

48 

ENERGY 

49 

ENERGY 

SO 

ENERGY 

SI 

ENERGY 

52 

ENERGY 

53 

ENERGY 

54 

ENERGY 

SS 

ENERGY 

56 

ENERGY 

57 

ENERGY 

58 
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L 


4 


noon 


74/74  OPT*? 


FTN  4.6.433b 


12/01/78  09.32.2S 


CHf 


SUBROUTINE  CM ( (N.NJ.NOVT .W.Q.U.L  INEW.DELTAE.CHIfiE  »OlI  IH ) 

CHI 

2 

c 

chi 

3 

c 

THIS  SUBROUTINE  CALCULATES  THE  REAL  AND  IMAGINARY  PARTS  OF  CHI 

FOR  CHI 

4 

c 

PARTICULAR  LOWER  state  vibrational  LEVEL 

chi 

s 

c 

Chi 

6 

c* 

OttJilG 

Chi 

7 

cs 

arrays 

chi 

8 

c* 

GUI  OS 

Chi 

9 

Rt  Al  !_  I  N£  *  1 5  «  5 ! 

CHI 

10 

0 1 MENS I ON  OELTAE (SO0> .U(5.5) .LINE* (5. 5) 

Chi 

11 

COMMON/ V 1 UCOM/wE . WE  At . WEYE . WEZE 

Chi 

12 

common /ro i com/ue • alpha  <  gamma • ok . be  t  a 

Chi 

13 

CUMmON/TCUM/TEMP 

CHI 

14 

irot*o 

CHI 

15 

I51 ORE  *0 

Chi 

16 

jSlORE=0 

CHI 

17 

IbRANCs-1 

Chi 

18 

SUMA*0. 

Chi 

19 

SUMB*0. 

chi 

20 

SUMC=0. 

Chi 

21 

SUMQsO, 

Chi 

22 

CH1R£*0. 

Chi 

23 

CHIIM*o. 

Chi 

24 

6  INOEX*IROT.ISTORE*1 

Chi 

25 

lE  < OELTAE  * 1NOEX ) . £Q . 0 1  GO  TO  16 

Chi 

26 

jnOeX*jSTURE»1 

Chi 

27 

Z*<OELT*E< INDEX)  .W)  ••2»0.2S,LIN£V((  JNUEx.Nl  ’••2 

CHI 

28 

SUM  13 (OELTAE I  INDEX) «WI/Z 

CHI 

29 

Z=‘OELTAEllNOEXI-W)»*2.0.2S»LINEW(JNOFx.N)»«? 

CHI 

30 

SUM?*<oELTAtlINOEX)-R)/Z 

CHI 

31 

SUMr,xSuMl  «SUM? 

CHI 

32 

Z= 'OELTAE  UNOEX 1  .WI«*2.0.2S»LINEwIJNOEX»N|»»2 

CHI 

33 

SUMi=-o.5«LINEW(JNOEX,N)/Z 

CHI 

34 

z=<oeltaei index >-w) **2.0.25^ inewijnoex.n>»«2 

CHI 

35 

SUM?* ♦ o .5*L INEW 1 JNOEX . N ) / Z 

CHI 

36 

SUMH=5uM1.SuM2 

CHI 

37 

IS  Z* lRO T ♦  1 

CHI 

38 

1EIibHaNC.EQ.-1)  Z=1H0T 

Chi 

39 

SUMA*Z*SUMG.SUMA 

CHI 

40 

SUMb=Z*SUMh.sumB 

CHI 

41 

lE ( IBRaNC .EG.- 1 )  GO  TO  ia 

CHI 

42 

Z*U(JNOEX.N)»*2 

CHI 

43 

SUMCsZ*SUMA.Siimc 

Chi 

44 

SUMO*Z»SUMB.SOM0 

CHI 

45 

IE1JST0RE.LE.N0VT)  GO  TO  20 

CHI 

46 

Z=IROT*<IMUT.l ) 

CHI 

47 

hE*2 . •()£•  1  12.»hE»#2*ALPHA*wEI/  (3.*wE##2> 

CHI 

48 

Z*<RE-ALPHA»0.5.GAMMA«0.25) •Z-(OE.BETA»0.5)»Z»*2.HE«Z»«3 

CHI 

49 

2*(?»IpOT.l)»EXPl-Z/B.V503F-0l/TEMR)/U 

CHI 

so 

CHlREJ,SUMC»Z.CHlRE 

CHI 

51 

ch!im*sumu*z.chiih 

CHI 

52 

IF  1  lR0T.Lt.NJ)  GO  TO  27 

CHI 

S3 

CHlREM .3*7E-1»ChIrE 

CHI 

54 

CHI IM*1 .35  7£- 1 *CHl I« 

CHI 

5S 

RE  I  l)RN 

CHI 

56 

16  SUMfisO. 

CHI 

57 

sumh«o. 

CHI 

58 
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Sample  Run 


c  ac  cut  Arrow  or  rwc  rm*r  o*t>€*  susccrTr*iurr  anO  aclatco  *wrc*fics 
oeuiciuuw  rvuoNioc 


THIS  CALCULATION  USFO  t  V  ISA AT  I ON AC  wCVCLS  ANO  30  *0 TAT I ON AC  LEVELS  WITH  1  OV CO TONES 


VIONATIONM.  CONSTANTS!  »€ 


WE  AC 

.*Ti*oo#«e»u 


wCvC  WCZC 

V3*0«0000C»06  .30 30*00 OE-02 


ROTATIONAL  CONSTANTS! 


•C  acpna 

•  iiooirooe«o2  .j«mhmc»oo 


GANNA  OC 

. J500000*C-02  .SOCOOOOOC-*3 


(IN  RCCtrNOCAL  CHI 
OCTA  (IN  OCCIRNOCaL  CM) 
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N«* 


tiMiriowu.  transition  "oncnts  (in  oea*e> 

M  N««  N.l  N.i  N.J 

«  .1RI44M*C**I  .»*10*0»#C-«1  -.MIMMIWJ 

TRANSITION  LINCRIOThS  >lN  4CC1RN0CAL  CNl 

N  N<«  N.l  N.2  N.J  N«4 

«  ,«i}IIUN’l«  .40*40080*. 08 

THC  INITIAL  V IRRATIONAL  LEVEL  •'ONULaTIONS 

•  1  2  J  a 

•  lOOOMOOC’Ol 

TENRCRATURC  TOR  THIS  CALCULATION  IS  . JROROROOC’OJ  KCLv 


THC  RUHR  FNCOUCNCICS  ARC  .AOQOOMRC’A*  .1S40000*C*04 

THC  RORULAT ION  TRANSFCRCO  RCTRCCn  tOJACCNT  lCVCLS  IS.  . 1880I0.0E-84 

RON  THC  FRCOUCNC Y  .*000000 «••* 

THC  ARSORRTION  COCriCIrNT  IS  .5»10*0»*C-04 

THC  RETRACT l VE  INOFI  HluUS  UNITY  IS  -.3RAA44RIC*** 

ROR  THC  FRCOUCNCY  .  ]S*SOOOOC«0* 

THC  ARSORRTION  COCFIClCNT  IS  . 34SA47ME-03 

THC  RETRACTIVE  INOfi  hImUS  UNITY  IS  -.JA21R2TAC*0J 


THC  RUHR  ERCOUCNCICS  aHC  .4OO00001E "0*  . 154S0000C-0* 

THC  PORULATION  TRAnsFErEO  ReTVCCn  ADJACENT  LCVCLS  IS.  .lOOOOOOOE-O* 

TOR  THC  TRCOUCNCT  .*0800000E*0» 

THC  ARSORRTION  COCRIClENT  IS  .5*lO*4«A£-«* 

THC  RETRACTIVE  1NOEI  nInUS  UNITY  Is  -.TRAA4AA1E-0* 


FOR  THC  TRCOUCNCT  . 1S«SRR0RC*«* 

THC  ARSORRTION  COCriClENT  Is  .3A373412C-03 

THC  RETRACT  IVC  INOER  Hlhus  UNITY  is  -.JSMOOMC-RJ 


THC  RUHR  TMCOUCNCIER  *HE  .60000000E.04  . 15500000E-0* 

THC  RORULATIOM  TRANSFCNEO  ReTNCCN  adjacent  LCVCLS  IS.  . IOOOOOOOC-04 

TOR  THC  TRCOUCNCT  ,*OOoOOOOC*04 

THC  ARSORRTION  COCriClENT  IS  .SRlOSORRC-O* 

THC  RCTRACTIvC  INOFI  nInoS  UNITY  Is  -.3RAA4RRIC-0* 

FOR  THC  TRCOUCNCT  . lSSo0000C*04 

THC  ARSORRTION  COCriClENT  Is  .JA1A4AORC-OJ 

-  -  . . .  .Hire  » i-  . 
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thc  pun*  raeouoreies  me  .»ooooaooc-*»  .isssgoooc.o* 

THC  POPULATION  TPANSPCMEo  U'XM  .OJACCNT  LEVELS  IS.  -Hal 

ro*  thc  r«eau*NCY  .gO0ao*»*c.«» 

thc  .0 so** Mom  cocrtcjCM  is  .saiaagiaE-** 

THC  0Cr*ACTIV<  INOP*  MlnUS  unity  is  •.imimii-n 

re*  thc  raeouENCY  .ilSjuuft. 

THC  .*SO*PTION  COCf ICItHT  Is  . 359S«3J&£-«3 

thc  »cr»Acrivc  inert  nines  unity  is 


thc  pu«p  racoucNciE*  me  -goaoaoooc.on  .issosoooc-o* 

THC  POPULATION  rSAMSrextO  scthCEn  aouacent  LEVELS  IS.  .1001 

ro*  thc  r*COue*cv  .»ao»oa»ac.o* 

THC  »*SO«PTION  tOCYICItHT  ts  .SaiOSOSOC-O* 

THC  KCCNACTIVC  INOP*  MINUS  UNITY  IS  -.]**«*«*  IE*** 

ran  tm€  r*C0UCNCY  .  iS*oooo«C*0a 

THC  A*S0*PTI0N  COCriCjCNY  IS  . 3STS»53«t-03 

THC  0CYNACTIVC  INOP*  MINUS  UNITY  Is  -. JSI 1 9820C-«3 


thc  pum*  FNeauCNCirs  me  .tootogooco*  .is*sooooe*g* 

thc  population  r»*NsrE*eo  sctnccn  .ojaCcnt  levels  is.  .looooeooc-ofc 

TOP  THC  EPCOUCNCY  . 48«0«00C£ •«* 

thc  absonpt ion  cocricieNY  is  .s*io**«ae-oA 

THC  PCPNACTIvC  IHOC>  MINUS  UNITY  IS  ••]M***S1E*«* 

ro*  thc  rpcauCNCY  .iS«s«agoe.o* 

thc  msoMprioN  coenciENr  is  .jsss32i9C-*3 

THC  PCPPACTIvC  INOP*  mImUS  unity  is  -.3»*S1326C-*3 


thc  pu«p  r»eaucNcips  i»£  .g«o#oggee-o»  .istoooooc.o* 

THC  POPULATION  TRaHSPCXCn  aclMCCN  .OJACCNT  LEVELS  IS.  .IOOO«goOE-»» 

ro*  THC  r*eouCNCY  .»o«oogg«£.o» 

thc  *bso*ption  cocftciCny  is  .s»ioaov*c-o» 

THC  PcrPACTIVE  INOP*  Minus  UNITY  is  -.3*M»aeiE~SN 

ro*  THC  rneouCNCY  .iSTgoOgOfO* 

thc  .asoMP'ioN  cocrinENT  is  .3S3S«3**e-a3 

THC  PCPPACTIVC  INOP*  minus  unity  is  3»5aSJlZE-03 


thc  pump  ppcoucncifs  **e  .gaoaaoaoc.o*  .i5Tsooaoe*o* 

the  population  t*anspchco  Sc tmeen  .ojaccmt  lcvclS  is.  . ioooo»aoe-o* 

ro*  thc  r*couCNcY  .gaageooac.o* 

thc  lasoaprioN  cocrictCNY  is  .saigaaVMC-a* 

THC  PCPPACTIVC  I  HOP*  Minus  UNITY  Is  -.jS6»i.*OIE-0* 


I 


re*  rut  frequency  . isJsaoaoc.o* 

THE  ABSORPTION  COEETCHMT  is  .MISIMMC**) 

tmc  refractive  inoei  kIwjs  uni t»  is 


tmc  pu«p  frequencies  are  .aoaaoaaac.a*  .isiaaaaoc.a* 

TMC  POPULATION  TRaa«5EC«10  SETveEN  aOJaCE"!  LEVELS  IS.  . 1 OaSOSaaC-a» 

re*  TMC  FREQUENCY  .»v«0«OO»E-»* 

tmc  AH soap I  ion  coc/tcieht  is  .saioao9ac-o* 

TMC  REFRACTIVE  INDEX  Ml >, us  uniiv  is 

EON  THE  FREQUENCY  . 

TMC  ABSORPTION  COCEICIEmT  IS  .3*9*3BB7E-a3 

THE  REFRACTIVE  INOCX  nInuS  unity  is  -.J»aS9991E*93 


TMC  PUMP  ENCOUCMCIES  a«E  .60000»»«E-«a  . 15SS009a£.0* 

TMC  POPULATION  TBANSEEPEO  8Ct»EEn  aOJACENI  LEVELS  IS. 

EON  TMC  FREQUENCY  .09000*0 OE.O* 

TMC  ABSORPTION  COEEICIENT  IS  .Stl0*0*»E-«* 

THE  RETRACTIVE  INDEX  nImuS  UNITT  [S  3****401£-0* 

FOR  THE  FREOUENCT  . lS»S000fE*«» 

THE  ADSORPTION  COEEICIENT  IS  .3A772130E-03 

THE  REERACTIVC  INDEX  nImiS  uNIIt  IS  -.33M0S11£*<3 


THE  PUMP  ENEOUCNClCS  aME  . AO«08000£*0 A  . 1590000oE*8a 

the  population  transeereo  setveen  adjacent  levels  is.  . laoaoooac-OA 


FOR  TMC  FREOUENCT  .6«000»00E-»» 

THE  ABSORPTION  COEEICIENT  Is  .SBl0«O9OE-a* 

THE  REFRACTIVE  INOEX  nInuS  uNITt  IS  -.3BAOAOSIE-4* 

FOR  Tm£  FREOUENCT  . 1PV«000«C*VA 

TME  ABSORPTION  COCFICItNT  IS  .3*S02799C-a3 

TMC  REFRACTIVE  INOEX  MINUS  UNITY  Is  -.33SA3B9*C*03 


the  pump  frequencies  are  .Aaaoaaaoc>a«  .iS’soooaE.o* 

THE  POPULATION  TRANSEEREO  BeTHLEn  ADJACENT  LEVELS  IS*  .lOOOOOOOC-O* 

eor  the  frcoucnct  .avaoooaac.a* 

tmc  absorption  coceicient  is  .saioao»aE-«* 

the  refractive  inoei  hInus  unitt  is  -.jBaaAsaic-a* 

For  the  frcoucnct  .isssaaaaca* 

THE  ABSORPTION  COEEICIENT  IS  .3AJ9S7R7E-03 

THE  REFRACTIVE  INOEX  minus  UNITY  Is  -.332I9192C-03 


the  pump  frequencies  ahe  .aaoaooaoE.o*  . iaooooooE.o* 

.  _  --»•••  ......  . «  ..««  i.»m.  i«.  .  1  onnnonOF-o*. 
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rod  tmc  racoucNCY  .a»soo8s«.vn 
rut  AdSoddTION  coericiCNT  is 

TMC  WdACTIVE  I  HOC*  MINUS  UNITY  is  -.JMIMICM 

roa  tmc  rMOucNcr 

TMC  UMWMW  COCFTCltXT  IS  .3N*lltdtE-S3 
me  x*r»*cTivc  me ex  minus  unity  is  •  .jj*)««?m*i3 


tmc  mum*  rocaucNcics  <«  .utiiMH-u 

TMC  MOdlJLAf  ION  TMAMSMCNCo  SCTNEEN  >OJtCCMl  LEVELS  IS.  .  lOSAOMK-SN 

TON  TMC  rdEOUCNCT  .Att'MStC.tN 

TMC  »OSO*OT|ON  OOCYiejENT  Is  .St1t**«SC>*N 

TMC  •CYNScrive  tMOri  MINUS  UNI  IT  IS  • .3«d«NtSlC-IN 

rod  tmc  ykcoucncy  .imssscsc.cn 
tmc  stseddfioN  cocrtcitNT  is 

TMC  dCTdACTUC  I  HOC  I  minus  UNIIT  Is  -.3EY|t3|TC*C3 


tmc  sunn  FNcauCHcirs  »(  .cssstssst'S*  .isiooosoe.o* 

tmc  momvil * i  ion  tramstcpCo  sc'mccn  .ojaCCMT  level*  IS.  - 1001 

rod  TMC  rdEOUCNCT  .SSSOOSSC'S* 

tmc  Aisoddf ion  cocftcicnt  is  .sdiodsssc-s* 

TMC  MCYMACTIVC  INOCT  minus  unity  is 

rod  TMC  YdCOUCNCY  • isioosssc.cn 

tmc  assonmtion  coericiCNT  is  .jssnsnstc-ss 

tmc  dCMdACTIVC  I  MOT*  minus  UNITY  Is  J233S*«SC-*C3 


tmc  mumo  mcoucMcirs  *nc  .«sosssc«c. sn  . i»isoso«.«* 

tmc  momulaTion  TdAMsrcMCo  ac'MCCM  adjacent  levels  is.  .iossvossc-sn 

rod  TMC  rdCOuCNCY  .  6S04000SC . SN 

TMC  AOSOddTION  cocr tciemt  is  .muhhi-1. 

TMC  »Cr»ACfI»C  INOFI  MINUS  U»|TY  Is  -OMNNdSIC'SN 

rod  tmc  rdeaucMCY  .i*isssosc.«n 

tmc  AdSoddTION  coericiCNT  IS  O3A7»Sd0E-»3 

TMC  BCMdACTive  !MOM»  nUus  unity  is  -oszsissic-ss 


tmc  mun*  mcoucMcirs  *«e  .nsoooosoC.sm  .iaJaoOoc- 

tmc  MOSUL* T ION  TdAN5rC*(0  SCTnEEn  adjacent  levels  is.  I 

rod  tmc  rdCouCHCv  .AtCDosooc.ts 

TMC  AdSoddTION  coericiCNT  IS  .SitOSSSdC-Ss 

TMC  "CMNACTIYC  I HOT I  MINUS  UNITS  is  •.JtdNNdSIC’** 

rod  tmc  rdCouCMcr  . i»2o»«»sC‘8» 

TMC  AdSoddTION  coericiCNT  is  .33N»Ndd»E-«J 

tmc  McrdACTive  iNocv  minus  unity  is  -.j2i»T«t7e-tj 


no 


«  |||  8**  ||  8**  8**  ||  8**  |||  ||  SI*  HI  ||  HI 


w»  frequencies  »<*t  .mhhim>m 

POPULATION  TRANfFC«tO  HlnUn  aOJaCENT  lCVClS  If,  .IttHHK-W 

THt  FREOuENCS  .*OOOOOOOE»»* 
absorption  coeficunt  is  .siuitwt-at 

RCFRACTIVE  I  HOT  l  Mlwuf  unjtc  is  •  ,]»44*00IC*04 
THt  FRCOUCNCS  .IMWIIKhU 

absorption  cocficiCnt  if  o.mi]97S-03 

REFRACTIVE  INOCI  HlNOS  UNITS  If  -.HUfltU-t] 


FREQUENCIES  A*t  .MtllUK-ft  .IIMIMK'!* 

POPULAT ION  TRAN5FCHtO  M'hCAN  AOJACENT  LlVCLl  |f,  .1001 

ihc  frcouEncy  .oootooooc><)« 
absorption  cocficiCnt  if  -smuMi'H 

»CFRa CTIVI  I  NOT  I  Hi  HUS  UNITS  If  -.MWtMIC-M 

THt  frcquCncs  . i*Soooooc,«h 
ABSORPTION  COCFIClCNT  If  OJIS0003C-03 
REFRACT IvC  INOCI  Mluus  units  if  •.31Mf040C*0) 


PUMP  FHCOUCNCItf  ARC  .400000001,00  . |4JSOOOO(«44 

POPULATION  T»»wt«lO  OCTUCCn  aOUACCNI  LtVCLS  If.  .1001 

THt  FRCOuCNCT  .4<OoaOOOC'0» 

ABSORPTION  COCFIClCNT  If  .MIOOO*OC>44 
REFRACTIVE  [NOCI  HlwJS  uN|Tr  |j  J«*4*0#lt-0» 

THt  FRCOUCNCS  .I«lf0000C>0* 

ABSORPTION  COCFIClCNT  If  OMOOTM-I] 

HCFRACTIVC  INOCI  H«nuS  unit*  IS  -.1U2TU0C-O] 


PUMP  FHCOUCMCIFf  AMC  .00000000C.0O  ,I4*00S00C,04 

POPUL»;'JN  THPNfFCNto  OCTuCCn  AOuACCNT  lESClS  IS,  .IOOOOOOOC-04 

THE  rHCOUCNCT  .oOOoOOoOC'O* 

AOSuMPT iOH  COCFIClCNT  If  .SOIOOOOOC-44 
HtFNACT  I  VC  |N0f«  MINUS  (P*l  Tf  JJ  -.MOOWllt^i 

THt  FRCOUCNCS  .l«NOOOOOC>P* 

ABSORPTION  COCFIClCNT  If  .»012«4«c-03 
HCFRACTIVC  I  NOT  I  MINUS  UNITS  If  -.JI09I  TWC-Of 


pump  FRcaucNcirs  *“t  .ooooooooc,o»  .iosooooot«i 

POPULATION  TRAnSFCMCo  ScTSCCn  ADJACENT  LCVClS  If,  . li 

THC  FRCOUCNCT  .k<IOoaOoOC,H4 
ABSORPTION  COCFIClCNT  If  .S01000*hC-04 
RCFRaCTIVC  INOCI  MUuf  OMlTr  If  -.J044A00IC-04 
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rxt  mowi  ion  cocrtcitsf  is  .iwiMncm 

rut  mctsactivc  inoci  minus  unity  is  i»33I7SOC-03 


rxt  NDW  TOSUCNCICS  lM  .UII«ltlC-l« 

tnc  m*w.»noN  T**Msri«to  H'Xbi  tojiCCNT  lcvcl*  is.  .1001001 

res  tmc  rsceucNCT  .«»hhk-» 

TNC  iOJOSST ION  COCf  ICICNt  IS  .lllll|1|C-l* 

TNC  MCTSSCTIVC  INOS0  nInUS  UNITY  is  -.JA***A01C-O» 

ros  TNC  rstauCNcr 

tnc  tasoNsrioN  cecrictCNT  ts  .**i»s**sc-«j 
tnc  scr««cTivc  iNDri  mI«us  unity  is  •.isitmtsc*o 


tnc  sunn  rsewcNctes  mi  .miithu-m  ,i»hiii£*m 

TNC  NONULSTION  TSSNSrCMtO  SCrwClN  iOJACCNT  LCYCLS  IS.  .3001 

ros  tnc  rSCOUCNCY  .|ll|MHC-l. 

TNC  ASSOSST ION  COCriClCNT  IS  ijlllllTtf-N 

TNC  scrsscrivc  INOfi  Minus  UNITY  |S  - ■ JOOiiOO IC-«» 

ros  TNC  rSCOUCNCY  .  1 »»tf  • «» 

TNC  iSSOSSTION  COCriClCNT  IS  ■MIUTllCTI 

TNC  RcrSACTIVC  I  Nor*  MlwtlS  UN|Tt  is  -.1*S2U17C~«3 


tnc  NUNN  r»COUCNClrS  Mt  (MSSSSStC.S*  ,]f*S«*0«C •»* 

TNC  NONW.STION  TNMsrCStO  SCTNCCn  iOJACCNT  LCYCLS  ISi  .1000001 

ros  tnc  rscaucNCT  .ootttsssc-s* 

TNC  iSSOSSTION  COCriClCNT  IS  .HllHTH-M 

tmc  scrsscT  IVC  I  NOT.  MINUS  unity  is  -OSSS*SSIC*SA 

ros  TMC  rSCOUCNCY  •|MSSSSS#C*S« 

TNC  'SU^SSTION  COCriClCNT  IS  .CS10IS3M-S1 

THt  J-.ru»CT|»£  INOri  MINUS  unity  is  -■  |IS*tH3C-«3 


TNC  »«NN  rSCOUCNCIF.  AMC  .MMlHlC-li  . | ITIIIIK •(. 

TMC  NONUIATION  TSUisrCSCO  SCTNCCN  iOJACCNT  LCYCLS  IS. 

ros  TNC  rSCOUCNCY  .aOSoosssc-m* 

TMC  ASSOMNT ION  COCriClCNT  is  .SSI0S*«*C-0» 

TMC  SCrSACTlYC  I NOT I  Mlius  UNITT  IS  •.JSSOASSIC-SA 

ros  TMC  rSCOUCNCY  . lY'OOOOOC-Oi 

INC  iSSOSSTION  COCriClCNT  ts  .«OAO«t  7C-03 

TNC  SerSACTlvC  INOri  MINUS  unity  is  -.isti**?k-«3 


tnc  suns  rscouCNcirs  anc  .aoooooooc.oa  .isfsooogc.o* 

TMC  SOSULATION  TSMsrCSCn  OCl»tCN  AOJACCNT  LCYCLS  is.  .1001 

ros  tmc  rscoucNcY  .*«so«s*sc.sa 
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,  - 5*Kd»[Wi !**• -  •**>  - *..:--■*** »**#* 


Tfcrf  ninirMOH  cocnctCftr  is  .wioiiwi-^ 

acroACTive  imom  »i«us  unity  is 

ro«  thc  fdtoucuer  .i»»so#ooe*«* 

rnt  a*0ddM«*  cwrti»>  tt  .  ^ 

THC  HCfdACTITt  1M0T1  HlnuS  UH|IT  IS 


ae^rsiiuro-s 

rrs 

THt  dCMACTIVC  I NOT I  hImuS  u»i>»  is  -.3»***s»ic  »» 

rod  ThC  rdCOUCNCT  . 

Tktf  aASOR#MQN  COCfTCI^NT  15  *2$S§ff4iC»«J 

thc  RcroAcrivc  (hotx  minus  uhitt  is  -,i»*«**mc*m 


tm*  Mil#  fMCOUCNCtCS  kNl  •MIHQMt*** 

THC  MO#ili.A r ION  THAMsrCBto  *cTyC£n  40JACCNT  LCVCtS  IS*  .UOOOOOOC 


ro#  THC  FBCQUCNCY  ,4«QOO®OOC**^ 

ThC  AAS6Mf»riOM  COCriClCNT  is  »5I|0M9*C*N 
thc  ncrddCTtvC  INOfi  hImts  uH|T»  is  -.J»***MIC-«» 

rod  THC  rdcaocdcr  _ 

THC  »dSO«dfiOH  cocricitHT  is  •»SS‘*‘»St-«3 
THC  dCOMCTlve  IMOri  HIHOS  uhitt  is  -,imshm*c 


thc  rOdU.«lfoNCTd*N4rCl|to  OCThCCH  tOJACCHT  LCVU.S  is, 

t2c  lasoHri^OH^cocricitHT*!8  *  .somosoc-o* 

thc  denucTITC  l HOT,  hImjS  uHITt  is  -.3S**»«S1£  «* 

rod  thc  rdcaoCMc r  .i»scoc#ee*<>* 

me  AosoddT ion  coericitHi  ;s  .*9si»**st-«J 

thc  dcrd»cTi»e  iMOe»  hihos  uhitt  is  ».i»i«J*hJC  ** 


,100 


-  n-.-a  rdrooCNCirs  *hc  •  1  iTTflOTfll* o 4 

THC  POdULAflON  TdAHSrtdCo  OcThICh  »OJ»C£Nf  LCVCLS  IS>  .IQOOOOOOC- 


tS  I5oI!MS1«nit£H?#l?  •.SOIOMSdC-O*,. 

thc  dcrdtcTivc  iHor«  hihus  uhitt  is  ..!•**•« ic  «* 

rod  thc  rdCguCHCT  .isssOOoOC-4* 

THC  »0S£>Hrr10H  cocricitHT  is 

thc  dcroACiive  ihoc«  hIhus  uhitt  is  -.it»sc*37c  « 
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rx  *utn*  rxuutNClis  »~i  •nunvvvvun'vm  .MVM«M*Mm.’wM 

rx  *o*ul*tion  toansteXo  K'XE*  »oj»cent  level*  is,  .laai 

iM  rx  xeouexy  .«**••***<.*« 

rx  umwiion  cocrieiCMr  is  .saioaa»*E-*» 

TX  *CT»»CTIvE  INOTI  hImUS  unity  IS  >.1MMMIM» 

ro »  rx  xcautMcr  .iicitiicx 

rx  xsoxhom  cocticient  is  -MsamsE-aj 

rx  R£T*»CTlvS  INOTI  minus  unity  is  •.iTSStMlCM] 


rx  nun*  nniueiii  xt  .mmhnix  .iiiioimC'H 
rx  *o*w.ation  nsNsreXs  oetxen  *ojacent  uras  is,  .loai 

ro*  rx  rosouEXY 

rx  xsoxrioH  coErtciCNr  is  .miimX-i* 
rx  *cr**crivc  inch  minus  unity  is  -.j**o*aaiE-«» 

ro*  rx  rxoxNcr  .iMniiK'it 

rx  ***OM*  MON  COETICtCNT  IS  .«OJ*l**C-*J 
rx  «ET»»CTIve  inoti  minus  unity  is  -.itosesmehs 


rx  MU**  r*C0UCNCIts  UK  .UtllMlC-l, 

rx  motulmiion  toansteoeo  ktwin  iojacint  levels  is,  .1 

ro*  rx  enequency  .»saao***C>a« 
rx  iosortmon  cocrieiENr  is  •mihiW-h 
rx  »ET*acTivE  inoti  minus  unity  is  -.I***«a*tE-«* 

ro*  rx  rxouexr  ,joioo»o*c-»* 

rx  ASSOXTION  COEriCtENT  is  .23*1M3SE>*] 

rx  Mcr**crivE  inoti  minus  unity  is  -.iTsoj**te-*3 


rx  »u«*  rucouENCirs  uk  •aaasooaoE.a*  .2«isasasE>a* 

rx  po*ulmtion  TNawirEXo  setxen  iojaccnt  level*  is.  .  u«oooooe-o» 

ro*  rx  r*coufXY  .»a**«om*«» 

rx  *oso**'ion  coErici£*r  is  .saiaaaax-a* 

rx  xr»»erive  inoti  nUus  unity  is  jaat,aaiE*a» 

ro*  rx  r»Eou€XY  .jaisaaoOE*** 

rx  *050**1  ion  coettcient  is  .2»in«*aE-«3 

rx  retmctive  inoti  minus  unity  is  -.irssaXOt-*} 


rx  *un*  r*EouEXirs  »«t  .***as**eE>*»  ,2o2ooooo£* 

rx  *0*UL»riON  r*MSTEMO  SeTXEN  iojaceni  level*  is.  -I 

rx  rx  rxouEXY  .oasaaosOE*** 

rx  aosooTi ix  coeticient  is  .soiaoaax-a* 

rx  *ET**CTIvE  INOTI  minus  unity  Is  -.)H«*SI1E-1* 

rx  rx  roEoxxY  .2a2oos**E*a* 

rx  *o*o»*r ion  cxrtctENr  is  .«tt3*tie-*3 

rx  x**»eri»E  inoti  minus  unity  is  i umiii-i! 
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thc  pump  mtuucus  u*  .uhihi(>h 

ThC  POPulaT  ION  TRansTCRCO  KTmCCn  aOjaCCnT  Llvb.1  is.  .lHIIHICtt 

roa  thc  rscaucNCT  .»Moiii<c<h 

THC  ABSORPTION  COCTICILNT  It  .IIIHMC'h 

THC  RCTRaCTMC  [HOC I  HluuS  UNITY  IS  -,]AAA»H IC-A* 

roa  thc  racoucNct  .2ACsso«oc>«* 
thc  AasoaaTioH  cocttcict  is 

THC  RCTRACTIVC  1N0C>  Minus  um|(T  IS  I I1H.W-U 


thc  w*  racouCHCtcs  ahc  .<«iihii('<4  ,2tsiiiasC*0A 

THC  POPULATION  TAAHSreaCD  BCTmaCN  AOJACCNT  LCVCLS  IS.  .IBOOBOOOC-OA 

roa  thC  racouCNcr  .ttitiMH*)* 

thc  absorption  cocticicht  is  .sbiabbbbc-aa 

thc  rctractivc  ihoti  huus  unity  is  ..jMteiitc-Tt 

roa  thc  racouCNcr  .cusoomac.a* 

thc  AasoaaTioN  cocticUht  is  .IVtlkMI-il 

THC  RCraaCTIvC  INOCI  MINUS  unity  IS  1SA12ASAC-AJ 


THC  PUN*  rRCOUCNCirS  ARC  .HMIMt(’). 

thc  population  transtcRCo  hthcCn  aOjaccnt  lcvcls  is.  .Iliggiii(-t4 

roa  thc  racoiiCigcY  .a*b«a0aAC.*» 

THC  ABSOHPT ION  COCrteiCNr  IS  -UllUtW-IY 
thc  acratcTivc  IHOTI  Minus  unity  is  •.]SMa«oic*4a 

roa  thc  racouCNCY  .cotsoooac-o* 

THC  ABSORPTION  COCPTCICnT  IS  .2ST3lAS*C-iJ 
thc  rctractivc  ihoti  nIhus  unity  is  •.iotassaic**! 


THC  PUMP  rHCOUCNCICS  ahc  .AAAAAAAAC.Ca  .20*OOAA«C>BA 

THC  POPULATION  TNAHsrCMCo  BCTM.CN  aOJACCMT  -LCVCLS  is.  .Ioasaabbc-al 

roa  thc  racouCNcr  .tHggfggc-g. 

THC  ABSORPTION  COCTICICHT  is  .SBlaSBBBC-A* 

THC  RCTBACTIVC  INOTI  MINUS  unity  is  -.SBAAaASIC'AA 

roa  Thc  racouCNCY  .2o»oo#«ac-y» 

THC  ABSORPTION  COCriCItHT  Is  .23T2iaBlC-4) 

thc  rctractivc  ihoti  hihus  unity  is  •.iaaisboic'A) 


thc  pump  racouCNCirs  arc  .AAAMoeoC-4*  .igisoiigl-gt 

THC  POPULATION  TRANSTECH  BCThCCh  aOJACCMT  LCVCLS  is.  .1001 

roa  thc  racouCNCY  .aosaoaaac.oa 

THC  ABSORPTION  COCTICICNT  Is  .SBISS*«AC-«A 

*  "■  TAINP.  Ml.  .*■  .I.|f»  .  lA*Ar.«rtte«*A 
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yon  roc  ntcauCNcr  ,^*hhk>u 

tmc  mmwmoii  eocrictCMT  is  .*stisiwc-sj 

tmC  DC'NACTIvC  INOTi  Ml»uS  UNIT*  IS  •  .  lt*?*CMC**3 


tmc  w*  nnwNcin  »»t  .miihiM'U  .csmscmc*** 

tmc  aonulaiion  T**««rc*co  k'«u>  ojaccnt  li«li  is.  .imhiik-m 

TON  TMC  FMCOUCMCY  .illalllK-l. 

TMC  ACSONNMON  COCrtClCOT  Is  .SSItMMC*** 

Tmc  NCfNACTIVC  INOYI  Mlwus  UNITY  IS  I 

TON  TmC  TKOUCMCT 

roc  *mon#tion  cocriciCMT  is 

TMC  «Cr»CCTlvC  INOrt  MINUS  unity  is  -.ICUCIMC"*) 


THATAS  ALL 
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